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FOREWORD 


In  this  report  we  will  discuss  the  progress  that  has  been  made 
in  the  area  of  understanding  and  quantifying  noise  spectra  that 
have  a  component  with  a  characteristic  frequency  dependence 
proportional  to  1/f.  Fluctuations  in  the  time  domain  of  the  voltage 
across  a  semiconductor  device,  or  of  the  current  through  it,  are 
ever  present.  With  the  down  scaling  of  the  size  of  electronic 
devices  the  signals  are  becoming  smaller  as  well.  Hence,  one  of  the 
limiting  factors  determining  the  device  performance,  the  noise  the 
device  itself  produces,  becomes  important.  This  becomes  apparent 
first  when  we  consider  detectors  which  we  use  to  distinguish  useful 
signals  from  "background"  signals. 

Understanding  the  noise  produced  by  a  device  is  to  be 
understood  as  (i)  being  able  to  identify  the  physical  mechanism 
that  gives  rise  to  the  fluctuations  and  (ii)  to  know  how  a  locally 
generated  fluctuation,  in  let  us  say  the  current,  will  contribute 
to  the  noise  as  it  is  measured  at  the  terminals  of  the  device. 
Whereas  the  first  has  to  do  with  the  basic  physical  forces  that 
work  on  a  carrier  as  it  traverses  a  device  the  second  has  to  do 
with  the  operating  principles  of  the  device  which  describe  the 
behavior  of  average  quantities  such  as  constant  current,  etc. 
Another  way  of  putting  it  is  considering  the  basic  carrier 
interactions  on  a  microscopic  level  and  evaluate  how  these 
processes  contribute  to  the  observations  on  macroscopic  level. 

Ideally  one  would  like  to  have  a  description  of  the  noise 
produced  by  the  device  as  a  function  of  the  operating  conditions. 
This  would  allow  one  to  predict  the  noisiness  of,  for  instance,  a 
transistor  in  a  particular  circuit  application.  One  idea  in  this 
direction  was  formulated  first  in  1975  by  P.  H.  Handel.  Since  1984 
A.  van  der  Ziel  has  been  taking  noise  measurements  on  devices  to 
either  support  or  refute  this  theory.  Some  of  the  experimental 
evidence  seems  to  support  the  theory.  In  the  last  two  years,  1989 
and  1990,  more  experimental  data  was  gathered  and  some  of  the 
aspects  of  the  theory  have  been  considered.  The  outcome  of  these 
investigations  is  collected  and  presented  in  this  report. 

Under  the  grant  the  principal  investigators  were  A.  van  der 
Ziel  (deceased,  January  1991)  and  A.  D.  van  Rheenen,  both  at  the 
University  of  Minnesota.  The  theory  part  was  subcontracted  with  P. 
H.  Handel,  University  of  Missouri-St . Louis . 
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I.  EXPERIMENTAL  RESULTS 


1 .  Introduction 


In  the  early  history  of  che  study  of  noise  in  electronic 
devices  investigators  were  particularly  intrigued  by  the  1/f 
spectrum  and  several  theories  were  developed  to  explain  the  unique 
spectral  shape.  Some  of  the  more  obvious  difficulties  with  a  pure 
1/f-shaped  power  spectral  density  is  that  integration  of  the 
spectrum  over  all  frequencies  yields  infinities  at  both  integration 
limits.  A  realistic  physical  system  is  of  course  properly  behaved, 
with  none  of  these  infinities.  In  practice  the  noise  measurements 
are  band  limited:  on  the  high-frequency  side  of  the  spectrum  by  the 
frequency  response  of  the  measuring  system  and  on  the  low-frequency 
side  by  the  time  the  equipment  is  on. 

One  attempt  to  describe  the  1/f  noise  was  to  assume  that  a 
number  of  traps  with  a  distribution  of  trapping  time  constants  is 
responsible  for  the  1/f  shape.  A  trapping  center  that  captures  and 
releases  again  a  charge  carrier  gives  rise  to  a  Lorentzian  shaped 
(1/[1  +  )  contribution  to  the  noise  spectrum.  It  is  possible 

to  simulate  1/f  behavior  of  the  spectrum  over  decades  of  frequency 
by  adding  Lorentzians  with  proper  time  constants  x.  Since  the 
theory  for  generation-recombination  (g-r)  noise  (the  noise  related 
to  trapping  and  de-trapping  of  carriers)  is  well  known  this 
approach  would  give  a  physical  interpretation.  It  also  eliminates 
the  mathematical  problems  mentioned  earlier.  A  sum  of  Lorentzians 
will  result  in  a  frequency  independent  part  at  low  frequencies  and 
a  part  that  rolls  off  as  1/f^  at  high  frequencies.  The  intermediate 
region  exhibits  the  1/f  spectrum.  Integration  of  this  spectrum  over 
the  frequency  gives  a  finite  result.  These  attributes  give  this 
theory  its  appeal.  The  only  issue  would  be  to  identify  traps  with 
exactly  the  right  distribution  of  time  constants,  which  by  the  way 
can  be  related  to  a  distribution  of  activation  energy  levels.  The 
first  person  to  propose  such  a  model  was  McWorther  [1]. 

In  1969  it  was  observed  by  Hooge  [2]  that  the  essence  of  1/f 
noise  in  a  wide  variety  of  devices  could  be  captured  by  a  single 
equation  relating  the  power  spectral  density  (SJ  as  a  function  of 
frequency  (f)  to  the  DC  operating  current  (I)  and  the  number  of 
free  carriers  (N)  in  the  device.  It  was  found  that  the 
proportionality  constant  (now  named  after  Hooge)  Oh  had  a  constant 
value  of  2xl0*\  independent  of  the  material  or  size  of  the  device. 
This  relation  is  expressed  as  follows: 

a„  I' 

s,(f)  - - .  (1) 

fN 

The  advantage  of  the  formulation  of  Eq.  1  is  that  by  assuming  that 
the  Hooge  parameter  is  a  constant  one  can  predict  the  noise  the 
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device  will  produce  under  certain  operating  conditions. 

With  the  advances  made  in  the  technology,  smaller  and  smaller 
devices  were  fabricated  and  from  noise  measurements  smaller  and 
smaller  values  for  the  Hooge  parameter  were  deduced.  Possibly  the 
improved  processing  techniques  resulted  in  cleaner  samples 
resulting  in  less  noisy  devices-  With  the  observation  that  the 
Hooge  parameter  was  no  longer  a  constant  the  interpretation  of  Eq. 
1  and  of  the  Hooge  parameter  changed.  Now  Eq.  1  would  be  used  to 
normalize  the  measured  current  noise  spectrum  to  find  the  parameter 
Oh-  This  then  allows  one  to  compare  the  noise  behavior  of  different 
devices  with  different  carrier  numbers  under  different  operating 
conditions.  The  Hooge  parameter  would  simply  be  a  measure  of  the 
noisiness  of  the  device;  a  small  parameter  value  would  indicate  a 
quiet  device  and  a  large  value  would  show  that  the  device  produces 
a  lot  of  noise. 

It  was  also  observed  that  the  value  of  a„  depended  on  the 
electric  field,  see  Bosman  et .  al.  [33.  This  made  the  case  stronger 
to  let  the  Hooge  parameter  be  an  adjustable  parameter  in  stead  of 
a  constant.  Continuing  this  line  of  thought  it  became  apparent  that 
there  was  a  need  for  a  theory  describing  how  the  noisiness 
parameter  depends  on  external  parameters  such  as  the  temperature, 
the  electric  field,  etc.  An  attempt  to  describe  the  noisiness  from 
first  principles  was  made  by  Handel  in  1975  [4]  .  The  physical 
mechanism  invoked  is  the  emission  of  Bremsstrahlung  when  electrons 
are  accelerated  and  it  was  shown  that  this  would  result  in  a  1/f^'® 
shaped  spectrum  with  e  a  positive  number,  small  compared  to  1.  This 
formulation  does  remove  the  infinity  at  f  =  0  when  integrating  the 
spectrum.  The  magnitude  of  the  noise  is  expressed  in  terms  of  the 
velocity  changes  the  electron  undergoes  when  it  interacts  with 
phonons,  impurities,  other  electrons,  etc.  This  formulation  then 
would  enable  one  to  predict  the  noisiness  of  the  device  from 
knowledge  of  the  dominant  scattering  mechanism  under  certain  device 
operating  conditions.  In  principle  the  theory  would  be  a  powerful 
tool  to  study  and  predict  the  magnitude  of  the  1/f  noise.  However, 
this  theory  has  been  severely  criticized,  both  on  the  grounds  of 
the  mathematical  detail  and  the  underlying  physical  assumptions. 

In  order  to  shed  some  light  on  the  controversy  the  Noise  group 
at  the  University  of  Minnesota  under  the  guidance  of  the  late  Dr, 
A.  van  der  Ziel  set  out  to  measure  noise  in  a  variety  of  devices 
and  compare  the  resulting  values  of  the  noisiness  (still  called 
Hooge)  parameter  with  predictions  made  by  Handel's  theory.  The  idea 
was  to  either  confirm  or  refute  the  theory  based  on  these 
measurements.  The  body  of  this  report  is  concerned  with  the  results 
that  were  obtained  over  the  period  of  the  grant  from  January  1, 
1989  through  March  31,  1991.  A  large  number  of  different  devices 
were  studied:  n*-p  and  n-i-p  Hgj.^Cd^Te  infrared  detector  diodes, 
AlGaAs/GaAs  laser  diodes,  InGaAs/InP  p-i-n  photodiodes,  silicon 
bipolar  transistors,  AlGaAs/GaAs  hetero junction  bipolar 
transistors,  silicon  Junction  Field  Effect  Transistors,  and 
AlGaAs/GaAs  Double  Barrier  Resonant  Tunneling  Diodes.  Most  of  the 
results  and  interpretation  of  the  data  has  been  published  in  the 
open  literature.  What  follows  is  a  recapture  of  these  results  in 
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the  form  of  a  section  on  each  subject.  At  the  end  the 
accomplishments  on  the  theory  part  of  the  project  will  be  high 
lighted. 


2 .  Results  of  research  as  published. 

2.1  Extensions  of  Handel's  1//  noise  equations  and  their 

semiclassical  theory  (see  p.  B2  for  full  paper) 

Handel's  expression  for  the  Hooge  parameter  (Qh)  which 

quantifies  the  1/f  noise  can  be  modified,  allowing  for  bunches  of 
charge  qe,  to  the  following  expression: 

4a„  q  F,(0) 

(Xh  = -  (-)2  ( - )2.  (2) 

37C  e  c 

In  this  expression  is  the  fine-structure  constant  (1/137),  c  the 
speed  of  light  and  Fa(0)  the  low-frequency  Fourier  transform  of  the 
acceleration,  a(t),  of  the  carrier.  The  differences  with  Handel's 
original  expressions  are  that  now  bunches  of  charge,  as  they  appear 
in  vacuum  devices,  can  be  accounted  for  and  second,  that  the 
vectorial  velocity  change  that  appears  in  Handel's  equation  can  be 
derived  from  semi-classical  (no  quantum  mechanics)  arguments  as  the 
Fourier  transform  of  the  acceleration. 


2.2  1/f  noise  characterization  of  n*-p  and  n-i-p  Hgi.^Cd^Te 

detectors  (see  p.  B6  for  full  paper) 

In  this  experimental  study  the  noise  behavior  of  Hgo.7Cdo.3Te  n^- 
p  diodes  supplied  by  Rockwell  International  was  investigated  in  the 
temperature  range  from  80-300  K  under  small  reverse  bias 
conditions.  The  theoretical  values  for  Oh  were  estimated  as 
follows : 

(i)  Coherent  state  1/f  noise  a„  =  4.6x10"^, 

(ii)  Umklapp  1/f  noise  a„  =  5x10"^  for  T  =  273  K, 

(iii)  Normal  collision  1/f  noise  a„  =  1x10'^  for  T  =  77  K. 

In  this  bias  regime  the  noise  spectral  density  of  the  current 
fluctuations  is  given  by: 

el 

Si(f)  =  a„ - F(V)  .  (3) 

ft 

Here,  F (V)  is  a  factor  dependent  on  the  bias  voltage  V  that  comes 
about  when  one  integrates  the  locally  generated  noise  along  the 
device  length  to  find  the  noise  at  the  terminals.  It  has  been  shown 
earlier  that  F  (V)  has  the  following  form  for  a  reverse  biased 
diode : 
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F(V)  —  +  —  -  —  ln(l+a)  (4) 

3  2a  a^  a^ 

where  a  =  exp(qV/kT)  -  1. 

The  time  constant  t  was  extracted  from  measurements  of  the  complex 
admittance  of  the  diode  and  fitted  to: 

y(0))  =  g^d+jcot)'*.  (5) 

This  yielded  values  for  the  lifetime  that  were  commensurate  with 
the  values  of  the  Honeywell  lifetime  tables. 

The  experiments  resulted  in  the  following  observations: 

(i)  in  the  temperature  range  from  113  -  193  K  the  values  for  the 
Hooge  parameter  (Ch)  were  (3-5)xl0“^  indicating  the  presence  of  the 
coherent  state  noise. 

(ii)  around  T  =  273  K  one  device  exhibited  a  value  of  5x10“^  as  one 
would  expect  for  the  Umklapp  noise. 

(iii)  normal  collision  1//  noise  has  not  been  observed.  This  noise 
contribution  is  probably  masked  by  the  other  two  sources. 

The  second  part  of  this  study  concerned  the  noise  performance 
of  Hgo,78Cdo.22Te  p-i-n  diodes  supplied  by  Santa  Barbara  Research 
Center.  These  diodes  were  tested  at  80  K.  In  an  earlier 
contribution  [5]  measurements  of  the  admittance  of  a  silicon  power 
diode  were  used  to  obtain  the  lifetime  of  the  carriers.  This 
lifetime  is  essential  when  one  wants  to  determine  the  Hooge 
parameter  from  Eg.  3.  This  expression  follows  from  substituting 
qN/x  for  one  of  the  factors  I,  into  Eq.  1.  This  analysis  resulted 
in  lifetimes  of  the  order  of  lO"’  s,  comparable  to  the  Honeywell 
lifetime  tables.  The  value  of  a„  was  found  to  be  4x10"^,  and  this 
value  agrees  with  the  predictions  made  according  to  the  coherent 
state  1/f  noise,  2a/K,  which  is  equal  to  4.6x10'^. 


2.3  Secondary  emission  1/f  noise  revisited  (see  p.  Bll  for  full 

paper) 

A  case  is  made  here  that  noise  measurements,  taken  in  1960,  of 
the  1/f  noise  in  secondary  emission  pentodes  conform  to  predictions 
made  by  Handel's  1/f  noise  theory.  The  analysis  is  carried  out  with 
one  adjustable  parameter,  the  distance  between  the  dynode  and  the 
anode,  and  assuming  Handel's  predictions  are  correct.  Two  sets  of 
measurements  taken  on  each  of  two  devices:  one  set  involved  a 
constant  anode  current  and  varying  anode-to-dynode  voltage,  and  the 
other  set  involved  a  constant  anode-to-dynode  voltage  and  varying 
anode  current.  From  a  fit  of  the  1/f  noise  measurements  the  dynode- 
to-anode  distance  was  found  to  be  4.8  ±  0.3  mm  for  one  device  and 
7.0  ±  0.6  mm  for  the  other.  The  physical  distance  is  estimated  at 
6  mm.  In  light  of  the  fact  that  the  dynode  and  anode  are  not 
parallel  pJ ates  but  curved  surfaces  one  can  argue  that  the  electron 
beam  might  traverse  a  path  length  that  differs  slightly  from  the 
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estimated  physical  distance  between  the  electrodes.  It  seems 
reasonable  to  assume  that  the  values  deduced  from  the  noise 
measurements  agree  with  the  physical  dimension.  This  then  supports 
Handel's  prediction  for  the  magnitude  of  the  1/f  noise  in  a  vacuum 
device  where  electrons  do  not  suffer  collisions  while  in  transit. 


2.4  Generation-recombination-type  l/f  noise  in  n-i-p  diodes  (see 
p.  B14  for  full  paper) 

For  devices,  such  as  resistors,  that  have  homogeneous  carrier 
numbers  (N)  and  current  densities  (I)  it  is  postulated  that  Hf 
contribution  to  the  spectral  intensity  of  the  current  fluctuations 
(Si)  can  be  described  by: 

a„  I' 

Si  - - .  (6) 

f  N 

In  the  case  that  the  device  is  not  homogeneous  a  different  approach 
must  be  used.  The  procedure  is  as  follows:  the  device  is  divided  in 
slices  of  thickness  Ax  and  the  contribution  (ASi)  of  that  slice  to 
the  total  noise  is  calculated.  In  particular,  it  is  assumed  here 
that  the  noise  is  associated  with  the  generation-  and/or 
recombination  (g-r)  processes  of  carriers  interacting  with  traps. 
The  development  of  the  model  then  proceeds  by  replacing  I  with 


(dI/dx)Ax,  N  with  (dNT./dx)Ax,  and  realizing  that  dl  =  q  dN/x,  so 
that  we  find: 

q  ttH  dl 

ASi  - - -  Ax.  (7) 

ft  dx 

Upon  integration  of  Eq.  7  we  find: 
q  ttH  I 

ASi  - - .  (8) 

f  X 


This  expression  was  used  to  calculate  the  Hooge  parameter  from 
noise  measurements  taken  on  Hgo.8Cdo.2Te  p-i-n  diodes  at  liquid 
nitrogen  temperatures.  This  analysis  resulted  in  values  varying 
from  2x10"®  to  2x10'^. 

The  reason  this  noise  was  interpreted  as  g-r-type  1/f 
noise  follows  from  the  observation  that  the  frequency  dependence  of 
the  best  fit  to  the  measured  spectra  behaves  as  1/f^,  with  y  larger 
than  1.  Therefore  the  values  we  give  for  the  Hooge  parameter  need 
to  be  viewed  as  upper  limits;  the  1/f  noise  is  certainly  not  larger 
than  the  values  given.  As  we  reported  earlier  (Sect.  2.2)  it  is 
expected  that  Umklapp  noise  gives  rise  to  a  value  of  10"^  for  the 
Hooge  parameter.  We  find  smaller  values  here.  Collision  1/f  noise 
would  result  in  values  of  10'®,  two  orders  of  magnitude  smaller  than 
the  upper  limits  we  found.  Possibly,  this  collision  1/f  noise 
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contribution  is  masked  by  the  g-r  noise. 


2.5  Burst-type  noise  mechanisms  in  bipolar  transistors  (see  p. 

B17  for  full  paper) 

When  one  observes  the  noise  signal  in  the  time  domain,  let  us 
say  by  using  an  oscilloscope,  it  shows  an  average  level  on  which 
both  upward  and  downward  spikes  are  superimposed.  The  occurrence  of 
the  spikes  as  well  as  their  amplitude  are  "random".  In  some  devices 
deviations  from  this  time  domain  signal  have  been  observed.  Pulses 
of  varying  duration  seem  to  be  superimposed  on  the  ordinary  noise 
signal,  raising  the  average  locally,  and  then  dropping  to  the 
"normal"  average.  The  occurrence  of  the  pulses  is  "random"  in  time 
but  their  amplitude  is  fixed.  This  type  of  noise  is  called  burst 
noise.  One  can  hear  a  distinct  popping  sound  when  one  feeds  this 
amplified  signal  to  a  loud  speaker.  This  explains  why  this  kind  of 
noise  is  sometimes  referred  to  as  popcorn  noise. 

The  occurrence  of  the  pulse  in  the  time  domain  translates  to 
a  Lorentzian,  1/ (l+co^t’),  shaped  contribution  to  the  noise  spectrum, 
in  addition  to  the  usual  1/f  noise.  Here,  measurements  of  the 
collector  current  noise  are  taken  on  two  different  kinds  of 
devices,  a  discrete  bipolar  junction  transistor  (BJT)  and  a  BJT  on 
a  chip.  Two  distinct  noise  spectra  were  observed,  best  described  by 
the  theoretical  fits  to  the  spectra: 

At^  +  B/f 

discrete  BJT  Sj.  - ,  (9) 

1  +  co^t^ 

AX^ 

chip  BJT  Si  = - +  B/f.  (10) 

1  +  CO^T^ 

By  inspecting  Eq.  9  we  see  that  for  large  frequencies  S.^l/f^  for 
the  discrete  device  whereas  Eq.  10  indicates  at  large 

frequencies  for  the  chip  BJT,  At  the  low-frequency  side  of  the 
spectra  both  equations  predict  identical  behavior. 

The  observation  of  two  distinctly  different  burst  noise 
spectra,  mixed  with  1/f  noise,  as  measured  in  two  different  devices 
is  explained  by  assuming  either  a  modulation  of  the  1/f  noise  by 
the  burst  noise  (Eq.  9)  indicating  a  common  source  for  the 
fluctuations  or  a  superposition  of  the  two  noise  sources  indicating 
unrelated  origins  for  the  fluctuations. 


2.6  1/f  noise  in  double-hetero junction  AlGaAs/GaAs  laser  diodes  on 

GaAs  and  Si  substrates  (see  p.  B21  for  full  paper) 


In  this  study  a  comparison  is  mav^e  between  the  noise  in  the 
electrical  drive  current  of  double-hetero junction  AlGaAs/GaAs 
lasers  built  on  a  silicon  substrate  and  ones  built  on  a  GaAs 
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substrate.  It  might  be  advantageous,  for  purposes  of  easy 
integration,  to  fabricate  lasers  for  optical  communications  on  the 
silicon.  Of  course  the  materials  issues  of  growing  GaAs  on  Si  need 
to  be  addressed,  but  also  the  electrical  and  optical  performance  of 
the  devices  need  to  be  characterized  and  possibly  explained. 

Often  the  laser  behavior  is  modelled  in  terms  of  rate 
equations  describing  the  evolution  in  time  of  the  number  of 
carriers  and  the  number  of  photons.  These  rate  equations  are 
coupled  and  it  is  expected  that  fluctuations  in  the  drive  current 
will  give  rise  to  fluctuations  in  the  optical  output  of  the  laser. 
Yamamoto  [6]  speculated  that  the  spectrum  of  the  frequency 
deviations  (phase  noise)  of  the  laser  output  is  limited  at  small 
frequencies  by  the  fluctuations  in  the  electrical  current. 

Measurements  of  the  low-frequency  noise  in  the  drive  current 
of  the  laser  diodes  as  a  function  of  the  drive  are  the  result  of 
this  investigation.  Both  types  of  devices  showed  the  magnitTide  of 
the  noise  divided  by  the  current  to  be  independent  of  the  current 
for  currents  smaller  than  0.1  mA  and  to  decrease  for  larger 
currents.  The  quantity  S^f/I  is  proportional  to  0.^/1,  see  Eq.  8.  If 
anything,  one  would  expect  the  lifetime  of  the  carriers  to  be 
reduced  when  more  carriers  are  injected  as  the  current  increases, 
due  to  an  increased  recombination.  This  would  increase  the  ratio 
o.^/X.  Therefore  we  must  conclude  that  Oh  decreases  as  the  current 
increases.  A  similar  effect  was  observed  in  Si  n*-v-n*  and  p'’-71-p* 
structures  [7] . 

The  two  curves,  see  figure  4  on  p.  B23,  of  S^f/I  versus  the 
current  for  the  two  types  of  devices  run  parallel,  with  the  noise 
of  the  devices  fabricated  on  the  Si  substrate  about  50  times  as 
noisy  as  the  ones  on  the  GaAs  substrate.  Because  the  current  flows 
through  the  substrate  it  i.s  possible  that  the  misfit  dislocations 
due  to  the  lattice  mismatch  at  the  Si-GaAs  substrate,  give  rise  to 
excess  noise.  This  would  explain  the  noisier  behavior  of  the  lasers 
on  the  Si  substrate. 


2.7  Low-frequency  noise  in  small  InGaAs/InP  p-i-n  diodes  under 
different  bias  and  illumination  conditions  (see  p.  B25  for 
ful]  paper) 

With  the  advances  made  in  the  1.3-1.55  )im  wavelength 
communications  it  is  natural  that  high-speed  photodetectors  are 
considered  to  allow  for  high  data  rates.  The  InGaAs/InP  p-i-n  diode 
is  such  a  fast  detector.  It  is  able  to  detect  signals  up  to  a  few 
Gbits/s.  It  is  important  to  characterize  the  noise  performance  of 
the  detectors  as  signal  levels  become  smaller. 

Here  the  noise  is  stuoied  in  devices  that  are  reverse  biased 
and  illuminated  and  the  results  are  compared  with  noise  data  talcen 
on  devices  that  are  forward  biased  and  not  illuminated.  The  most 
remarkable  result  that  the  frequency  exponent  of  the  spectra  of 
four  devices  was  the  same  under  the  two  operating  conditions  we 
just  mentioned,  although  it  differed  from  device  to  device.  Two 
devices  exhibited  a  slope  of  -1  in  the  noise  spectrum  whereas  the 
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other  two  showed  a  slope  of  -0.8.  The  fact  that  the  spectral  shape 
is  the  same  when  the  devices  are  forward  biased  and  reverse  biased 
tends  to  indicate  that  the  mechanisms  responsible  for  the  noise  are 
present  under  both  bias  conditions.  Second,  the  extrapolated  value 
of  Sid  Hz) /I^,  averaged  over  at  least  five  current  values,  for  the 
devices  under  reverse  bias  and  illumination  was  found  to  be  0.7x10' 
“  and  2.6x10'“  /Hz  for  two  devices  (spectral  shape  f'“°)  and  about 
9.5x10'“  /Hz  for  the  devices  with  a  spectral  shape  of  f'°-®.  See 
table  I  on  p.  B26.  Although  there  is  some  difference  between  the 
noise  magnitudes  of  these  two  sets  of  devices,  they  are  within  one 
order  of  magnitude.  This  again  suggests  that  the  same  mechanisms 
govern  the  noise  production.  Third,  under  forward  bias,  no 
illumination  the  value  of  S^l  Hz) /I,  again  averaged  over  at  least 
five  current  values,  is  virtually  the  same  for  all  devices:  5x10'“ 
A/Hz.  This  in  spite  of  the  differences  in  the  shape  of  the  noise 
spectra  for  the  different  devices.  See  table  II  on  p.  B28. 

It  is  speculated  that  fluctuations  in  the  quantum  efficiency 
of  the  detector  can  result  in  the  I^  dependence  of  the  noise  when 
the  diodes  are  reverse  biased  and  illuminated.  This  is  the  subject 
of  a  following  paper.  Section  2.9. 


2.8  Extension  of  the  Hooge  equation  and  the  Hooge  parameter 
concept  (see  p.  B30  for  full  paper) 

The  Hooge  equation  should  be  applied  strictly  to  noise  spectra 
that  are  proportional  to  f'“°.  However,  many  spectra  exhibit  a 
frequency  exponent  y  that  deviates  from  -1.0.  It  is  proposed  here 
to  modify  the  Hooge  equation  to  read: 

««  I' 

Si(f)  - -  .  (11) 

N  fT 

When  we  allow  this  a^,  is  not  a  dimensionless  parameter  anymore,  but 
has  units  Hz'’^'^.  This  formulation  would  allow  the  comparison  of  the 
noise  parameter  of  different  devices  by  comparing  the  value  of  the 
spectral  intensity  at  1  Hz.  Since  we  now  allow  for  the  spectral 
shape  to  be  different  that  comparison  is  not  easily  extended  to 
different  frequencies. 


2.9  Noise  in  quantum  efficiency  of  p*-n  diodes  due  to  fluctuations 
in  the  surface  recombination  of  carriers  (see  p.  B32  for  full 
paper) 

In  a  brief  note  it  is  shown  that  spectral  intensity  of 
fluctuations  in  the  quantum  efficiency  of  the  p*-n  diodes  due  to 
fluctuations  in  the  surface  recombination  of  the  carriers  scales 
with  the  quantum  efficiency  squared.  Since  the  current  due  to 
incident  light  on  the  detector  is  proportional  to  the  quantum 
efficiency  the  noise  in  the  current  will  be  proportinal  to  the 
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current  squared. 


2.10  Current  fluctuations  in  double-barrier  quantum  well  resonant 

tunneling  diodes  (see  p.  B33  for  full  paper) 

Resonant  tunneling  diodes  are  some  of  the  devices  with 
possible  operating  frequencies  close  to  1  THz .  Applications  that 
make  use  of  these  devices  are  for  instance  mixers.  The  noise 
performance  of  the  diodes  in  such  applications  is  important  because 
the  low-frequency  noise  will  appear  as  phase  noise  in  the  mixed 
signal.  In  this  study  measurements  of  the  noise  taken  on  a  device 
with  AlAs  barriers  are  compared  with  data  taken  on  identical 
devices  that  have  Alo.3Gao.7As  barriers  in  stead  of  the  AlAs 
barriers.  The  measured  spectra  are  made  up  of  a  1/f  noise 
contribution  and  contributions  related  to  trapping  of  carriers,  g-r 
noise . 

In  both  devices  the  magnitude  of  the  1/f  noise  did  not  change 
with  temperature  in  the  temperature  range  from  78-400  K  when  the 
diodes  were  biased  close  to  the  peak.  Since  the  tunneling  component 
of  the  current  is  independent  of  the  temperature  it  seems  to 
suggest  that  the  1/f  noise  is  related  to  the  tunneling  current. 
Also,  the  magnitude  of  the  1/f  noise  is  the  same  for  both  types  of 
devices  in  spite  of  the  different  barrier  materials.  When  the 
diodes  are  biased  at  voltages  larger  than  the  valley  voltage  the 
devices  are  less  noisy  at  the  same  operating  current. 

From  the  g-r  noise  measurements  the  activation  energies  of  the 
traps  were  extracted.  We  observed  five  traps  in  the  AlAs  barrier 
device  and  two  in  the  AlGaAs  barrier  device  with  activation 
energies  equal  to  two  of  the  activation  energies  of  the  AlAs 
device.  This  seems  to  suggest  that  increasing  the  A1  mole  fraction 
of  the  barrier  material  introduces  more  traps 


3 .  Results  not  vet  published. 


3.1  Noise  in  BJTs  and  HBTs  (Ph.  D.  Thesis  Alister  Young) 

(for  Abstract  see  p.  C2) 

It  was  proposed  by  van  der  Ziel  that  it  would  be  possible  to 
identify  the  physical  location  of  the  sources  responsible  for  the 
noise  in  bipolar  transistors  by  measuring  the  noise  in  the 
collector  current  and  the  base  current  when  the  transistor  is 
biased  in  the  common  emitter  configuration  and  when  the  transistor 
bias  configuration  includes  a  resistor  between  the  emitter  and 
ground.  Based  on  a  small  signal  equivalent  circuit  model  that 
includes  possible  independent  noise  sources  the  analysis  proceeds 
by  fitting  the  magnitude  of  the  noise  sources  as  well  as  the  base 
spreading  resistance  to  the  four  noise  measurement  data.  This 
procedure,  although  developed  earlier,  has  never  been  applied  to 
study  the  noise  behavior  of  actual  devices.  In  this  thesis,  the 
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first  experimental  test  of  the  procedure  is  presented.  Noise 
measurements  were  taken  on  commercial  high-frequency  silicon 
bipolar  transistors  (Motorola  MRF941  and  Motorola  MRF951)  as  well 
as  on  experimental  AlGaAs/GaAs  hetero junction  bipolar  transistors 
(HBT)  fabricated  by  Rockwell  International  and  by  Tektronix. 

Silicon  Bipolar  Transistor  (MRF941  and  MRF951) 

The  small  signal  equivalent  circuit  model  includes  three  noise 
sources:  (i)  a  source  between  base  and  collector,  ibc/  due  to  the 
generation  in  the  base-collector  spce  chaz-ge  region,  (ii)  a  source 
between  the  base  and  the  emitter,  i^*,  resulting  from  recombination 
in  the  base,  and  (iii)  a  source  between  the  collector  and  emitter, 
i^af  associated  with  the  diffusion  of  carriers  from  the  emitter  to 
the  collector.  By  varying  the  resistance  in  the  emitter  lead  we 
manipulate  the  relative  contributions  of  these  sources  to  the 
terminal  noise  current.  The  following  observations  are  made: 

a.  The  magnitude  of  i^c  is  much  smaller  than  that  of  ibe-  This  can 
be  expected  from  improvementsthat  have  been  made  in  the  processing 
of  the  materials. 

b.  The  collector  current  noise  is  dominated  by  the  contribution  of 
source  ice* 

c.  Depending  on  the  bias  current  values  between  10  and  100  Cl  were 
found  for  the  base  spreading  resistance,  reasonable  for  the  device 
structure  and  doping  levels. 

d.  The  values  of  the  Hooge  parameter  were  estimated  as  follows: 


MRF941 

base  current: 

0.1 

-  55  jXA 

ttn: 

6x10-’ 

-  5x10-' 

MRF951 

base  current: 

0.01 

-  15  ^A 

CCm: 

7x10-® 

-  IxlO-** 

MRF941 

collector  current: 

0.1  -  1650 

\IA 

cCh  : 

2x10-®  - 

2x10-® 

MRF951 

collector  current: 

0.1  -  1250 

\IA 

a„: 

2x10-®  - 

2x10-® 

In  general  the  value  of  ttw  seems  to  increase  with  the  current 
and  the  MRF951  appears  to  be  quieter  than  the  MRF941.  At  the  moment 
more  careful  analysis  is  under  way.  We  expect  to  publish  the 
results  of  this  more  detailed  study  shortly. 

Hetero junction  Bipolar  Transistors  (AlGaAs/GaAs  from  Rockwell  and 
Tektronix) 

Both  types  of  devices  are  considered  experimantal  and  no 
detailed  information  is  available  to  us  concerning  the  precise 
device  structure.  Both  are  vertical  n-p-n  transistors  and  have  a 
cut-off  frequency  of  200  GHz.  The  maximum  value  of  hp-E  is  2.5  for 
the  Rockwell  device  and  60  for  the  Tektronix  device. 

The  base  current  noise  spectra  as  well  as  the  collector 
current  noise  spectra  show  a  frequency  dependence  described  by  f"', 
with  Y  between  0.6  and  1.4.  This  indicates  to  us  that  there  is  a 
considerable  amount  of  generation-recombination  noise  present.  It 
is  well  known  that  the  exposed  base  area  is  likely  to  have  many 
surface  states  that  can  act  as  trapping  centers  for  the  carriers. 
It  is  not  useful  therefore  to  try  to  extract  values  of  the  Hooge 
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parameter  from  these  measurements.  The  magnitude  of  the  base 
current  noise  is  proportional  to  the  current  squared,  deviating 
from  a  junction  noise  theory  developed  by  Kleinpenning  [8]  that  did 
apply  to  the  silicon  BJTs.  These  observations  are  reminiscent  of  a 
surface  type  1/f  noise  [9] .  The  Rockwell  device  is  noisier  than  the 
Tektronix  device  at  the  same  base  current. 

The  collector  current  noise  spectra  show  again  that  the 
Rockwell  device  is  much  noisier,  about  100  times,  than  the 
Tektronix  HBT.  Due  to  the  relative  small  current  gain  of  the  HBTs, 
when  compared  with  the  BJT,  it  is  difficult  to  supply  sufficient 
feedback  to  reduce  the  contribution  of  the  noise  sources  at  the 
input  of  the  transistor.  As  a  result  we  cannot  estimate  the 
magnitude  of  the  source  ice- 

The  HBTs  are  inherently  noisier  than  the  BJTs  we  measured  on. 
The  1/f  noise  sources  are  generally  masked  by  g-r  noise  processes. 


3.2  Noise  in  JFETs  (Master's  Thesis  loannis  Stephanakis) 

(for  Abstact  see  p.  C5) 

Measurements  of  the  channel  noise  in  JFETs  are  analyzed  using 
a  model  in  which  the  channel  is  split  up  in  an  ohmic  part,  where 
the  crrier  velocity  is  proportional  to  the  electric  field  strength, 
and  a  saturated  part,  where  the  carrier  speed  is  the  saturation 
velocity.  In  this  model  the  value  of  Oh  is  assumed  to  be  constant 
throughout  the  channel  and  the  noise  measurements  are  used  to 
extract  this  parameter.  Long-channel  JFETs,  both  p-  and  n-channel, 
are  studied  as  well  as  short  n-channel  devices. 

The  30-iim  long  n-channel  JFETs  exhibited  values  for  the  Hooge 
parameter  on  the  order  of  10'®  for  drain  currents  of  about  1  mA. 
These  values  are  consistent  with  normal  collision  noise  and  Umklapp 
1/f  noise. 

The  30-p,m  long  p-channel  JFETs  showed  values  of  (5-10)  xl0‘®  for 
the  Hooge  parameter  for  drain  currents  of  1-2  mA.  These  values  are 
much  larger  than  expected  for  incoherent  state  1/f  noise  but  much 
smaller  than  the  values  expected  for  coherent  state  1/f  noise. 
Possibly  the  extracted  values  point  to  an  operating  regime  where  a 
transition  is  observed  from  one  mechanism  to  the  other. 

The  2-nm  long  n-channel  JFETs  featured  Hooge  parameter  values 
in  the  range  from  (3-11)  xl0"^°,  about  an  order  of  magnitude  smaller 
than  the  smallest  theoretically  expected  values.  During  these 
measurements  the  devices  were  operated  at  large  drain  currents,  20 
mA.  There  have  been  other  reports  in  the  literature  where  very 
small  values  of  the  Hooge  parameter  have  been  observed. 

It  is  expected  that  these  findings  will  be  published  shortly. 
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By  replacing  the  change  in  velocity  Av  by  the  low-frequency  Fourier  transform  F^IO)  of  the  elec¬ 
tron  acceleration  a(/).  Handel's  equations  for  the  Hooge  parameter  are  put  in  equivalent  forms 
that  are  not  only  applicable  to  collision  1//  noise  in  semiconductors  but  also  to  acceleration  I// 
noise  in  long  semiconductor  resistors.  To  prove  these  expressions  semiclassically,  one  evaluates 
first  the  bremsstrahlung  energy  dE  of  a  single  radiation  pulse  in  a  frequency  interval  df,  and  then 
defines  dn  —dE/hfas  the  number  of  photons  of  a  single  radiation  pulse  in  a  frequency  interval  df 
and  finally  dr  =dn  /t„  =dE  /hfr,  as  the  rate  of  photon  emission  in  a  single  pulse  in  a  frequency  in¬ 
terval  df.  It  is  then  found  that  the  expression  for  dr  already  contains  the  Hooge  parameter  a^.  It 
thus  seems  that  the  Hooge  parameter  depends  only  on  the  bremsstrahlung  emission  process  but  not 
on  the  details  of  the  electron-photon  interaction.  This  may  explain  why  Handel’s  expressions  for 
Off  so  often  agree  with  experiment.  One  must  now  bear  in  mind  that  the  elementary  current  event  is 
described  by  a  current  pulse  Ht)  of  duration  having  a  Fourier  transform  F,(0).  If  one  next 
defines  S"{f)=dr/df,  then  the  current  noise  spectrum  is  obtained  by  multiplying  S"(f)  first  by 
2Ff(0),  to  obtain  the  effect  of  a  single  elementary  event  per  second,  and  then  multiply  by  X,  the  num¬ 
ber  of  elementary  events  per  second.  This  leads  immediately  to  the  Hooge  equation  and  to  the 
Hooge  parameter  Off.  Van  Rheenan 

Reprinted  with  permission  of  the  publisher. 

Pages  1806-1809,  from  Physical  Review  B 

Vol40#3  by  A.  D.  Van  Rheenan  ,  ©  1989  The  American  Physical  Society . 


I.  INTRODUCTION 

According  to  Handel's  quantum  1 //-noise  theory'-^ 
the  \/f  current  noise  in  electronic  devices  is  associated 
with  the  low-frequency  bremsstrahlung  emitted  by  ac¬ 
celerated  or  decelerated  electrons.  This  I//  noise  is 
characterized  by  the  Hooge  parameter  a„  which,  in  turn, 
is  defined  by  the  Hooge  equation^ 

S,(f)  =  I^a„/fN  ,  (1) 

where  S,( /)  is  the  spectrum,  /  the  current,  /  the  frequen¬ 
cy,  N  the  number  of  carriers  in  the  system,  and  the 
Hooge  parameter.  Handel's  quantum  I //-noise  theory 
gives  expressions  for 

van  der  Ziel  et  al.  have  experimentally  verified 
Handel’s  expressions  for  for  a  number  of  cases.*  To 
understand  what  is  involved  here,  we  start  from  Handel’s 
equation  for  collision-frct  devices  (vacuum  tubes).  In  that 
case. 


from  experimentally  verifying  or  refuting  Eq.  (2) 
Verification  means  that  Eq.  (2)  has  heuristic  validity  but 
does  not  guarantee  the  validity  of  its  derivation. 

van  der  Ziel’  *  recently  gave  a  semiclassical  derivation 
of  Handel’s  expression  for  Off  and  introduced  a  new 
I //-noise  source  in  semiconductors,  the  so-called  ac¬ 
celeration  \/f  noise.  This  is  the  semiconductor  counter¬ 
part  of  1  //  noise  in  collision-free  devices  and  can  easily 
be  incorporated  into  Handel’s  schematic. 

It  is  the  aim  of  this  note  to  extend  and  generalize 
Handel's  equations  to  new  cases  and  to  establish  the  de¬ 
gree  of  validity  of  the  various  steps  in  van  der  Ziel’s  semi¬ 
classical  theory. 

II.  EXTENSION  AND  GENERALIZATION  OF  HANDEL’S 
EXPRESSIONS  FOR  an 

We  now  give  various  generalized  expressions  for  the 
Hooge  parameters  of  electron  devices.  Then  we  apply 
Eq.  (2)  to  old  and  to  new  situations. 

A.  Cases  where  the  electrons  flow  in  bunches 


Here,  c  is  the  velocity  of  light,  Av  is  the  change  in  veloci¬ 
ty  along  the  carrier  path,  and  ao=  is  the  fine-structure 
constant  for  electrons.’ 

For  collision-dominated  devices  involving  a  single 
collision  process,  Au’  must  be  replaced  by  Au^,  where  the 
averaging  is  carried  out  in  wave-vector  space  over  all 
scattering  events.  For  multiple<ollision  processes  the 
discussion  is  easily  extended.*  Often  alternate  equivalent 
formulations  are  necessary. 

Handel’s  theory  of  quantum  1  //  noise  is  not  generally 
accepted.  However,  this  should  not  prevent  anybody 


In  cases  where  the  charges  q  carrying  the  current  differ 
from  the  electron  charge  e,  Eq.  (2)  has  been  extended  as 
follows.  The  charges  flow  through  the  device  as  a  unit 
and  hence  generate  bremsstrahlung  as  a  unit.  Const 
quently,  the  fine-structure  constant  Oq  for  electrons  must 
be  replaced  by  the  fine-structure  constant  a  for  charges  q, 
since  a  is  proportional  to  q^,  one  expects  a=ao(q/e)’ 
Hence  Eq  (2)  should  be  written 

\i 

(3! 


4ao 

1 
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where  Uj=fi^E  and  r  is  the  average  time  between  col¬ 
lisions, 

•  (8) 

Consequently, 

^0  F(0)^  ^0  d  ^ 

=  - - i— =  ~r"  — 

3tr  3ir  cr 

We  therefore  see  that  acceleration  \/f  noise  fits  into 
Handel’s  extended  schematic,  and  that  it  does  not  require 
the  semiclassical  theory. 

Experimentally,  t~  10“'^  sec,  in  reasonable  agreement 
with  Eq.  (8).  For  d  >200  /im.  Off  may  saturate  to 
Handel's  coherent-state  value,  (2-5)X  10“^.''  The  experi¬ 
ments  were  performed  on  large  numbers  of  MOSFET’s 
operating  in  the  linear  mode,  but  not  all  MOSFET  types 
show  the  effect. 

III.  BREMSSTRAHLUNG  NOISE  SPECTRA 

We  now  make  a  Fourier  analysis  of  the  bremsstrahlung 
pulse  emitted  by  an  accelerated  charge  conglomerate  of 
charge  q.  As  shown  in  textbooks  on  electromagnetic  ra¬ 
diation,  the  bremsstrahlung  pulse  P(t)  of  such  a  charge 
conglomerate  is'^  (Larmor  equation) 

F(f)=^|i{/)|^  (9) 

3c  ^ 

for  0<t  <T,,  where  is  the  duration  of  the  pulse  and 
alt)  is  the  acceleration. 

The  Fourier  transform  of  a  (t)  is 


/''[a(t)]'dt=2  r\F,{j<o)\^df 


C 

Fa0<u)=J  a(/)exp(  — 7a»t)dt 

•  n 


=  FJ0) 


F.(;«u) 
FJO)  ’ 


where 


F«<0)=  f'MiOdt 

»  n 


=  v(t,)-v(0) 


is  the  low-frequency  value  of  F^ija)  and  Fjj(o)/FJ0) 
is  the  complex  "form  factor”  of  the  pulse.  It  is  carried 
along  to  make  f  g\Fg{ja))\^ df  converge  at  the  upper 
limit;  for  most  frequencies  of  practical  interest,  the  form 
factor  is  unity  and  can  be  omitted. 

We  now  evaluate  the  total  energy  E  in  the  radiation 
pulse  by  applying  Parseval’s  theorem.  According  to  the 
definition  of  E 

E=  f  '  Pi t)dt 
■'0 

=  ^  r*|F.(yai)l'd/ ,  (11) 

3c 

since  according  to  Parseval’s  theorem 


The  energy  dE  in  a  frequency  interval  d/  in  a  single  pulse 
is  therefore 

dE  =  ^\FJjco)\^df=SEi/)df,  (11a) 

where  S^l  f)  =  dE /df  ={4q^/'ic^)\FJja})\^  is  the  energy 
spectrum  of  photons  of  frequency  f  in  a  single  pulse  of 
duration  t^. 

Dividing  ( 1  la)  by  h/  yields,  for  the  number,  dn,  of  pho¬ 
tons  in  a  frequency  interval  d/  for  a  single  pulse  of  dura¬ 
tion  T., 


• = j/- = c' 


hf  3c' 


d/=5;(/)d/  , 


where  S’„(f)=^dn  /hf  =  iAq^ /‘ic^)[\Fg(j(o)\^  /hf]  is  the 
number  spectrum  of  photons  of  frequency  f  in  a  single 
pulse  of  duration  r^. 

Dividing  (11b)  by  yields,  for  the  rate,  dr,  of  photon 
emission  in  a  frequency  interval  df  for  a  single  pulse  of 
duration  t^, 

,  dn  dE 
dr  =  — =-r7~~ 

_  4q^  \FJj(o)\^ 

3c'  hfr,  ^ 

=s;(/)d/ ,  (lie) 

where  S,'(/)  =  dr/d/=(4q'/3c')[|F,(;(u)|'/h/T^]  is  the 
spectrum  of  the  photon-emission  rate  for  photons  of  fre¬ 
quency  /  in  a  single  pulse  of  duration  r,. 

The  spectra  S£{f)  =  dE/df,  Sl,if)  —  dn/df,  and 
S'( /)  =  dr/df  have  no  further  meaning  than  that  given 
by  the  definitions.  Those  who  object  to  the  symbol  S  may 
replace  it  by  any  symbol  of  their  liking.  Usually 
|F,(  )|'  =  F,'(0)  for  all  frequencies  of  interest. 

So  far,  we  are  on  firm  ground.  We  now  try  to  connect 
Sl,i  f)  and  (or)  S'^i  f)  to  the  current  noise  spectrum  S,( /). 
To  that  end.  we  introduce  the  Hooge  parameter  Off, 


where  aQ  =  2ne^ /he,  so  that  4ao/3ir  =  (4/8ff)(2rre'/lic) 
=  8e'/31ic.  In  this  notation. 

Off 


B3 


s;i/i= 


5;(/)=-^  =  — ^2e' 


2/r,  ’ 

S/if)  _  Off 


40 


EXTENSIONS  OF  HANDEL’S  1//-NOISE  EQUATIONS  AND  .  . . 


180V 


The  occurrence  of  the  factor  iq/e)^  has  been  verified  ex¬ 
perimentally,  both  in  vacuum  photodiodes’  and  in 
secondary-emission  multipliers,'®  so  that  it  has  at  leas, 
heuristic  validity. 

B.  Replacing  Av 

We  next  turn  to  the  change  Av  in  velocity  along  the 
electron  path  in  a  semiconductor,  introduced  in  Eq.  (2), 
to  further  elucidate  its  meaning.  If  a(t)  is  *he  accelera¬ 
tion  of  the  electrons,  and  if  r,  is  the  carrier  transit  time, 
then 

Au  =  f  *a(r)d/  =  F(0)  (4) 

■'o 

is  the  low-frequency  Fourier  transform  of  a(f).  Hence 
Eq.  (2)  becomes 

[F(0)1'  , 

We  shall  see  that  this  expression  has  many  applications. 

C.  Collision  I  / /  noise 


We  first  apply  it  to  collision  1//  noise.  Here,  is  the 
time  between  collisions  and  i(t)  =  {h  /m*)dk/dt,  where 
k  is  the  wave  vector,  so  that  F(0)  =  (h /m  *  )Ak.  This 
must  now  be  applied  to  the  various  collision  processes.* 
Process  1:  normal  collision  process  due  to  electron- 
phonon  interactions  in  the  conduction  band.  Here,  F^(0) 
must  be  evaluated  by  averaging  contributions  Ak  over  all 
collisions. 

Process  2;  umklapp  l/f  noise  by  transitions  of  elec¬ 
trons  into  the  next  Brillouin  zone.  The  lattice  thereby  ex¬ 
changes  an  average  momentum  h/a  with  the  electron, 
where  a  is  the  lattice  spacing.  This  process  can  only 
occur  in  degenerate  semiconductors  or  in  narrow-gap 
semiconductors.  Hence, 


4ao 

h 

"17 

m  *ac 

Process  3:  umklapp  and  intervalley  scattering  l/f  noise 
in  n-type  Si.  In  n-type  Si  the  conduction  band  has  six 
equivalent  valleys.  This  leads  to  the  following  collision 
processes:  (a)  normal  collisions  involving  transitions 
within  the  same  valley;  (b)  there  are  also  transitions  possi¬ 
ble  to  an  opposite  valley  [they  are  also  called  umklapp 
transitions  and  give  rise  to  an  expression  similar  to  Eq. 
(6)];  and  (c)  in  addition,  there  are  transitions  possible  to 
an  adjacent  valley  (they  are  called  intervalley  transitions 
and  give  rise  to  an  expression  similar  to  Eq.  (6)]. 

It  should  be  noted  that  each  valley  has  one  opposite  val¬ 
ley  and  four  adjacent  valleys  and  that  each  adjacent  val¬ 
ley  contributes  to  the  intervalley  1//  noise  equally  and 
independently. 

Since  there  can  now  be  multiple-collision  processes,  the 
definition  of  must  be  extended.*  To  that  end  one  at¬ 
tributes  a  mobility  ^i,  and  a  Hooge  parameter  «//  to  each 
process,  and  writes 


I 


=I 


5fi, 


Mr 


iL 

M, 


8^ 


au  = 


- - —  y  - - - 

3fl- 


(7a) 


In  process  2  or  3(b)  one  finds  /i/^,  =exp(  -  Qg/AT), 
where  0^  is  the  Debye  temperature.  In  n-type  Si  both 
processes  3(b)  and  3(c)  contribute;  this  has  not  been  prop¬ 
erly  analyzed  so  far  and,  therefore,  the  accurate  compar¬ 
ison  between  theory  and  experiment  still  needs  to  be  car 
ried  out.  The  experiments  yield  =2  5X10  *  at 

r  =  300  K  and  the  temperature  dependence  of  a„  is 
compatible  with  the  factor  exp(  —Q^/IT).  In  p-type  Si 
there  is  only  one  central  valley,  so  that  only  normal  col 
lision  1 // noise  contributes.  Experimentally,  O;,  ^  10  * 

at  7'  =  300  K,  in  reasonable  agreement  with  theoretical 
estimates.’ 

We  thus  see  that  the  introduction  of  the  parameters 
F(0)  gives  a  generalized  description  of  all  collision  \/f 
noises. 


D.  Accelerarion  1  / /  noise 

Next  we  turn  to  acceleration  \/f  noise  in  semiconduc¬ 
tor  resistors.  The  duration  of  a  collision  is  of  the  order  of 
10~'’  sec  and  the  average  time  between  collisions  is  of 
the  order  of  10“'^  sec  in  most  semiconductors.  During 
the  time  between  collisions,  the  electrons  arc  accelerated 
by  the  applied  electric  field  E  and  emit  bremsstrahlung 
During  the  much  shorter  collision  times  the  electrons  are 
accelerated  or  decelerated  by  the  collision  field  and 
also  emit  bremsstrahlung.  The  electron  interacts  with  its 
emitted  bremsstrahlung,  and  produces  current  l/f  noise 
in  either  case.  Note  that  E  has  a  constant  magnitude  and 
direction,  whereas  E^  has  a  fluctuating  magnitude  and 
direction. 

The  bremsstrahlung  components  due  to  individual  col¬ 
lision  processes  are  independent  and  add  quadratically  so 
that  is  independent  of  the  device  length  d,  whereas 
the  bremsstrahlung  components  generated  between  col 
lisions  are  fully  correlated  and  add  linearly,  so  that 
vanes  as  d^  (see  below),  where  d  is  the  length  of  the  de¬ 
vice  in  the  direction  of  the  field.  For  long  devices  the 
latter  process  therefore  predominates;  the  noise  is  called 
acceleration  1  //  noise  and  was  observed  for  long  Si 
metal-oxide-semiconductor  field-effect  transistors 
(MOSFET's),  both  p  and  n  type.  The  collision  noise  itself 
predominates  for  short  devices,  since  it  is  independent  of 
d. 

We  now  calculate  F(0)  for  the  acceleration  I  //  noise 
If  E  is  the  applied  field,  the  acceleration  aiti  —  eE/m”, 
where  m  *  is  the  effective  mass  Since  Uj  =  Ax/At  —n„E 
is  the  drift  velocity  of  the  carriers,"  we  can  write 


Fioi- r 


eE  At  , 

- - dx 

0  m  *  Ax 


=  f 

'  d 

eE 

J 

0 

-  f 

‘d 

eE 

'  J 

0  m*p^E 

dx  —d/r  , 
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FIG.  1.  The  current  pulse  Ht)  =  ev(t)/d  modulated  by  the 
I//  noise  caused  by  the  random  emission  of  photons  by  the 
electron. 


and  5;'(/)  is  the  noise  for  the  case  of  one  electron  per 
second  where  A.  =  / /e.  Consequently, 

s;{f)  =  s;{/)2e^  .  <13) 

s}{/)  =  s;{f)—  .  (13a) 

We  shall  see  that  (13)  has  physical  meaning,  but  that 
(1 3a)  has  not. 

Due  to  the  random  emission  of  photons  by  the  elec¬ 
tron,  the  pulse  /(f)=«»(r)/d  has  1  //  noise  superimposed 
on  it  (Fig.  1).  The  ensemble  average  at  the  instant  f, 
lii)=ev{l)/d  gives  shot  noise  5,( f)  =  2e^k,  whereas  the 
1//  noise,  described  by  modulates  it.  Therefore, 

Sjif)  =2c^S'( /),  and,  hence, 

5,(/)  =  2e^A5;(/)=a„-^  , 

/  'a 

as  had  to  be  proved. 

On  the  other  hand,  2e^/r„  represents  the  shot  noise  for 
the  case  in  which  there  is,  on  the  average,  one  electron  in 
V  the  system;  however,  the  combination  (13a)  makes  no 
physical  sense. 


IV.  CONCLUSIONS 
We  conclude  the  following. 

(1)  Handel's  expressions  for  the  Hooge-paramt'er  (Off ) 
1  / /  noise  can  be  generalized  to 


4ao 

1 

2 

FJO) 

3Tr 

e 

C 

where  the  current  is  carried  by  charge  conglomerates  and 
F^IO)  is  the  low-frequency  Fourier  transform  of  the  car¬ 
rier  acceleration  a(t).  This  holds  for  collision-free  and 
collision-dominated  devices  including  acceleration  1// 
noise  in  semiconductor  resistors. 

(2)  This  formula  can  also  be  derived  semiclassically  in  a 
more  rigorous  manner  than  in  van  der  Ziel’s  original  pa¬ 
pers.  To  that  end,  one  first  derives  the  energy  dE  of  a 
single  radiation  pulse  in  a  frequency  interval  df  with  the 
help  of  Parseval’s  theorem.  Next,  one  defines 
dn  —dE/hf  as  the  number  of  photons  of  a  single  pulse 
in  a  frequency  interval  df  and,  finally,  defines 
dr  =dn  /T^=dE  /hf  as  the  rate  of  photoemission  in  a 
single  pulse  of  duration  in  a  frequency  interval  df.  It 
is  then  found  that  the  expression  S',’{f)  =  dr /df  already 
contains  the  Hooge  parameter  a^. 

(3)  If  S,{  f)df  is  the  current  noise  in  a  frequency  inter¬ 
val  df  and  A.  the  number  of  pulses  per  second,  then 
S'j'{f)df  =[S,(/)/A.]d/  is  the  current  noise  fo*-  a  single 
pulse  in  a  frequency  interval  df.  Since  each  emitted  pho¬ 
ton  contributes  both  to  S''(  f)df  and  to  S','(  f)df,  the  two 
icrms  are  proportional.  Hence,  if  occurs  in  S'fif)  it 
will  also  occur  in  5/''( /)  and,  hence,  in  S,(/). 

(4)  We  can  evaluate  S/(/)  directly.  An  elementary 
event  gives  rise  to  a  radiation  pulse  Pit)  that  leads  to 
S'/if)  and  to  a  current  pulse  lit)  both  of  duration  r^. 
Let  i{t)  have  a  low-frequency  Fourier  transform  F,(0); 
then,  S;'(/)  =  S;'(/)F,^(0)  and,  hence,  S,{f)  =  S;'if)k 
—  S'/if )F^{0)k.  This  leads  immediately  to  the  Hooge 
equation. 
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I.  IK  I RODUCTION 


In  electronic  devices  the  relative  1  //  noise  spectra  S,{f)/l^ 
can  usually  be  represented  by  a  Hooge-type  equation': 


Si(f)  ^  Off 

~P  7^’ 


(1) 


where  the  exponent  y  is  close  to  unity,  IV  is  the  number  of 
carriers  in  the  system,  and  is  a  constant  parameter,  the 
so-called  Hooge  parameter;  it  has  the  dimension  (Hz)’'"  '  It 
can  always  be  defined,  and  can  be  used  as  a  measure  for  the 
noisiness  of  the  system  under  investigation,  but  may  not  al¬ 
ways  have  physical  meaning.  This  is,  e.g.,  the  case  for  classi¬ 
cal  1// noise  caused  by  surface  traps.  In  other  cases  can 
often  be  identified  with  distinct  physical  processes,  each  re¬ 
sulting  in  a  characteristic  value  for  a„,  dependent  on  the 
process  involved. 

For  example,  Hooge'  found  a„  =2x  10“’  for  relatively 
long  semiconductor  resistors,  independent  of  the  material 
and  the  doping,  as  long  as  the  doping  was  not  too  large.  This 
may  be  identified  with  Handel’s  coherent  state  1// noise  pro¬ 
cess  with 


=  2a(/rr  =  4.6X  10  ’,  (la) 

where  On  =  1/(137)  is  the  fine-structure  constant.^ 

In  short  devices  a„  can  be  much  smaller.'*’  It  can  be 
identified  with  collision  processes  (normal  collision  I// 
noise,  umklapp  I//  noise,  and  intervalley  scattering  I// 
noise)  or  hole-electron  pair  generation  and  recombination 
1//  noise  processes.  For  a  single-collision  process  Handel 
has  proposed  the  quantum  1// noise  equation'’: 


^  ‘*^0  _ 

}rr 


(2) 


where  Xv  is  the  vectorial  change  in  velocity  during  the  colli¬ 
sion  process.  If  several  collision  1  If  noise  processes  are  oper¬ 
ating  simultaneously,  Kousik  et  aV  proposed 


Oh  =  <2a) 

where  a,„  is  associated  with  the  ;th  collision  process,  and  p., 
is  the  mobility  associated  with  that  process  In  generation- 
recombination  1// noise  processes  Air  can  be  calculated  as 
lE/m*.  where  E  is  the,  energy  associated  with  the  pro¬ 
cesses*’  (see  the  Appendix). 


If  normal  and  umklapp  processes  are  both  present,  Eq. 
(2a)  yields' 


5tr  \m*ac)  V  IT/ 


(21  > 


Here  m*  is  the  effective  mass,  a  is  the  lattice  spacing,  and  Ojy 
is  the  Debye  temperature  of  the  material,  Eq.  ( 2b)  takes  into 
account  that  in  an  umklapp  collision  the  lattice  takes  up  or 
gives  up  a  momentum  h /a  and  that  the  weight  factor 
(it/H, )'  =  exp(  —  9[f/2T).  The  small  contribution  of  nor¬ 
mal  collision  1// noise  has  been  neglected. 

We  now  discuss  how  can  be  measured  experimental¬ 
ly.’  If  the  noise  can  be  represented  by  a  “lumped”  model,  as 
in  vacuum  tubes,  Eq.  ( 1 )  can  be  used  directly,  and 


f^S,(f)S  S,(\)N 

r  n 


(3) 


as  long  as5,(/)  =  5,(1)//’'  with  constant  y.  There  is  one 
difficulty  left,  however.  Equation  ( 1 )  implies  iV>  1,  but  in  a 
vacuum  tube  or  solid-state  diode  N can  have  a  value  less  than 
unity.  In  that  case  we  put  /  =  Ne/r,  where  r  is  a  transit  time 
or  a  lifetime;  consequently  Eq.  ( 3 )  becomes 


Oh 


•y,(  l)r 
le 


or  5,  ( 1 )  =  a 


(3a) 


This  equation  remains  valid  for  N<  1  and  hence  Fxj.  (3a) 
should  replace  Eq.  ( 3 )  in  that  case. 

In  many  diodes,  however,  Eq.  ( 1 )  must  be  written  in  the 
distnbuted  form’ 


S,(x,f)  g„ 

l\x)  f’'Mx)Xx  ' 

for  a  section  Xx  at  x,  where  N{x)  is  the  number  of  carriers 
per  unit  length  at  x  and  I(x)  is  the  minority-carrier  current 

at  X. 

One  now  introduces  a  random  noise  course  //(.t,t ).  evalu 
ales  Its  spectrum  S„(xj(',f)  from  Eq.  (4),  and  so  deter¬ 
mines  the  spectrum  S,  ( / )  in  the  external  circuit  by  integra 
lion.  This  yields 

5,  ( / )  =  a  ff  E{y),  ( 5 ) 

where  Fisthe  (back)  bias  of  the  diode  and /"(  F)  is  evaluated 
in  the  integration  process.  The  time  constant  r  can  be  mea- 
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sured  by  measuring  the  device  admittance  as  a  function  of 
frequency.'” 

If  all  parameters  are  known,  a„  can  be  evaluated. 


II.  n  Hg,_,Cd,Te  PHOTODIODES 

We  now  turn  to  various  Hg,  _,Cd,Te  n*-p  diodes  sup¬ 
plied  by  Rockwell  International.'”  They  are  relatively  long 
diodes,  with  a  cadmium  content  of  30%;  they  have  con¬ 
tacts  of  1 50-/xm  diameter.  The  length  of  the  p  region  { u;, )  is 
larger  than  the  diffusion  length  {L„).  The  current  flow  is  by 
diffusion,  as  follows  from  the  diode  admittance 

=  (6) 

The  low-frequency  admittance  go  is  equal  to  dl /dV. 

Measuring>'(  jeo)  near  zero  bias  yielded  r~2  X  10“^  s  at 
193  K  in  good  agreement  with  the  Honeywell  lifetime  ta¬ 
bles.'*  Far  from  zero  bias  the  experimental  data  could  be 
fitted  to  the  function  g„(  1  -f-ytur) with  7-0.9.  In 
that  case  the  values  of  r  deviated  from  the  Honeywell  tables 
and  were  considered  unreliable. 

Since  the  noise  seems  to  be  diffusion  \/f  noise,  the  follow¬ 
ing  Hooge  parameters  are  expected"; 

( i )  Coherent  state  l//noise:a'„  =  4.6  X  10~’  [Eq.  (la)]. 

(ii)  Umkiapp  1//  noise  to  next  Brillouin  zone: 
Off  =  5X  10~’  at  T  =  273  K  (Eq.  (2b) ]  and  jc  =  0.30. 

(hi)  Normal  collision  1// noise:  =  10“^  at  r=:  77  K 

and  X  =  0.20  (estimated). 

For  normal  elastic  collisions^ 


_  402  ^^7* 

~  3tr  mV  ’ 


(6a) 


corresponding  to  o^=  0.34x10"’  at  T  =  77  K;  for 
X  —  0.20,  m  =  0.0073  m^,  where  m^  =  9. 1 1  X  10"’'  kg,  the 
mass  of  the  free  electron.  If  the  inelasticity  of  the  collisions 
was  taken  into  account  in  the  case  of  electrons  in  silicon  o^ 
was  found  to  be  —3.5  times  larger.  Assuming  the  same  cor¬ 
rection  factor  for  Hg,  .,Cd,Te  yields  a„  —  1.2x  10"’  in 
close  agreement  with  our  estimate  of  10"’. 

When  we  applied  Fq.  (5)  and  plotted  5,  (/)/[/F(  V)  j  at 
a  fixed  frequency / vs  V,  we  found  that  S,  ( / )/[IF{  V)  ]  was 
constant.  The  theory  predicted  that 


3 


J_ 

2a 


- r  ln(l  +a); 

a-* 


a  =  expieV/kT)  —  1, 


(7) 


so  that  f(0)  =0,  as  follows  from  a  Taylor  expansion  of 
Inf  1  -b  a)  for  small  a  up  to  the  a’  term.  This  function  was 
used  in  Fig.  1.  Since  S,(\)/[IF{V)\  was  independent  of 
bias,  a^/r  was  apparently  independent  of  bias  (Fig.  l),even 
when  Yo  had  values  of  0.7  to  0.9. 

The  spectra  were  of  the  form  I//’'  with  y  very  close  to 
unity  (Fig.  2),  as  expected  for  quantum  I// noise.  Figure  3 
shows  a„aai  function  of  back  bias,  indicating  that  a  high  Yo 
and  a  high  a„  go  together.  Since  a^/r  is  constant,  the  error 
in  Oh  is  due  to  a  corresponding  error  in  r. 

Measuring  5,( /  )  at  1 13  and  at  193  K,  we  found  in 
all  units  near  zero  bias  a^,  =  (3-5)  X  10"’,  indicating  that 
the  1//  noise  was  cohennt  state  l/f  noise.  Similar  values  of 

£7 


Fig.  I.  The  relation  between  [S,( /)  -  S, (  »  )  ],'1/F(  10 1  and  the  bias 
voltage  F  for  n  *  -p  diode  3-454  is  represented  by  a  horizontal  bar. 


otff  were  obtained  for  y  =  0.7-0.9  by  using  the  r  values  of  the 
Honeywell  lifetime  tables. 

At  273  K  one  unit  had  a„  =  5x  10"’  as  expected  for 
umkiapp  l/f  noise  to  the  next  Brillouin  zone.  This  unit  had 
r  =  1.2  X  10"’  s  according  to  the  lifetime  tables. 

What  happens  if  In  that  case  the  admittance  is'’ 


tanh(2y<yrj ) 


(8) 


where  =  wl/lD., ,  is  the  length  of  the  p  region,  and  D„ 
is  the  diffusion  constant  of  the  electrons  in  thep  region.  Cal¬ 
culating  5,  (/ )  and  applying  Eq.  (5)  yields 


S,(f)  =a„ 


elF(V) 


F(  F)  =  J  In 


_ _ 

exp(  -  eV/kT) 


(9) 


where  is  the  recombination  velocity  of  electrons  at  the  p 
contact.  For  most  materials  is  of  the  order  of  10’  cm/s. 
Note  that  F(V)  ~\eV /kT  for  large  negative  V.  Also, 


Fig  2  The  spectrum  S,>  f)  for  diode  5. 1 28M  at  F  =  0  002  V  is  of  the  form 
I//’  with  yvery  close  to  unity 
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Fig.  3  Hooge  parameters  a„  vs  applied  voltage  K  Note  the  effect  of  the 
parameter  f,,. 


>  Z),  /u;^  in  most  cases.  ( For  D„  =  250  cm  Vs  and  =  5 
fxm,  D„/w^  =  0.5 X  10“  cm/s,  which  is  an  order  of  magni¬ 
tude  smaller  than  . ) 

How  can  the  operation  of  the  detector  be  improved?  To 
that  end  one  must  make  or  [  S,  (/)]''  ^  as  small  as  possi¬ 
ble.  This  means  that  one  must  make  (a„/r)  and  IF(  F)  as 
small  as  possible.'^ 

The  highest  value  of  (5x  10“  V  occurs  for  coherent- 
state  1//  noise;  this  fundamental  noise  source  should  there¬ 
fore  be  avoided. 

The  next  lower  value  of  (5  X  10“’)  occurs  for  umk- 
lapp  1//  noise.  It  might  be  achieved  by  making  the  device 
dimensions  somewhat  smaller,  and  would  give  a  value  of  a„ 
t  hat  is  ~  1 00  times  smaller.  Hence,  if  t  does  not  change  very 
much,  as  our  data  seem  to  indicate,  could  be  lowered  by  a 
factor  of  10.  This  needs  further  investigation,  especially  as 
far  as  the  value  of  r  is  concerned.  At  273  K  the  p  region  was 
nearly  intrinsic,  and  this  might  make  the  theory  of  Sec.  Ill 
applicable  when  =  w\/2D„ 

Can  the  normal  collision  1// noise  be  reached?  We  esti¬ 
mated  a„  2  10“’,  but  it  can  probably  only  be  reached  for 
very  small  e.g.,  w^<,5  pm.  But  in  that  case  is  very 
small  and  hence  CLff/r^  may  not  be  much  improved.  For 
example,  our  unit  at  T  =  273  K  has  a„  =  5x  10“’  and 
r=  1.2  X  10“ ’s,  so  that  a„/r  =  400  s  “'.  Suppose  now  that 
we  could  achieve  0^  =  10“’  at  =  5  pm.  Then 
Tj  =  w^i,/2D„  =  5x  10"  '“s  (ifwetakeZ),  =  250cmVs),so 
thata„/7-j  =  200  s  “'.  This  is  an  insignificant  improvement. 
So  we  should  search  for  the  normal  collision  1//  noise  in 
order  to  complete  our  catalog  of  noise  sources,  and  to  see 
whether  our  pessimistic  views  about  P^  arc  correct,  but  we 
should  not  keep  our  hopes  up  too  high.  Fortunately,  better 
ways  have  been  found’^of  achieving  lower  values  of  P^  (see 
Sec  III). 

The  diffusion  I// noises  in  n  *~p  Hg,  .  ,Cd,Te  diodes  are 
so  high  because  the  effective  mass  of  the  electrons  is  so  small 
(  =0.01  my);  as  a  consequence  (a„)j,*  >5x  10“’  Trap¬ 
ping  1  // noises,  however,  when  described  by  a  Hooge  param¬ 
eter,  usually  havea^  =  10““-10"’  so  that  those  noise  pro- 

Bg 


cesses  will  be  completely  masked  by  diffusion  1 //noise.  For 
that  reason  we  can  understand  why  only  one  Rockwell  unit 
showed  burst  noise  and  none  showed  trapping  1// noise. 


III.  Hg,_,Cd,Tep-Z-rt  DIODES 

In  order  to  study  the  p-i-n  structures  we  first  investigated 
long  p-i-n  silicon  power  diodes  as  far  as  admittance  and 
noise  were  concerned.'*  If  the  i  region  extends  from 
—  d<x<d  then  a  p-i-n  diode  is  called  “long”  if  d>L, 
where  L  is  the  diffusion  length  and  “short”  if  d  <L.  The 
admittance  had  a  complicated  structure  and  showed  two  re¬ 
sonances.  The  impedance  of  the  device  could  be  written  as 
1/K^  Z„,  where  =gj  +jci>Cj  is  the  admittance  of 

the p-i  junction  and  is  the  modulation  impedance  due  to 
the  modulation  of  the  carrier  density  in  the  /  region  by  the  ac 
current.  Here 


1 


H - +j(i}C, 

jo)Lo 


(10) 


where  gives  the  first  resonance  and  LqC  the  second;  g,, 
Cj,  and  Z-o  b*  functions  of  frequency. 

Next  we  solved  the  Fletcher  equations  for  ac  and  found 
this  yielded  a  term  for  the  high-frequency  admittance  of  the 
junction”: 


Y  —  Z.  tanh(<//Z.') 

’  I'  tanh(d/L)  ’ 


(11) 


where  go  =  dl  /dV,  L  ‘/L  =  ( 1  -(-  ytur)  “  and  r  is  the  car¬ 
rier  lifetime. 

For  long  diodes  (d  /L>  1 ),  tanh(<//Z.)  s  1,  and  tanh(d/ 
Z, ')  1,  so  that  at  lower  frequencies 


Yj=go-^  =  «o(  •  +>-r) =  g,  +;wC^,  (11a) 

corresponding  to  Eq.  (6);  here  g^  and  C,  are  frequency  de¬ 
pendent;  this  corresponds  to  Eq.  ( 6 ) .  For  shorter  diodes  ( d  / 
L<1), 

Y,  =go(Z,/Z.')^  =go(l  +yiyr)  =  go -I- /tuCp.  (11b) 

where  go  and  Cp  are  independent  of  frequency. 

In  addition  one  can  evaluate  Z„  as  a  function  of  frequency 
for  higher  frequencies.  This  leads  to  a  complicated  reso¬ 
nance  ford /Z  >  1,  as  described  by  Eqs.  (6)  and  ( 10)  and  to  a 
simple  resonance  described  by  Eq.  (lib)  and  Z„  for  d/ 
Z  <  1,  where 


1  juLn 

2  +ja)L„ 


(12) 


The  simple  resonance  is  now  determined  by  the  frequency 
<u„  =  (LoCp)~''^  (Fig.  4). 

Calculating  all  circuit  elements  numerically  we  can  evalu¬ 
ate  r  either  from  Y^,  when  Y^  predominates,  or  from  Z„ 
when  Z„  predominates.  Applying  this  to  the  long  Si  p  ^  -i-n 
power  diode  yielded  1.4x  10“’  s  at  T  =  400  K,  as  ex¬ 
pected  for  such  diodes,  whereas  for  the  Hg,  _  ,Cd,Te  diodes 
with  X  =  0.22  and  T  =  80  K  all  Santa  Barbara  units  gave 
r  =  10"’  s,  comparable  to  the  Honeywell  lifetime  tables.'" 

We  conclude  from  the  latter  result  that  the  Santa  Barbara 
diodes  are  relatively  short  (d  /Z  <  1 )  p-i-n  diodes  and  not 
n*-p  junction  diodes. 
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Fraquancy  (MHz) 


Table  I.  a„  values  for  different  diodes  at  T  =  80  K 


No. 

yAmV) 

/j(fiA) 

5;(/)  (aVHz) 

Off  ( r  =  10  ^  s) 

M 

-0.8 

-  3.5 

2.25X10"” 

2.0  X  10  ' 

G 

-  3.0 

-  3.5 

0.95  X  10” 

15  vIO  ' 

F 

-  20 

-  6.0 

4.3  xIO"” 

2.25X  10  ' 

3 

-  25 

-  7.0 

0.9  X  10  ” 

8.0  ,<  10  * 

K 

-  30 

-  10.0 

1.0  XlO"” 

3.0  X  10  ‘ 

H 

-40 

-  12.0 

0.8  XlO"” 

2  0  X  10  ‘ 

Fang”  found  much  lower  values  for  a„  than  either  the  co- 
hsrent-state  theory  (4.6x10“’)  or  the  umklapp  theory 
(1.7x  10“^)  would  predict  for  the  case."  Measurements  by 
Fang  and  He"  indicated  that  |F(  V)  |  was  independent  of  the 
back  bias  ( y\  (Fig.  5).  We  therefore  assumed /"(  V)  =  1  and 
obtained  the  results  shown  in  Table  I. 


Fig.  4.  The  conductance  G(  /)  vs  frequency  for  Santa  Barbara  Research 
Center  devices  at  80,  1 1 3.  and  193  K.  Crosses:  experiment;  full-drawn  line: 
theory  (r  =  10"’  s).This  representative  figure  seems  to  indicate  that  the 
devices  are  n~i-p  diodes. 


We  now  turn  to  the  1// noise.  The  p-i-rt  diode  Wf  noise 
has  not  been  evaluated  theoretically  but  by  analogy  with  the 
previous  discussion  one  would  expect  a  formula  of  the  type 


where  F{  10  is  usually  a  function  of  Kbut  F{  F)  may  not  be 
equal  to  zero  for  K  =  0.  Measuring  S,(/)  and  r  vs  V,  we  can 
thus  evaluate /  (  V)  from  the  data. 

For  long  Si  p-i-n  power  diodes  (<//£,>!)  Fang'*  found  at 
T  =  400Kthat5,  (/)/|/|  was  independent  of  the  back  bias 
voltage  V;  this  means  that  /(  F)  was  independent  of  F. 

Assuming  /(F)  =  I  yielded  a„  =  4. 1  X  10~’  in  good 
agreement  with  what  the  coherent-state  1 //noise  theory  pre¬ 
dicts  (4.6X  10“’).  A  more  accurate  theory  of  S,(f)  is,  of 
course,  needed  for  this  case,  but  the  simple  results  obtained 
so  far  look  plausible. 

For  shorter  Hg,  _  ^.Cd.Te  diodes  with  x  =  0.22  at  80  K, 


Voltaga(V) 


Fk».  i.  S,(/)/e\I\  — vtbackbiaf  for  diode  G 1 4,  indicating  that 


ft  F)  ia  independent  ofbiaa  (biaa  in  mV,  not  V). 


J.  Vac.  Sci.  TBCitnol.  A,  VoL  7,  No.  2,  Mar/ Apr  1M9 


The  spectra  were  of  the  form  1//’',  where)' deviates  signif¬ 
icantly  from  unity  (Fig.  6).  The  latter  is  incompatible  with 
quantum  1//  noise,  but  can  be  explained  by  superimposed 
trapping  processes  or  classical  modulation  effects. 

How  have  such  low  values  of  cr„  been  achieved?  Only 
because  diffusion  I// noise  processes  are  inactive  in  this  case 
(see  the  Appendix).  What  is  then  left  are  various  trapping 
processes,  probably  at  the  surface  of  the  p  region  (or  the 
surface  of  the  junction  space-charge  region),  or  at  other 
crystal  defects.  The  noise  is  not  fully  describable  in  terms  of 
quantum  1// noise  processes,  though  they  may  play  a  part; 
classical  modulation  effects  might  enhance  it,  just  as  is  the 
case  in  Si  bipolar  junction  transistors  (see  the  Appendix). 

The  a„  values  of  Table  I  vary  by  a  factor  of  10.  This  is 
explainable  in  terms  of  trapping  effects,  for  S,(f)  should  be 
proportional  to  the  trap  density  and  can  therefore  vary  from 
unit  to  unit  even  in  a  single  array. 

But  there  could  be  a  lower  limit  set  by  direct  hole-electron 
pair  generation  or  recombination.  We  can  make  the  follow¬ 
ing  estimate.  If  a  hole-electron  pair  is  generated,  the  average 
energy  of  the  free  electron  is  3A:  r/2  and  the  average  energy  of 


Fraquancy  (Hz) 

Fig.  6.  [5,(/)  -  5,(00  )|  vs  rrequency  for  diodes  G 14  indicating  that  the 
spectrum  is  proportional  to  \/p  with  I .  A  is  the  parameter  characteriz¬ 
ing  the  computer  fit. 
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the  free  hole  IkT /2,  so  that  the  total  energy  involved  is 
E—lkT and  hence  Ap  =  3E/m*,  or*® 


4«o/ 

2E  > 

1  _  4ao 

6k  T 

377-  ' 

V  m*c^  > 

'  3jr 

m*c^ 

(14) 


corresponding  to  the  elastic  theory  of  normal  collision  l/f 
noise.  Van  der  Ziel®  has  explored  the  possibility  of  adding  a 
term  to  E,  where  £,  is  the  gap  width.  It  is  unlikely  that 
this  conjecture  can  be  justified  theoretically. 

It  should  be  understood  from  Table  I  that  the  lower  limit 
for  a„  has  not  yet  been  reached.  Hence  Eq,  ( 14)  has  not 
been  verified  and  can  only  be  used  as  a  guide  for  interpreting 
data. 

Jones  and  Radford^”  found  a„  =  5x  10“’  for  similar 
units.  They  also  found  that  the  temperature  dependence  of 
the  1 //noise  was  proportional  to  exp(£'j/2AD,  as  expected 
for  g-r-type  l/f  noise  processes.  On  the  other  hand,  the  l/f 
noise  in  n*-p  Hg,_  ^Cd.Te  photodiodes  has  an 
exp(  E^/kT)  temperature  dependence,  as  expected  for  diffu¬ 
sion-type  \/f  noise  processes. 
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APPENDIX 
1.  diode 

We  can  elucidate  the  admittance  and  the  noise  properties 
of  (he  p-i-n  diode  by  considering  the  dc  solution  of  the 
Fletcher  equations.  Here />(x)  =«(  <:)  in  the  /  region.  Let  the 
i  region  extend  from  —  d  to  +  d  and  let  L  be  the  diffusion 
length  of  the  carriers.  Then,  according  to  Fletcher,” 


p(x)  =p(!) 


coshfjc/L) 
cosh(d  /L)' 


(  —d<x<d). 


(Al) 


with /?(ac)  =/»(/)  at  a:  =  ±d. 

For  long  diodes  (d^L)  p(x)  depends  strongly  on  jt;  there 
is  a  earner  density  gradient  at  the  end  points  ±  d  and  hence 
there  is  diffusion  at  the  boundaries.  It  is  therefore  not  sur- 
pnsing  that  Eq.  ( 1  la)  agrees  with  Eq.  (6).  For  short  diodes 
(d/L4,\)  p(x)  is  independent  of  x  (/>(x)  =/>(/)  for 
-  d  <x  <d\\  there  is  no  carrier  density  gradient  and  hence 
no  diffusion.  For  Eq.  (lib)  one  would  expect  go  =  /dV. 
Furthermore,  for  lack  of  diffusion  one  would  not  expect  dif¬ 
fusion  l/f  noise. 

According  to  Handel,  the  quantum  l/f  noise  due  to  bulk 
or  surface  traps  can  usually  be  written  in  the  form* 


_  4ao _ 

3ir  m*c^ 


(A2) 


where  E always  contains  a  term  (3kT /2)  or  ikT.  In  addi¬ 
tion,  there  can  be  another  energy  term  E'  due  to  potential 
barriers,  built-in  voltages,  etc.  But  it  is  doubtful  whether  val 
ues  of  a„  >  10“*  can  be  achieved  in  this  manner. 

Since  our  measurements  of  a„  are  in  the  range  2  x  lO  *- 
2  X  10“’,  it  seems  that  they  cannot  be  interpreted  by  quan¬ 
tum  effects.  We  believe  that  classical  modulation  of  (sur¬ 
face)  recombination  centers  by  the  fluctuating  occupancy  of 
traps  in  the  surface  oxide  might  explain  the  discrepancy. '®  In 
addition  the  frequency  dependence  of  most  spectra  also  rules 
out  a  quantum  l/f  noise  interpretation. 
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Fang  and  van  der  Ziel's  analysis  of  secondary  emission  1//  noise  in  secondary  emission 
pentodes  is  extended  to  16  individual  data  points.  It  is  found  that  by  proper  choice  of  the 
secondary  electron  path  length  ,  good  agreement  between  the  experimental  values  of  the 
Hooge  parameter  a„  and  the  theoretical  values  predicted  by  Handel’s  quantum  \/f  noise 
theory  is  obtained  for  14  data  points.  Handel’s  equations  for  the  Hooge  parameter  of  the  tubes 
thus  seem  to  have  heuristic  validity.  The  values  of  the  secondary  electron  path  length  differ 

somewhat  from  tube  to  tube  and  may  also  depend  slightly  on  bias;  these  effects  are  attributed 
to  the  electron-optical  system  formed  by  screen  grid,  dynode,  and  anode. 


In  I960  Schwantes  and  van  der  Zeil'  observed  second¬ 
ary  emission  1// noise  in  Philips  EFP60  secondary  emission 
pentodes  but  they  had  no  satisfactory  explanation  of  the  ef¬ 
fect.  In  1987  van  der  Ziel^  applied  tlie  predictions  made  by 
Handel’s  quantum  1// noise  theory’  to  the  problem  and  ob¬ 
tained  good  agreement  with  experiment  aOer  a  slight  modifi¬ 
cation  of  Handel’s  expression  for  the  Hooge  parameter  of  the 
system. 


*'  Now  M  Ihc  EE  Depaftment,  Univenity  of  South  Ftorida,  Tampa,  FI. 


The  I// noise  is  generated  in  the  region  between  the 
secondary  emission  dynode  ( d)  and  the  anode  ( a )  collecting 
the  secondary  electrons.  According  to  Hooge’s  equation^ 

=  a„{ll//N)  =  a„{eI,/fT^ )  =  AkTR„gl,,  ( I ) 

where  is  the  Hooge  parameter,  and  S  ■=  I,T^/e  the 
number  of  electrons  in  the  system.  Furthermore,  is  the 
electron  transit  time  from  dynode  to  anode,  /.  is  the  anode 
current,  =  3/,/dV^  the  transconductance  of  the  tube, 
and  R„  the  secondary  emission  noise  resistance  measured 
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by  Schwantes  and  van  der  Ziel. '  (The  paper  of  Fang  and  van 
der  Ziel^  has  the  obvious  misprint  g„„  =  dl^/dV^.)  Han¬ 
del’s  modified  expression  for  the  Hooge  parameter  a„  is^ 

a„  =  (4a„/3ir)6^(A«Vc^),  (la) 

where  a„  =  1/(137)  is  the  fine  structure  constant  for  elec¬ 
trons  in  the  system,  S  is  the  secondary  emission  factor,  Au  is 
the  velocity  change  of  the  secondary  electrons  along  the  elec¬ 
tron  path  from  dynode  to  anode,  and  c  is  the  velocity  of  light. 
The  modification  consists  in  adding  a  factor  (5^  to  (la);  it 
comes  about  because  the  secondary  electrons  are  emitted  in 
bunches  of  charge  Se,  so  that  the  fine  structure  constant  of 
the  bunches  should  be The  factor  <$Mn  (la)  was  first 
introduced  empirically  and  subsequently  justified  by  the  ar¬ 
gument  given  above. 

If  the  secondary  electrons  are  emitted  with  near-zero 
velocity,  {V,  —  V^)  is  the  potential  difference  between 
dynode  and  anode  and  is  the  length  of  the  electron  path 
in  the  dynode-anode  system,  we  obtain 

Ai;^=  (2e/m)(F,  -  FJ, 

’■d.  =  2</,./Au  =  ld^/[le(  V,  -  V^)/mY'\  ( lb) 

van  der  ZieP  found  good  agreement  between  theory  and 
experiment  by  choosing  =  0.50  cm,  which  was  compati¬ 
ble  with  the  device  dimensions. 

By  substituting  ( lb)  and  ( la)  into  ( 1 ),  obtain 

=a„  Au 

2/dd. 


4go  (2e/m)^^^ 
llT 


2/ad. 


Fang  and  van  der  Ziel  measured  the  spectral  intensity  of 
the  fluctuations  in  the  anode  current  in  two  EFP60  tubes. 
They  determined  5,  (/)/( K,  —  at  constant  S  and 
constant  /,  as  a  function  of  (F«  — Fj),  and  also 
•S'/,  ( /  )  at  constant  (  F,  —  Fj )  as  a  function  of  the 

screen  grid  voltage  F^;  they  found  horizontal  lines  in  both 
cases.  They  also  deduced^ 


•S,.(/)/[(F. 

from  their  data,  plotted  it  versus  (F.  —  F^)’'^5*/,  and 
found  again  horizontal  lines.  The  average  values  were 

5,_(10)/[(F.  -  F^)^'W.] 

=  1.27x10-“  AHz-' 


for  tube  No.  1  and0.79x  10““  A  Hz"'  for  tube  No. 
2  both  at  to  Hz.  The  average  over  both  tubes  is  therefore 
1.03x  10-“  A  Hz-^'*  at  10  Hz.  Assuming  Eq.  (2)  to  be 
valid  they  found  =  (d^, ),..  =  0.56  cm  from  this  average; 
individual  tubes  give  d^,,  =  0.45  cm  and  dj,j  =  0.73  cm, 
respectively. 

Actually,  Fang  and  van  der  Ziel  did  eight  independent 
measurements  on  each  tube  and  each  measurement  can  yield 
additional  information.  To  retrieve  this  we  introduce  the 
theoretical  a„,  denoted  by  (a„),h«o  and  the  experimental 
a„,  denoted  by  According  to  Eqs.  ( lb)  and  (2) 


)|h«» 


_  522 

3ir  ’ 


(««) 


't%p 


S,.(/) 

(e/.Au) 


2d,J. 


(3) 


TABLE  I.  Device  No.  l./=  10  He. 


(a)  /„  =  lOmA,  5=  3  6. 

K  -  y. 

(V) 

K,  =  250  V 

S,  (lOHz) 

(t0-'“  AVHz) 

(mm) 

125 

2.32 

4  48 

100 

1  56 

4  77 

75 

1  06 

4  56 

50 

0.48 

5  48 

(b)  F,  -  y. 

=  125  V 

y. 

/, 

5,  (10  Hz) 

d,.' 

(V) 

(mA) 

S  (10  '“AVHz) 

(mm) 

250 

10  00 

3  60 

2.30 

4  52 

225 

5  88 

304 

0  94 

4  64 

200 

2.94 

2.93 

0.27 

5  16 

175 

1.44 

1  97 

0  15 

2.99 

*</„.=  4  8  ±  0.3  mm  when  the  Iasi  measurement  is  left  out  of  the  averaging 
procedure. 


We  see  that  is  independent  of  whereas 

(Off  ),,p  varies  linearly  with  d^, .  By  proper  choice  of  d^,  one 
can  thus  bring  the  two  a„  values  to  coincidence;  equating 
=  (ot/y).h«,  yicldsd,j,  for  each  measurement.  Actu¬ 
ally,  due  to  errors  in  the  measurement,  the  individual  values 
of  d^.  will  scatter. 

The  aim  of  this  communication  is  to  point  out  that  the 
averaging  process  eliminates  a  large  amount  of  valuable  in¬ 
formation.  Also,  it  neglects  the  possibility  that  the  secondary 
electron  path  length  between  the  anode  and  dynode  can  be 
different  for  the  two  tubes.  Therefore,  we  calculate  the  value 
ofdj,  from  the  measurements  for  each  tube  separately  using 
Eq.  (3).  For  each  tube  two  runs  of  four  measurements  are 
available;  one  run  for  which  the  anode- to-dy node  voltage 
( F„  —  Fj )  is  varied  between  50  and  1 25  V  while  keeping  the 
grid  voltage  constant  and  one  run  where  the  grid  voltage 
(Fj)  is  varied  between  175  and  250  V  when  the  anode-to- 
dynode  voltage  is  kept  constant.  The  results  are  presented  in 
Table  I.  Similar  experiments  were  done  on  a  second  tube  and 
these  results  are  supplied  in  Table  II.  The  values  of  d^.  tend 


TABLE  II  Device  No  2./=  10  Hz. 


(a)  /,  =  8.0  mA,  S  = 

y.-y. 

(V) 

3  1.  K,  =  250  V 

S,  (lOHz) 

(IO'-aVHz) 

(mm) 

125 

8.7 

709 

100 

6.1 

7  24 

75 

4.3 

667 

50 

2.6 

600 

(b)  y.  -  y. 

=  125  V 

yj 

/. 

S,  (lOHz) 

(V) 

(mA) 

<5 

(10- '’AVHz) 

(mm) 

250 

800 

3  1 

8  7 

709 

225 

50 

2  6 

3.4 

7  98 

200 

2.5 

2  1 

13 

6  81 

175 

1.2 

1.7 

0  3 

9  28 

•rf.i,  =  t  o  ±  0.6  mm  when  the  list  measurement  is  lefi  out  of  the  averaging 
procedure. 
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to  deviate  from  the  average  value  for  small  bias  voltages  Kj 
and  (y,  —  Vj)  or  small  currents  in  both  devices.  Omitting 
these  values  from  the  averaging  procedure  yields  d,. 

=  4,8  +  0.3  mm  for  tube  No.  1  and  =  7.0  ±  0.6  mm  for 
tube  No.  2  in  agreement  with  Fang  and  van  der  Ziel’s  data.  ’ 
Apart  from  the  fact  that  these  are  reasonable  values  for  the 
anode-dynode  distance,  the  root-mean-square  deviation 
from  the  average  is  less  than  8%  for  seven  measurements  on 
each  device.  The  eighth  measurement  in  each  set  deviates 
strongly  from  the  other  seven.  These  are  the  measurements 
at  the  smallest  bias  voltages  and  currents;  they  do  not  deviate 
systematically  cither  to  the  high  or  the  low  side.  Hence  we 
consider  them  poor  measurements. 

We  conclude  that  Eq.  ( 3 )  accurately  describes  the  mea¬ 
sured  noise  in  secondary  emission  tubes.  Not  only  is  the 


magnitude  correctly  predicted  but  the  dependence  on  the 
applied  voltages  and  the  current  as  well 

The  screen-grid-dynode-anode  system  guides  the  pri¬ 
mary  electron  beam  towards  the  dynode  and  the  secondary 
electrons  to  the  anode  The  point  where  the  primary  beam 
hits  the  dynode  will  vary  from  tube  to  tube  due  to  slight 
variations  in  the  screen  geometry,  screen  grid  distance,  grid 
dynode  distance,  and  may  also  depend  somewhat  on  bias 
Accordingly,  the  path  the  electrons  travel  from  dynode  to 
anode  will  vary  from  tube  to  tube,  as  already  demonstrated. 

'R  C  Schwantcs  and  A  van  der  Ziel.  Physica  26,  1 162  f  1*^82 ) 

’A  van  der  Ziel.  Physics  144  B.  205  ( 1987) 

*P  Fang  and  A  van  der  Ziel.  Physica  147  B.  3 1 1  (  1987 ) 

"F  N  Hooge.  Phys  Lett  A-29.  139  ( 1969) 
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GENERATION-RECOMBINATION-TYPE  1//  NOISE  IN 

n-i-p  DIODES 

A.  VAN  DEB  ZiEL,  L.  He.  a.  D  van  Rheenen  and  P.  Fang* 

Electrical  Engineering  Department,  University  of  Minnesota.  Minneapolis,  MN  SS45S,  U.S.A. 

{.Rectwed  13  March  1989;  in  revised  form  6  April  1989) 

Abatract — It  is  shown  for  p-i~n  diodes,  in  which  the  current  flow  is  by  hole-electron  pair  generation  and 
(or)  recombination,  that  the  I//'  noise  is  due  to  generation-recombtnation  processes  involving  traps  and 
(or)  recombination  centers  and  that  the  spectrum  may  be  written  as  5,(/)  >  aHe|/|/l/'r),  where  Sh  is 
the  Hooge  parameter,  e  the  electron  charge,  |/|  the  absolute  current,  r  the  time  constant  associated  with 
the  pair  generation  and  pair  recombination  process, /the  frequency  and  y  is  the  exponent  of  the  spectrum. 

This  is  studied  experimentally  and  the  Hooge  parameters  of  various  devices  are  determined. 


I.  THEORY 


According  to  HoogefI]  the  relative  I//  noise  spec¬ 
trum  of  a  long  uniform  resistor  is  given  by  the  Hooge 
equation 

S,(/)//'-«„/(/\)  (1) 


where  ixh  <8  the  Hooge  parameter,  /  the  frequency, 
and  N  the  number  of  carriers  in  the  system.  It  was 
found  that  Hh  bad  a  value  of  2  x  10*’  for  long 
resistors,  decreased  with  decreasing  device  length(2] 
and  was  many  orders  of  magnitude  smaller  for  small 
Junction  diodes,  bipolar  junction  transistors,  and  field 
effect  transi$tors(3). 

It  should  be  understood,  of  course,  that  in  the 
small  devices  eqn  (I)  must  be  properly  modified.  For 
example,  since  the  devices  are  non-uniform,  eqn  (I) 
must  be  written  in  differential  form.  Moreover,  since 
the  spectrum  is  often  not  exactly  I//.  /  must  be 
replaced  by  /^  The  contribution  of  the  noise  pro¬ 
duced  in  a  slice  of  thickness  At  located  at  x  to  the 
total  spectrum  is 


5,(.t./) 


A.t)^Xh 

pN{x)Ax 


(2) 


where  is  now  not  a  dimensionless  constant  but  has 
the  dimension  (Hz)' ' '  N{x )  is  the  number  of  earners 
per  unit  length  and  /(.x )  is  the  minority  earner  current 
at  X.  This  generalized  formula  is  usually  valid.  For 
long  n*-p  diodes  S,(x,f)  and  /(.x)  decrease  with 
increasing  x. 

The  calculation  of  the  spectrum  of  the  current 
fluctuations  in  the  external  circuit  is  done  by  properly 
integrating  eqn  (2).  For  an  n'-p  H, ,  ,Cd,Te  diode 
in  which  the  current  flow  is  by  diffusion,  this  proce¬ 
dure  yields(4] 


(3) 


where  the  function  F{  V)  follows  from  the  integration 
process.  For  long  diodes  r  is  the  life  time  of  the 
carriers,  which  can  be  evaluated  from  the  h.f.  device 
admittance,  and  for  shorter  diodes  t  =  r ^  =  wJ/2Z),  is 
the  electron  diffusion  time,  which  can  be  evaluated  if 
the  electron  diffusion  constant  D,  and  the  length  w, 
of  the  p-region  are  knosvn.  The  Hooge  parameter  Zh 
can  then  be  determined;  for  long  diodes(4] 
Zh  » (3-S)  X  I0’\  comparable  to  Hooge’s  value 

2  X  10-^ 

For  n-i-p  diodes  we  postulatedfS]  by  analogy 

....  ZH'|/dl/(F) 

S,(f)  = - -pr - .  (4) 


where  t  is  again  the  carrier  lifetime  and  the  function 
/(F)  may  differ  from  the  corresponding  function  in 
eqn  (3).  Experimentally  Zh/CF)  was  found  to  be 
independent  of  back  bias.  Since  Zh  is  supposed  to  be 
a  constant. /(F)  is  apparently  also  a  constant.  As¬ 
suming  /( F)  —  I  we  obtained  very  reasonable  values 
for  a„(51 

It  is  the  aim  of  this  note  to  prove  the  validity  of  eqn 
(4)  and  to  show  that  /(F)  »  I.  We  do  so  as  follows. 
Van  der  Ziel  has  observed  that  in  the  cast  of  1// noise 
due  to  recombination  centers  the  number  S  of  elec¬ 
trons  in  eqn  (I)  must  be  replaced  by  the  number  Vy 
of  recombination  ceniers(6J.  Moreover,  for  the  dis- 
tnbuted  form  corresponding  to  eqn  (2)  l{x)  must  be 
replaced  by  the  net  current  A/(.x)  A.x  generated  in  the 
section  Ax  at  x.  We  then  obtain  for  the  noise 
generated  in  the  section  A.x 


AS,(z,/)  Ax 


[A/(.x)  Ax]^Zh 
/'/Vt(jc)Ax  • 


(5) 


where  Nfix)  A.x  is  the  effective  number  of  recombina¬ 
tion  centers  in  the  section  A.x  at  x.  But,  if  t  is  the 
time  constant  associated  with  the  generation- 
recombination  process 


|A/(.x)|  Ax -eA'T(z)  Ax/r  (6) 

is  entirely  due  to  g-r  processes  where  t  can  be 
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determined  from  the  h.f  device  admittance(5,7)  Sub¬ 
stituting  eqn  (6)  into  eqn  (S)  yields 

A5,(.r./)A.t=^|A/(;t)|A;t  (7) 

/'t 

and  hence  by  integration 

(8) 

Comparing  eqn  (8)  with  eqn  (4)  shows  that  eqn  (4) 
is  correct  and  that  f[V)  »  I. 

Since  the  above  approach  is  derived  by  slicing  the 
device  up  into  sections  A.t,  it  is  generally  valid, 
holding  both  for  volume  and  surface  g-r  processes. 

2.  EXPERIMENTS  ON  H|,..Ca,T*  DIODES 
WITH  t  -a.M  AT  77K 

Noise  measurements  were  performed  on  relatively 
short  Hg|  _  ,Cd,Te  p-i~n  diodes  with  cadmium  mo¬ 
lar  fraction  jr=0.20  at  77  K.  It  was  found  that 
(S,(/)  -  S,(«))]/|/|  was  nearly  independent  of  bias 
for  -  100  mV  <  K  <  0;  here  S,(oo)  is  the  shot  noise 
of  the  devices.  This  is  equivalent  to  what  was  found 
in  long  /i‘p  Hg,  _,Cd,Te  diodes,  where 
(S,(/)  -  S,(ao ))/{//( K)|  was  independent  of  bias. 
Applying  eqn  (8)  we  thus  see  that  Sh/t  is  independent 
of  bias.  For  F  <  —  100  mV  the  noise  increased  due  to 
pre-breakdown  processes;  for  F  >  0  the  noise  also 
increased  due  to  carrier  injection  into  the  i-region. 

The  time  constant  r,  determined  with  the  help  of 
the  diode  admittance  method,  was  found  to  be  I0~’  s 
and  was  nearly  independent  of  bias.  The  resulting 
values  of  3^  are  shown  in  Fig.  I  and  Table  I;  they  are 
also  nearly  independent  of  bias. 

We  now  turn  to  the  diode  admittance.  It  was 
shown(5,7)  that  for  relatively  short  p~i~n  diodes  the 
diode  itself  can  be  represented  by  a  (g,,  C,)  parallel 
circuit,  where  g,  >  d//  dF  and  C,  is  the  diode  capaci¬ 
tance.  In  addition  the  i-region  has  a  series  resistance 
/!,(/);  it  decreases  with  increasing  current  because  of 
carrier  injection  (holes  and  electrons)  into  the  i- 
regton.  If  an  a.c.  voltage  is  applied  to  the  p-i-n  diode, 
the  current  /,  and  hence  /?,(/).  are  a.c.  modulated. 
The  net  effect  is  that  A,(/)  must  be  replaced  by  an 
(R,L)  circuit  (modulation  impedance),  where 


Fig  I  The  Hooge  parameter  as  a  function  of  the  diode 
voltage  (diode  0-1 4) 


Fig.  2.  Conductance  of  Ihe  diode  as  a  function  of  frequency 
at  different  temperatures. 

R  ^  ARJ  d!  and  L  is  proportional  to  the  carrier 
lifetime  r.  As  a  consequence  the  device  shows  a  sharp 
series  resonance  (Fig.  2)  at  the  frequency 
0)0  =  (LC,)"'  ^  from  this  resonance  the  life  time  t  can 
be  determined 

Figures  3  and  4  show  15,(/)  —  S/(ao))  for  diode 
G-14  as  a  function  of  frequency  for  different  bias 
voltages.  The  spectra  are  of  the  form  I//’  with  y 
differing  significantly  from  unity.  Since  quantum  I// 
noise  has  a  y -value  that  is  very  close  to  unity(3],  it 
follows  that  the  noise  is  not  of  quantum  origin. 
Figure  5  shows  y  vs  bias  for  diode  G-14 

Since  in  a  p-i-n  diode  the  current  flow  is  by 
generation-recombination  processes,  the  I//  noise 
should  be  associated  with  these  processes.  Except 
perhaps  in  the  cleanest  samples,  the  noise  involves 
traps  or  recombination  centers  and  S,(/)  should  be 
proportional  to  the  density  of  these  centers.  Since  the 


Table  I 


# 

•"((mV) 

/.OiA) 

5,(/)(A'/Hi) 

1  (r  -  10  ’  5) 

M 

-08 

-3  5 

2  25  V  lO  ” 

2  0  X  10  ’ 

G 

-3  0 

-3  5 

0  95  X  10  ” 

1  5  X  10’ 

F 

-20 

-60 

4  3  X  10  ” 

2  25  X  10  ’ 

J 

-25 

-70 

0  9  X  10” 

80  X  10  * 

K 

-30 

-10  0 

1  0  X  10  ” 

3  0  X  10  * 

H 

-40 

-12.0 

0  8  X  10  " 

2  0  X  10  ‘ 

Fig  3  Noise  spectrum  under  different  reverse  bias;  10,  30 
and  80  mV 
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Fraquwicy  (Hi) 

Fig.  4.  Noise  spectrum  under  different  bias,  2.4,  1 14  mV 
(both  are  reverse  biased),  and  28  mV  (forward  bias). 


density  of  the  centers  will  differ  from  unit  to  unit, 
5/(/)  and  hence  3^  differ  from  unit  to  untt,  as 
experiment  indicates  (Table  I).  Since  the  spectrum 
has  y  #  I.  the  noise  is  not  of  the  quantum  I// 
noise-type  and  is  non-essential  noise.  For  the  I// 
noise  is  proportional  to  the  density  of  the  centers  and 
can  be  reduced  by  reducing  that  density. 

The  limiting  I// noise  is  obtained  when  all  recom¬ 
bination  centers  have  been  removed;  in  that  case  the 
noise  is  due  to  dirtct  band  to  band  transitions  and  has 
a  fundamental  nature.  It  is  a  form  of  quantum  I// 
noise  that  can  be  described  as  follows. 

In  the  generation  process  an  electron  and  a  hole  are 
created  simultaneously,  each  with  an  average  energy 


>0 


1  s 


I  0 


os 

ISO  <00  IS  0  so 

Vo4tag*(mV) 

Fig.  5.  The  exponent  y  as  a  function  of  the  bias  voltage  for 
diode  G-14 


ikT/2.  Therefore  the  electron  mean  square  velocity  is 
vl  =  3kT!m*  and  the  hole  mean  square  velocity 
is  r;*  ikT/m*,  where  m*  and  m*  are  the  effective 
masses  Since  m*  <  m*.  we  have  r;  >  rj  so  that  only 
the  electrons  give  a  significant  contnbution.  Hence(5] 


43,Ai  ;_43o  ikT 
3Tt  c’  3?r  m*c‘ 


(9) 


At  T  =  n  K  and  t  =  0.20  we  have  m*  =  0.0073m,, 
where  m,  =  9.  II  x  10""  kg  is  the  mass  of  a  free 
electron.  Hence  =  I  7  x  10  *.  which  is  a  factor  100 
lower  than  the  lowest  Xh  value  measured  so  far.  It  is 
therefore  doubtful  whether  the  quantum  limit  for 
p-i-n  Hg,  _,Cd,Te  diodes  will  soon  be  reached. 


J.  CONCLUSIONS 

1.  The  Ilf  noise  in  p-i~n  diodes  is  of  the  genera¬ 
tion-recombination  type,  involving  recombination 
centers.  The  noise  is  proportional  to  the  density  of 
these  centers. 

2.  The  best  unit  measured  so  far  has  Sh  =:  2  x  10"*; 
the  quantum  limit  lies  about  a  factor  100  lower. 

3.  A  formula  for  S,{f)  is  derived  that  allows  the 
evaluation  of  the  Hooge  parameter  3h  from  the 
measured  noise. 
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BURST-TYPE  NOISE  MECHANISMS  IN 
BIPOLAR  TRANSISTORS 

X.  L.  Wu.t  A.  VAN  OER  ZiEL,  A.  N.  Birbas  and  A.  D  van  Rheenen 
Department  of  Electrical  Engineering.  University  of  Minnesota.  Minneapolis,  MN  S54SS,  U.S.A. 

(Rtetived  27  March  1989) 

Abatract — Two  types  of  burst  noise  have  been  observed  tn  silicon  bipolar  transistors.  They  can  be 
characterized  by  the  typical  frequency  dependence  of  their  current  fluctuation  spectra.  Interestingly, 
observation  of  the  noise  signal  in  the  time  domain  gives  two  distinctively  diflTerent  pictures  of  the  bistable 
waveform.  Also,  amplification  of  the  noise  signal  yields  diflierent  sounds  when  fed  to  a  speaker.  One  of 
the  noise  spectra  is  the  superposition  of  I// noise  and  a  Lorentzian  component  (burst  noise)  and  the  other 
can  be  described  mathematically  as  I  //  noise  modulated  by  the  burst  noise  The  classification  of  those  two 
types  of  burst  noise  and  the  mathematical  explanation  will  lead  to  a  better  understanding  of  the  bipolar 
transistor  bunt  noise  itself 


I.  INTROOUCnON 

Bistable  current  waveforms  resembling  a  random 
telegraph  signal  known  as  burst  noise  have  been 
reported  for  the  last  30  years.  This  phenomenon  has 
been  observed  in  forward  and  reverse  biased  diodes, 
bipolar  junction  transistors,  optical  isolators  and 
detectors,  junction  field  effect  transistors  and 
metal-oxide-semiconductor  field  effect  transistors. 

The  origin  of  the  burst  noise  has  been  attributed  to 
physical  mechanisms  such  as  microplasmas  related  to 
the  dislocations  in  the  crystal,  emitter  edge  disloca¬ 
tions  of  bipolar  transistors,  crystallographic  defects, 
generation-recombination  (g-r)  centers  associated 
with  metallic  precipitates  or  extensive  defects  near  the 
surface  of  the  junction,  etc.  It  has  also  been  estab¬ 
lished  that  the  burst  noise  in  JFETs  is  related  to  point 
defects  in  silicon  and  the  burst  noise  in  MOSFETs 
has  been  attributed  to  charge  state  transitions  at  traps 
in  the  oxide  very  close  to  the  Si/SiOj  interface. 

The  mechanism  that  causes  the  burst  noise  in 
bipolar  transistors  is  not  well  understood  even 
though  much  work  has  been  done  in  the  area.  This 
is  due  to  the  poor  understanding  of  the  mechanisms 
for  current  generation  and  the  spatial  nonuniformi¬ 
ties  in  the  current  flow.  However,  it  is  well  established 
that  dislocations  at  the  emitter  as  well  as  trapping  are 
sources  of  bunt  noise  and  I  //  noise.  It  is  the  aim  of 
this  paper  to  shine  some  light  on  the  burst  noise 
originated  at  the  emitter  junction  of  a  bipolar  transis¬ 
tor  as  well  as  its  correlation  with  the  locally  produced 
I//  noise. 

L  MECHANISMS  OF  BURST  NOISE 

Much  work  has  been  done  in  the  area  of  low- 
frequency  noise  of  silicon  planar  devices.  Various 


tPresent  address;  Unisys  Corp.,  San  Diego,  CA  92127, 
USA. 


authors  have  attributed  the  low-frequency  noise  to 
dislocations  at  the  emitter-base  junction.  Morrisonjl] 
predicted  that  it  could  be  caused  by  the  dislocations 
in  the  bulk  whereas  Bess{2]  attributed  the  noise  to 
dislocations  that  communicated  with  the  silicon  sur¬ 
face  and  thus  allowed  impurity  atoms  to  diffuse  to  the 
silicon  surface.  Green{3]  and  Jordan  have  established 
that  the  low-frequency  noise  is  affected  by  ute  dislo¬ 
cations  only  when  they  emerge  at  a  silicon  surface 
near  the  junction. 

Emitter-phosphorus  diffusion  normally  used  for 
fabrication  of  silicon  planar  n-p-n  transistors  may 
create  two  kinds  of  dislocations;  (i)  diffusion-induced 
dislocations,  which  appear  inside  diffused  emitter 
areas,  and  (ii)  emitter-edge  dislocations,  which  ap¬ 
pear  around  the  planar  edges  of  the  diffused  emitter 
areas,  usually  after  the  emitter  oxidation  step  While 
researchers  in  the  past  have  attributed  the  low- 
frequency  noise  of  silicon  planar  n-p-n  transistors  to 
both  diffusion-induced  dislocations  and  emitter-edge 
dislocations!#],  studies  of  modem  planar  devices 
show  that  the  low  frequency  noise  in  n-p-n  transis¬ 
tors  is  affected  only  by  emitter-edge  dislocationsfS). 
The  effect  is  enhanced  rapidly  in  n  *  emitters  for  emit¬ 
ter  doping  densities  larger  than  4.3  x  10“  cm  ’(6]. 

We  can  now  discriminate  between  various  mecha¬ 
nisms  responsible  for  the  low-frequency  noise.  If  the 
emitter  doping  exceeds  the  4.3  x  I0“cm‘’  value, 
then  there  are  a  number  of  edge  dislocations.  It  is  also 
well  known  that  considering  the  regularity  of  the 
pulse  amplitude  (about  10'*  A)  associated  with  the 
burst  noise,  it  is  inconceivable  that  we  are  dealing 
with  the  statistics  of  a  large  number  of  carriers.  It 
seems  more  feasible  that  a  single  carrier  controls  the 
flow  of  10*  carriers  in  a  reasonably  periodic  fashion 
producing  a  pulse  duration  about  I  ms.  Indeed,  a 
defect  associated  with  the  edge  of  the  emit¬ 
ter-metallurgical  junction  can  produce  burst  noise 
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due  to  the  modulation  of  the  current.  The  modulation 
is  caused  by  the  change  in  occupancy  of  a  genera¬ 
tion-recombination  center  located  in  the  space- 
charge  region  adjacent  to  the  defect(7].  Furthermore, 
when  a  lot  of  dislocations  are  present  trapping  of 
carriers  can  occur  at  that  point.  If  there  is  a  uniform 
distribution  of  the  time  constants  associated  with  the 
trapping  and  detrapping  then  the  resulting  noise 
spectrum  is  proportional  to  I//.  In  conclusion  we  see 
that  in  tiie  case  of  a  heavily  doped  emitter  the 
dislocations  produce  low-frequency  noise.  The  trap¬ 
ping  effect  of  one  of  those  dislocations  (the  closest  to 
the  metallurgical  junction  with  the  appropriate  time 
constant)  may  give  burst  noise  and  the  superposition 
of  the  others  may  be  giving  I//  noise  if  there  is  a 
uniform  distribution  of  their  time  constants.  Both  the 
burst  and  the  I//  type  of  noise  are  due  to  the 
dislocations,  they  originate  from  the  same  region 
of  the  device  and  since  they  are  due  to  the  same 
physical  mechanism  they  are  related  and  the  one  type 
modulates  the  other. 

If  dislocations  are  not  present  at  the  emitter  edge 
then  the  burst  noise  can  be  produced  at  the  space- 
charge  region  of  the  device  if  a  trap  is  presentjTj. 
Then  the  burst  modulation  current  will  be  given  by 
^lf^iq^ll2ntkT(y/itA)iT\  where  all  the  quantities 
have  the  normal  meaning,  and  A  is  the  effective  cross 
sectional  area  at  a  forward  bias  junction  with  a 
boundary  at  the  surface  (planar  BJT).  In  addition  to 
the  burst  noise  in  that  case  we  can  have  Mf  noise 
associated  with  the  emitter-base  junction  recombina¬ 
tion  current  which  is  not  transmitted  to  the  collector. 
In  modem  transistors  we  rather  have  diffusion  Mf 
noise  and  in  good  low  noise  devices  mobility  fluctua¬ 
tion  Mf  noise(8).  In  this  case  the  burst  noise  and  the 
I //noise  come  from  the  same  region  but  they  are  due 
to  different  mechanisms;  so  they  are  unrelated,  act  in 
parallel  and  their  spectra  are  superimposed. 


The  rate  at  which  open  switches  will  close  is  given  by: 
dt  r. 


[iV,  -  A.V,  1  -  4- W(M  (1) 

t; 

where  H{t)  is  a  random  source  function  describing 
spontaneous  fluctuations  in  the  rate  at  which  switches 
open  or  close. 

Because  of  the  constraint  on  ;V,  and 
+  ^,  =  ,V  =  total  number  of  switches)  we 
see  that  AiV,  =  -AiV, .  Additionally  d;V|;dr  =  0  by 
definition.  Therefore  eqn  (1)  reduces  to; 

+  r-=r,-'  +  r;-  (2) 


A?i 


from  which  we  find 


N,  =  N- 


R,  =  M- 


'J 

(ti  +  tj)' 

since  we  have  a  binomial  distribution. 

Once  we  find  the  rate  equation  (eqn  2).  the  spectral 
intensity  of  the  fluctuations  in  the  number  of  open 
switches  is  computed  by  Fourier  transforming  eqn  (2) 
and  averaging.  The  spectrum  will  be 


In  this  equation  5«  is  the  spectrum  of  the  driving 
force  HU). 


3.  SWITCHING  MODEL  OF  THE  BURST  NOISE  *  BURST  NOISE  AND  Ilf  NOISE 


We  now  turn  to  the  formulation  of  the  two  cases 
of  the  joint  appearance  of  1  If  and  burst  noise  mech¬ 
anisms  described  in  Section  2.  Burst  noise  is  a  form 
of  g-r  noise.  It  consists  typically  of  random  current 
pulses  of  variable  duration  and  equal  height,  but 
sometimes  the  random  pulses  seem  to  be  supenm- 
posed  upon  each  other  The  random  switch  model[8j 
considers  that  the  burst  noise  is  caused  by  an  assem¬ 
bly  of  random  switches,  each  carrying  a  fixed  current 
when  closed  and  no  current  when  open.  During  the 
time  interval  di  let  an  open  switch  have  the  probabil¬ 
ity  di/r,  of  closing,  and  a  closed  switch  a  probability 
di/t,  of  opening.  Let  furthermore  HU)  be  a  random 
source  function  describing  spontaneous  fluctuations 
in  the  rate  at  which  switches  open  or  close. 

At  a  given  time  t  S,  switches  are  open  and  (Vj 
switches  are  closed  Develop  around  the  equilibrium 
values  we  find:  iV,  »  jV,  +  AAj  and  jV, «  +  A/V,. 


(a)  In  Section  J  eqn  (I)  describes  the  burst  noise 
produced  by  the  rate  fluctuations 

Let's  now  turn  to  the  application  of  this  equation 
to  an  electronic  device  where  carrier  number  fluctua¬ 
tions  (I// noise)  appear  along  with  the  burst  noise. 
We  consider  that  an  Mf  modulation  occurs  of  the 
number  of  carriers  (.V,  open  switches)  along  with  the 
burst  noise  modulation  H^(t)  and  let  this  be 
described  by  a  source  function  HfU)  Then  eqn  (I)  is 
modified  to: 


d(A;V,  +  ^,)  „  I 

- - -  =  -  [iV,  +  A;V|  +  H,u  )1  — 

df  Ti 


+  [,*7,-AA,-W, (I )]--)- //,(!)  (4) 
T; 


The  d.c.  part  of  this  equation  is  the  same  as  that 
descnbed  in  eqn  (2).  Then  the  time  dependent  part  of 
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eqn  (4)  may  be  written 


dAA^i 

dr 


A/V, 

To 


To 


«o{0; 


t 


-I 

0 


Tf'+Tf'  (5) 


Keeping  in  mind  that  //,(r)  and  H^(t)  are  indepen¬ 
dent  and  after  Fourier  transforming  we  get; 


SH(0)tl  +  S,if) 
I  -hro^ro 


Axl  +  B  If 
I  <U*Tj 


(6) 


since  S„(0)  =«  A  and  Bjf. 

■rVie  resulted  spectrum  is  a  modulation  of  the 
Lorentzian  (burst  noise)  by  the  I.//  noise. 


(h)  I  If  noise  and  burst  (Lorenlzian)  noise  occur  in 
parallel  and  have  independent  spectra  because  they 
come  from  fluctuations  in  independent  parameters 

The  time  dependent  Langevin  Equation  for  the 
burst  noise  is  then: 


dAAf, 

dr 


(7) 


where;  is  the  burst  noise  source;  after  Fourier 
transforming  and  averaging  we  get  the  first  term  of 
the  right  hand  side  of  eqn  (3)  where  we  replace  by 
Sht,  Furthermore,  if  we  consider  that  a  1//  noise 
source  which  corresponds  to  an  independent  process 
occurring  either  in  the  same  or  in  a  different  region 
of  the  device  then  a  I //source  H,  has  to  be  superim¬ 
posed  on  the  burst  noise.  This  I//  response 
{•S*<(/)  “  BIf)  is  independent  of  the  burst  noise  and 
will  give  the  second  term  of  the  right  hand  side  of 
eqn  (8). 

The  total  spectrum  will  be; 


S(f) 


A^t^  B 
I  /  ’ 


(8) 


In  this  case  the  resulting  spectrum  is  a  superposition 
of  the  burst  noise  and  the  I// noise. 


1  EXPCRIMCNT  AND  DISCt'SSION 

Measurements  of  the  electric  noise  m  two  different 
devices  illustrate  nicely  the  different  spectra.  The 
spectrum  of  the  collector-current  of  a  discrete  silicon 
BJT  IS  shown  in  Fig.  I.  The  solid  line  in  this  figure 
is  the  best  fit  using  eqn  (6).  At  the  high  frequency  end 
the  spectrum  reaches  the  thermal  noise  limit  of 
3.8  X  I0‘'*A^/Hz.  The  most  interesting  feature  of 
this  spectrum  is  that  it  decreases  as at  larger  fre¬ 
quencies  until  it  reaches  the  thermal  noise  floor  The 
effect  is  more  dramatic  when  one  subtracts  the  ther¬ 
mal  noise  floor  from  the  data.  This  is  done  for  the  last 
6  data  points  and  indicated  by  the  open  squares  in 
Fig.  I .  The /*'  dependence  is  even  more  pronounced. 
This  observation  is  in  accordance  with  the  expression 
in  eqn  (6)  At  higher  frequencies  the  right  hand  side 
of  eqn  (6)  is  dominated  by  the  first  term  and  falls  off 
asf-r 

Ht  1MI-H 


Fig.  I.  The  experimental  spectrum  (•)  of  the  discrete 
silicon  BJT  and  the  theoretical  spectrum  following  eqn  (6) 
( — )  The  open  squares  (□)  are  experimental  data  points 
with  the  thermal  noise  floor  subtracted.  At  the  high  fre¬ 
quency  end  the  spectrum  falls  off  faster  than  /  ”'. 


As  an  example  of  the  second  type  of  burst  noise 
(Section  4b)  we  present  the  fluctuation  spectrum  of 
the  base  current  of  a  BJT  on  a  chip,  in  Fig.  2.  The 
spectrum  behaves  as  predicted  by  eqn  (8).  By  inspect¬ 
ing  eqn  (8)  we  find  that  at  the  higher  frequencies  it 
should  fall  off  as  /“ '  We  observe  this  phenomenon 
in  Fig.  2. 

When  comparing  the  two  equations  (6  and  8)  there 
is  no  difference  between  the  spectra  at  the  low-fre¬ 
quency  end  (provided  that  we  choose  the  same  values 
for  the  constants  A,  B  and  t„).  However,  the  differ¬ 
ence  between  the  two  models  can  be  observed  at  the 
high-frequency  end  of  the  spectra.  Spectra  following 


Fig.  2  The  experimental  spectrum  (#)  of  the  BJT  on  the 
chip  and  the  theoretical  curve  following  eqn  (8)  At  the  high 
frequeiKy  end  the  spectrum  falls  off  as  / 
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Fig.  3.  The  oscilloscope  picture  of  the  burst  noise  of  type 
a.  The  duration  of  the  total  time  signal  is  40  ms. 


Fig.  4.  The  oscilloscope  picture  of  the  burst  noise  type  b. 
The  duration  of  the  total  lime  signal  is  40  ms. 


eqn  (6)  will  fall  off  as  and  spectra  following 
eqn  (8)  will  fall  off  as  /"'. 

In  the  time  domain  the  two  burst  noise  types  have 
two  distinct  patterns  as  is  shown  in  Figs  3  and 
4.  The  noise  signal  in  the  time  domain  (Fig.  3)  shows 
two  levels  with  noise  in  each  level.  The  duration  of  the 
total  time  trace  is  40  ms.  The  duration  of  the  up  pulse 
is  of  the  same  order  as  the  one  of  the  down  pulse.  This 
is  in  contrast  with  the  time  signal  of  the  chip  transis¬ 
tor  which  we  depict  in  Fig  4.  On  this  photograph  we 
see  a  broad  noise  base  line  and  short  pulses  on  top 
of  it.  When  this  signal  was  amplified  and  fed  to  a 
speaker  a  distinctive  popping  sound  could  be  heard. 
This  phenomenon  showed  itself  only  with  the  chip 
transistor  and  not  with  the  discrete  one. 

The  discrete  device  has  a  heavily  doped  emitter  and 
we  expect  a  lot  of  dislocations  in  the  emitter-base 
space-charge  region.  Therefore,  we  were  not  sur¬ 
prised  to  see  the  noise  spectrum  to  follow  eqn  (6)  (see 
Section  2),  apparently  the  chip  transistor  did  not 
appear  to  have  too  many  dislocations  and  therefore 
shows  a  different  noise  spectrum,  more  in  line  with 
the  predictions  of  eqn  (8). 


4.  CONCLL'SION 

The  observation  of  two  distinctly  different  burst 
noise  spectra,  mixed  with  I//  noise,  as  measured  in 
two  different  devices  is  explained  by  assuming  eithei 
a  modulation  of  the  1//  noise  by  the  burst  noise 
indicating  a  common  source  for  the  fluctuations  or  a 
superposition  of  the  two  noise  types  indicating  unre¬ 
lated  sources  for  the  fluctuations. 

Acknowledaement — The  auttu'rs  would  like  to  thank  M. 
Mihaila  for  providing  one  of  the  devices. 

REFERENCES 

1.  S.  R  Morrison.  Phvs  Rev.  104,  619  (1956). 

2.  L.  Bess.  Phys.  Rev  m.  72  (1956). 

3.  D  Green  and  A.  G  Jordan,  Inl  J  Electron  27,  159 
(1969) 

4  Ntshida.  IEEE  Trans  ED-20,  221  (1973). 

5.  N  D  Slojadinovic.  Electron  Lett.  15,  340  (1979) 

6.  M.  Mihaila  and  K.  Ambenadis.  Solid-St.  Electron  26, 
109  (1983) 

7  S.  T  Hsu.  R  J  Whiltier  and  C.  A.  Mead,  Solid-St. 
Electron.  13,  1055  (1970). 

8.  A.  van  der  Ziel.  Soise  in  Solid  State  Devices  and 
Circu.is  Wiley  Inlersaence.  New  York  (1986) 

9.  X.  L  Wu.  Ph.D  thesis.  University  of  Minnesota  (1987). 


BIO 


1// noise  in  double-heterojunction  AIGaAs/GaAs  laser  diodes  on  GaAs 
and  on  Si  substrates 

R.  Z.  Fang,*’  A.  D.  van  Rheenen,  A.  van  der  Ziel,  and  A.  C.  Young 

Electrical  Engineering  Department.  University  of  .Minnesota,  Minneapolis.  Minnesota  55455 

J  P.  van  der  Ziel 

A  TA  T  Bell  Laboratories,  Murray  Hill.  New  Jersey  07974-2070 
(Received  15  September  1989;  accepted  for  publication  6  June  1990) 

Low-frequency  electrical  current  noise  measurements  are  reported  on  double-heterojunciion 
AlGaAs/GaAs  laser  diodes  fabricated  on  GaAs  and  on  Si  substrates.  The  noise  spectra  show  a 
frequency  dependence  proportional  with /  '  with  y  close  to  unity.  The  spectral  intensity  is 
proportional  to  the  current  for  smaller  currents  (  <  0  1  mA )  and  levels  off  at  larger  currents 
(  >  I  mA).  The  diodes  built  on  the  GaAs  substrate  are  50  times  less  noisy  than  the  ones  built 
on  the  Si  substrate.  This  effect  is  attributed  to  the  fact  that  the  density  of  dislocations  at  the  Si 
interface  is  much  larger  than  at  the  GaAs  substrate/device  interface. 
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I.  INTRODUCTION 

The  data  rates  needed  in  communications  systems  have 
increased  to  the  point  where  the  only  viable  links  are  optical. 
The  interest  in  these  systems  has  grown  enormously  in  the 
last  several  years.  These  links  consist  of  semiconductor  laser 
with  associated  driver  and  perhaps  multiplexing  and  demul¬ 
tiplexing  circuitry,  optical  fibers,  and  semiconductor  detec¬ 
tors  with  their  associated  amplifier  and  possibly  multiplex¬ 
ing  or  demultiplexing  circuitry.  The  level  of  noise  generated 
in  the  components  is  crucially  important  to  know  particular¬ 
ly  in  long-distance  systems  where  signal  levels  can  be  low 
and  distance  between  repeaters  is  determined  by  the  noise. 

With  the  emergence  of  lasers,  and  semiconductor  laser 
diodes  in  particular,  as  useful  light  sources  there  has  been  an 
enormous  interest  in  the  description  of  the  operation  of  the 
lasers.  The  noise,  specifically  in  the  light  output,  has  been  the 
subject  of  a  large  body  of  publications  that  appeared  starting 
in  the  1960’s.  Different  formalisms  to  describe  this  quantum 
noise  were  developed.  The  two  basic  methods  are  the  Lange- 
vin  method  and  the  density  matrix  method.  Exponents  of  the 
former  are  quantum-mechanical  rate  equations  with  fluctu¬ 
ating  terms,'  ^  the  Fokker-Planck  equation  for  the  photon 
amplitude  probability  density, '  and  the  van  der  Pol  equation 
with  a  noise  source  as  the  driving  force.^  ’  The  density  ma¬ 
trix  method  has  been  studied  by  Scully  and  Lamb*  and  Lax’ 
who  established  that  both  formulations  yield  the  same  result 
when  the  photon  number  (i.e.,  pump  rate)  is  large  enough. 
Corrections  for  small  photon  numbers  were  discussed  by 
Smith.* 

Because  the  Q  factor  is  relatively  low  in  a  semiconductor 
laser  as  compared  to  gas  lasers  and  solid-state  lasers  the 
quantum  noise  will  be  of  enormous  importance  in  such  ap¬ 
plications  where  semiconductor  lasers  are  an  integral  part  of 
optical  fiber  communication  systems.  Components  of  these 
systems  include  transmitters,  moiiulators,  optical  ampli¬ 
fiers,  local  oscillators,  and  detectors.  Yariv'’  identifies  three 
noise  quantities,  the  field  spectrum,  the  intensity  spectrum, 
and  the  frequency  deviation  spectrum  (phase  noise),  that 
describe  the  fluctuations  in  the  laser  oscillation  field. 


Almost  all  of  these  studies  deal  with  the  laser  field  and 
hardly  any  work  is  done  to  relate  the  fluctuations  in  the 
electrical  current  to  fluctuations  in  the  laser  output.  Yama¬ 
moto'"  made  the  observation  that  the  spectrum  of  the  fre¬ 
quency  deviations  seems  to  be  limited  by  the  electrical  cur¬ 
rent  noise  spectrum  at  small  frequencies. 

For  purposes  of  easy  integration  of  silicon  electronics 
with  optical  devices  people  are  investigating  the  possibilities 
to  fabricate  GaAs  base  structures  on  Si  substrates.  With  this 
in  mind  we  set  out  to  study  the  electrical  noise  properties  of 
identical  laser  diode  configurations  composed  on  Si  and 
GaAs  substrates. 

II.  EXPERIMENTS  AND  RESULTS 

The  device  structure  we  investigated  is  presented  in  Fig. 
1  where  we  indicate  the  layer  configuration  and  layer 
thicknesses.  These  double-heterojunction  injection 
AlGaAs/GaAs  devices  are  designed  to  have  a  low  threshold 
voltage,  a  high  quantum  efficiency,  and  to  operate  at  room 
temperature.  Diodes  No.  1  and  No  3  were  grown  on  silicon 
substrates  by  melalorganic  chemical  vapor  deposition 
whereas  diodes  No,  5,  No.  6,  and  No.  7  were  fabricated  on 
gallium-arsenide  substrates  using  the  same  growth  tech¬ 
nique. 

The  current-voltage  characteristics  of  ail  the  diodes  are 
shown  in  Fig.  2.  All  of  these  curves  can  be  represented  by 

/=/,  exp[(K-//?J/nK,  ].  (1) 

In  Eq.  ( 1 )  /  is  the  current,  /,  the  reverse  saturation  current, 
y  the  voltage,  /?,  the  diode  series  resistor,  n  the  ideality  fac¬ 
tor,  and  Vj  the  thermal  voltage  (k„T/q).  Table  1  gives  a 
listing  of  the  values  of  the  ideality  factor  and  the  series  resis¬ 
tance  of  each  diode 

During  the  electrical  noise  measurements  the  diodes 
were  connected  in  series  with  a  lO-kfl  resistor  which  allowed 
us  to  measure  the  current  noise  produced  by  the  lasers.  The 
signal  is  fed  into  a  low-noise  amplifier  and  passed  on  to  a 
dynamic  signal  analyzer  The  measurement  technique  is  de¬ 
signed  to  allow  for  the  noise  produced  by  the  amplifier  and 
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FIG.  i  Structure  of  the  laser  diodes. 


possible  nonflatness  of  the  gain  of  the  amplifier  over  the  fre¬ 
quency  band.  The  current  noise  spectra  are  measured  in  the 
range  from  10  Hz  to  100  kHz. 

T  ypical  examples  of  the  spectral  intensity  of  the  current 
fluctuations  of  both  types  of  diodes,  diode  No.  I  on  a  Si 
substrate  and  diode  No.  5  on  a  GaAs  substrate,  are  presented 
in  Fig.  3.  The  noise  spectra  arc  measured  when  the  currents 
are  as  indicated  in  the  figure.  Even  when  we  measured  at 
frequencies  up  to  100  kHz  did  we  not  see  any  significant 


TABLE  I.  Collected  values  of  the  senes  resistance  (S, )  and  ideality  factor 
(/I)  of  the  diodes. 


Diode 

R 

(ID 

n 

1 

II 

20 

3 

21 

20 

5 

III 

}  5 

6 

31 

20 

7 

14 

20 

leveling  off  of  the  spectra.  This  indicates  that  the  shot  noise 
associated  with  the  carriers  crossing  the />-«  junction  and  the 
thermal  noise  associated  with  the  series  resistance  /? ,  of  the 
diodes  are  well  below  the  measured  excess  noise  in  this  fre¬ 
quency  range.  One  would  exptect  the  shot  noise  and  thermal 
noise  contributions  to  add  accorv  ’n^  ’o 


In  this  equation  is  the  small  signal  resistance  of  the  diode: 


dV  nVr 
dl  ~  I 


(3) 


For  small  currents  (R^  is  large )  S,  =  2ql  whereas  for  large 
currents  [R^  is  small)  S,  =  ‘^kg  T/R,.  The  measured  excess 
noise  was  much  larger  than  this  white  noise. 

The  low-frequency  parts  of  the  spectra  exhibited  a  fre¬ 
quency  dependence  of/  “  '  with  j'close  to  1,  within  10%.  In 
Fig.  4  we  present  the  value  of  the  spectral  intensity  extrapo¬ 
lated  to/=  1  Hz  divided  by  the  current  for  both  laser  diodes 
No.  I  and  No.  5.  Two  interesting  observations  can  be  made; 
(i)  the  value  of  5,  ( 1  Hz)//  for  the  laser  built  on  the  GaAs 
substrate  (diode  No.  5)  is  much  smaller  (factor  of  50)  than 
the  value  of  5,  ( 1  Hz )  //  for  the  diode  built  on  the  Si  substrate 
over  the  whole  range  of  currents  ( 1  /iA-lO  mA);  and  (ii) 
the  relative  noise  is  constant  up  to  ~  I  mA  and  then  starts  to 
decrease  for  both  diodes. 


FIG  2  Current-voltage  chancteriMics  of  the  diodes 


III.  DISCUSSION 

In  order  to  be  able  to  compare  the  magnitude  of  the  1// 
noise  produced  by  different  types  of  devices  the  following 
expression"  is  often  used: 

S,(/)  =a„/V/7V.  (5) 

In  this  equation  N  is  the  number  of  carriers  in  the  device  and 
a„  is  a  dimensionless  parameter  that  measures  the  noisiness 
of  the  device.  This  parameter  is  defined  by  Eq.  (5).  In  Ref. 
12  it  is  shown  that  this  expression  can  be  rewritten  as 

S,{f)=a„qI/fT,  (6) 

where  it  is  assumed  that  the  noise  is  produced  in  the  vicinity 
of  the  lasing  medium.  In  this  expression  the  1// noise  is  as¬ 
sumed  to  be  due  to  a  distribution  of  recombination  centers 
and  T  is  the  lifetime  of  the  minority  carriers.  Equation  (6) 
predicts  the  current  dependence  that  we  found  for  the  low- 
frequency  noise  in  our  diodes  (see  Fig.  4),  at  least  for  the 


J  Appi  Phys  ,  Voi  68.  No  8, 15  October  1990 


4088 


Fartg  etal 


4088 


FIG  3  Spectral  intensity  of  the  current  fluciiiniinns  of  ;i 
typical  diode  on  a  Si  substrate  and  of  a  typical  duKle  on  a 
GaAs  substrate  for  different  current 


10'  102  103  10^ 

Frequency  [Hz] 


smaller  currents  (  <0. 1  mA).  Another  way  oflooking  at  it  is 
to  say  that  a„/T  (see  Eq.  (6)  ]  is  independent  of  the  current 
for  small  currents  and  decreases  for  larger  currents  (  >  I 
m  A ) .  If  anything,  one  would  expect  the  lifetime  of  the  carri¬ 
ers  to  be  reduced  when  more  carriers  are  injected  as  the  cur¬ 
rent  is  increased.  This  would  increase  the  ratio  a„/r.  There¬ 
fore.  we  must  conclude  that  a„  is  decreasing  as  the  current 
increases.  A  similar  effect  was  observed  in  Si  n  -yn  *  and 
p  *  -iT-p  *  structures." 

More  significant  is  that  our  data  indicate  that  diodes 
built  on  GaAs  substrates  are  50  times  quieter  than  identical 
structures  built  on  Si  substrates  It  is  known  that  the  disloca¬ 
tion  density  at  the  Si/GaAs  interface  is  much  larger  than  at 
the  GaAs/GaAs  interface.  Since  the  current  is  flowing  verti¬ 


cally  through  the  device  (Fig.  I)  it  is  not  surprising  that 
generation-recombination  processes  show  up  in  the  noise. 
However,  this  observation  cannot  be  quantified  in  terms  of 
Eq.  (6)  since  the  expression  used  for  the  current  /  =  qS /t  '\% 
applicable  to  the  diode  region,  not  to  the  device-substrate 
interface  region 

IV.  SUMMARY 

We  have  reported  here  electrical  current  noise  measure¬ 
ments  that  were  taken  on  identical  double-heterojunction, 
low-threshold  voltage,  high-quantum-efficiency 
AIGaAs/GaAs  laser  diodes  fabricated  on  Si  and  GaAs  sub¬ 
strates. 


FIG  4  F.\lrnp(ilalcd  valuenrthe  iKiise  spectra  ul  I  ll/diMdcd 
by  the  current  a>i  .1  function  of  the  current  both  for  a  diode  on  a 
Si  substrate  and  for  a  diixle  on  a  GaAs  substrate 
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The  main  observations  are  as  follows:  (i)  The  spectral 
intensity  of  the  current  fluctuations  scale  with  the  current 
for  smaller  (  <  0. 1  mA )  values  of  the  current  whereas  it  lev¬ 
els  olf  for  larger  currents.  This  is  true  for  devices  built  on 
both  substrates,  (ii)  The  diodes  built  on  the  Si  substrate  are 
50  times  noisier  than  the  ones  built  on  GaAs  substrates  in  the 
range  of  measured  currents  ( 1  /zA-lO  mA). 
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Low-frequency  noise  measurements  are  reported  on  small  InGaAs/InP  p-i-n  photodiodes 
under  different  bias  and  illumination  conditions.  In  the  first  experiment  measurements  were 
taken  when  the  diodes  were  reverse  biased  and  illuminated  with  incandescent  light.  At 
high  frequencies  the  noise  is  shot  noise  and  at  low-frequencies  the  noise  has  a  spectrum 
proportional  to/’'  with  y  =  —  1.0  for  one  set  of  two  devices  and  y  =  —  0.8  for  another  set 
of  two  devices.  At  low-frequencies  the  spectral  noise  intensity  is  proportional  to  the 
current  squared.  In  a  second  experiment  the  diodes  were  forward  biased  (no  illumination) 
and  the  spectral  intensity  of  the  low-frequency  noise  was  proportional  to  the  current. 

Under  these  bias  conditions  it  was  possible  to  extract  the  parameter  The  obtained  values 
for  this  parameter  are  not  compatible  with  quantum  1//  noise  but  do  seem  to  coincide 
with  values  related  to  1//  noise  due  to  recombination  centers. 


I.  INTRODUCTION 

In  this  paper  we  report  electrical  noise  measurements 
on  small  (7S-^m-diam)  p-i-n  InGaAs/InP  photodiodes 
designed  to  be  used  in  optical  communication  systems.  In 
certain  applications  the  sensitivity  of  the  photodetectors 
might  be  limited  by  the  noise  and  it  is  therefore  of  interest 
to  quantify  the  noisiness  of  the  detectors  to  be  able  to 
compare  the  noise  performance  of  different  device  struc¬ 
tures.  The  magnitude  of  the  noise  might  for  instance  de¬ 
termine  the  number  of  repeater  stations  necessary  to  send  a 
signal  a  certain  distance.  In  addition,  measurements  of  the 
noise  of  devices  with  different  structures  or  as  a  function  of 
the  bias  will  reveal  the  preferred  device  configuration  or 
the  optimum  operating  conditions.  With  this  in  mind  we 
have  set  out  to  investigate  the  bias  dependence  of  the  cur¬ 
rent  noise  of  two  batches  of  two  different  devices. 

II.  EXPERIMENTS  AND  RESULTS 
A.  Illuminated,  reverse  biased 

In  this  first  set  of  experiments  the  diodes  were  sub¬ 
jected  to  reverse-bias  voltages  of  —  0.5  and  —  3  V.  The 
diodes  were  illuminated  through  the  InP  substrate  with 
incandescent  light.  The  magnitude  of  the  current  (I-I4 
fiA)  was  controlled  by  the  amount  of  light  we  shined  on 
the  diode.  To  be  able  to  measure  the  voltage  noise  due  to 
the  current  fluctuations  in  the  diode  a  resistor  Rg  ( 10  kfl) 
was  connected  in  series  with  the  diode.  This  voltage  signal 
is  fed  into  a  low-noise  amplifier  and  passed  on  to  a  spec¬ 
trum  analyzer.  When  calculating  the  spectra  allowance 
was  made  for  the  noise  produced  by  the  amplifier. 

In  Fig.  1  we  present  the  current  noise  spectra  of  the 
illuminated  diode  biased  at  y  =  —  3  V  The  spectra 
were  measured  for  different  currents:  I,  3,  5,  7,  and  lOfiA. 
At  low  frequencies  the  spectra  are  proportional  to  /’'  with 
y  very  close  (within  6%)  to  —  1  The  magnitude  of  this 
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component  of  the  noise  is  proportional  to  the  current 
squared.  These  results  are  summarized  in  Table  1.. 

At  higher  frequencies  the  spectra  are  independent  of 
the  frequency  and  are  very  well  described  by 

S,(oo)=2ql  -k-^ksT/RB-  (1) 

In  this  equation  5/(  oo )  is  the  spectral  intensity  at  high 
frequencies,  q  the  electronic  charge,  I  the  current,  feg  is 
Boltzmann's  constant,  and  T  is  the  ambient  temperature. 

In  Fig.  2  the  current  spectra  of  diode  Z.)  are  plotted. 
The  diode  was  also  biased  at  —  3  V  and  the  light  intensity 
was  controlled  so  that  the  currents  were  0.5,  1.0,  5.0,  and 
10.0  pA.  The  spectra  seem  to  have  the  same  shape  as  the 
spectra  of  diode  Z)|  with  the  main  difference  being  that  the 
slopes  of  the  low-frequency  components  are  equal  to 
—  0.83  ±0.01.  The  high-frequency  parts  of  the  spectra  are 
again  well  described  by  Eq.  ( 1 ).  In  Fig.  2  the  thermal  noise 
of  resistor  Rg  has  already  been  subtracted  from  the  spectra 
for  /  =  0.5  and  1.0  pA  but  not  for  the  other  two  current 
values.  This  allows  us  to  better  estimate  the  magnitude  of 
the  low-frequency  part  of  the  spectra.  Again,  the  magni¬ 
tude  of  the  low-frequency  part  of  the  spectra  scales  with 
the  current  squared. 

From  Fig.  3  it  is  clear  that  the  influence  of  the  bias 
voltage  in  the  range  from  —  0.5  to  -  3  V  on  the  spectra  is 
very  small  This  result  presented  here  for  diode  L^  is  rep¬ 
resentative  of  all  four  diodes. 

In  Table  I  we  summarized  the  significant  results  we 
obtained  from  measuring  the  noise  spectra  of  all  four  di¬ 
odes  under  reverse  bias:  D^,  Di,  L,,  and  L,.  For  each  diode 
we  measured  the  spectra  for  at  least  five  diflTerent  currents 
in  the  range  from  0.5  to  14  pA.  For  each  current  we  cal 
culated  the  normalized  magnitude,  5,(1  Hz)//^  and  the 
slope  y  and  then  averaged  these  values  over  the  currents. 
The  most  significant  observation  to  be  made  is  the  differ¬ 
ence  in  the  slopes  of  the  low-frequency  parts  of  the  spectra 
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FIG.  I  Current  noise  spectra  of  the  reverse-bia.sed  and  illuminated  diode 
D,  for  different  currents. 

The  exponent  of  the  spectra  of  both  diodes  D^  and  are 
very  close  to  —  1  whereas  the  exponents  for  diodes  and 
£,3  are  both  close  to  —  0.8.  Additionally,  the  normalized 
magnitude  of  the  spectra  is  smaller  for  diodes  Z),  and  by 
about  a  factor  of  4  than  it  is  for  diodes  and  Z.3.  This  is 
also  obvious  from  Fig.  4  where  we  present  the  spectral 
intensity  at  /  =  1  Hz  versus  the  current  for  the  four  de¬ 
vices.  The  solid  line  indicates  a  proportionality  with  the 
current  squared.  The  spectral  intensities  of  diodes  L|  and 
£,3  are  almost  identical  and  larger  than  the  noise  of  diodes 
£)|  and  ZJj.  The  noise  of  diode  is  about  an  order  of 
magnitude  smaller  than  that  of  diodes  £.,  and  £.3. 

B.  Not  Illuminated,  forward  biased 

In  this  second  set  of  experiments  the  diodes  were  kept 
in  the  dark  and  a  forward  voltage  was  supplied  so  that  the 
currents  would  be  equal  to  the  ones  during  the  first  exper¬ 
iments.  In  Figs.  5  and  6  the  current  spectra  of  diodes  £)| 
and  L^  are  presented,  respectively,  for  different  currents. 
The  general  shape  of  the  spectra  is  very  similar  to  what  we 
found  when  the  diodes  were  reverse  biased.  The  high- 
frequency  parts  of  the  spectra  follow  Eq.  ( I )  and  at 
smaller  frequencies  we  observe  a  frequency  dependence 
proportional  to  P  where  y  =  —  1.0  for  the  spectra  of  di¬ 
ode  Z?|  and  Y—  —0.8  for  the  spectra  of  diode  £.,.  The 
magnitude  of  the  low-frequency  noise  is  linearly  propor- 

TABLE  I  The  measured  average  noiK  intensity  at  I  Hr  [5, ( I  Hz)]  of  the 
diodes  under  rtvtnt  bias  and  illumination,  normalized  by  the  square  of 
ihe  current  The  slope  of  the  measured  low-frequency  spectra  is  y  The 
enfnes  in  the  table  are  averaged  over  at  least  live  different  currents  rang¬ 
ing  from  I  to  14  ^/A 


Diode 

No 

"  Hz  ') 

r 

D, 

0.7*01 

-0  99  *  0  03 

D, 

2  6*04 

-  0  99*006 

L, 

8  9*0  8 

-0  83*0.01 

L, 

9  8*08 

-  0  80  *  0  03 

10®  10'  10*  10*  10*  10* 


Frequency  (Hz) 

fig  2.  Current  noise  spectra  of  Ihe  reverse-biased  and  illuminated  diixle 
L,  tor  different  currents.  In  order  to  be  able  to  better  estimate  the  m.ig- 
niiude  of  the  low-frequency  part  of  the  spectra  it  was  necessary  to  sub¬ 
tract  the  contribution  of  the  thermal  noise  of  Ihe  bias  resistor  only  from 
Ihe  two  small-current  spectra. 

tionai  with  the  current.  This  is  in  contrast  to  what  we 
observed  for  the  reverse-biased  diode  spectra.  Figure  7 
shows  this  current  dependence  of  the  noise  spectra  for 
three  of  the  diodes.  By  the  time  we  took  the  forward-biased 
data  diode  Di  had  expired.  The  spread  in  the  magnitude  of 
the  noise  is  smaller  when  compared  with  the  previous  ex¬ 
periment.  Some  of  these  observations  are  summarized  in 
Table  II  where  we  give  the  average  values  of  the  current 
spectra  at  1  Hz  divided  by  the  current  as  well  as  the  aver¬ 
age  of  the  slope  of  the  low-frequency  parts  of  the  spectra. 

III.  DISCUSSION 

The  most  remarkable  result  is  that  the  exponents  ( y) 
of  the  low-frequency  parts  of  the  spectra  are  —  1.0  for 
diodes  £>|  and  both  under  forward  and  reverse  bias  and 
this  exp"':cnf  is  -  0.8  for  diodes  £,,  and  £,3,  again  both 
under  forward  and  reverse  bias.  Intriguing  as  well  is  the 
fact  that  the  magnitudes  of  the  noise  spectra  when  the 
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FIG  3  Current  noise  spectra  of  the  reverse-biased  and  illuminated  diode 
L,  for  different  currents  and  different  voltages 
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FIG.  4  Currenl  dependence  of  the  noise  intensity  at  I  Hz  for  the  four 
diodes  under  reverse  bias  and  illumination.  The  slope  of  the  line  is  ^  2 

diodes  are  forward  biased  are  almost  equal  for  diodes  D^, 
Z.i  and  Lj  (see  Table  II).  Although  there  is  some  difference 
between  the  two  sets  of  devices  when  the  diodes  are  reverse 
biased  the  spectral  densities  are  within  one  order  of  mag¬ 
nitude  (see  Table  I).  Both  these  observations  seem  to  in¬ 
dicate  that  the  mechanisms  that  are  responsible  for  the 
noise  are  present  whether  the  diodes  are  forward  or  reverse 
biased.  This  ts  surprising  when  one  considers  that  the  dc 
current  flow  is  governed  by  different  mechanisms  under  the 
two  bias  conditions.  This  statement  is  supported  by  the  fact 
that  the  current  dependence  of  the  magnitude  of  the  noise 
IS  different  for  the  two  bias  conditions  (see  Figs  4  and  7). 

Another  method  that  has  been  used  to  separate  out  the 
different  noise  producing  mechanisms  is  to  evaluate  the 
parameter  a ^  as  it  is  defined  by  Eq.  (2): 

S,(/)//^=a«//iV.  (2) 

In  this  equation'  S,(  f  )  is  the  spectral  intensity  of  the 
current  fluctuations,  /  the  dc  current,  /  the  frequency,  and 
iV  the  number  of  earners  in  the  system.  If  nothing  else. 
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Fraquancy  (Hz) 

FIG  5  Current  noise  specira  of  ihe  forward-biased  (no  illuminaiion) 
diode  D,  for  different  currents. 
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Fraquancy  (Hz) 

FIG  6.  Current  noise  specira  of  the  forward-biased  (no  illuminaiion) 
diode  Z.,  for  different  currents. 


an  is  a  dimensionless  measure  for  the  noisiness  of  the  de¬ 
vice  and  permits  as  such  a  comparison  between  devices.  By 
comparing  the  derived  values  of  a n  with  theoretical  pre¬ 
dictions  made  by  Handel,^  which  are  only  applicable  to 
true  1//  spectra  and  many  of  which  have  been  verified,^ 
the  noise  processes  can  be  identified.  However,  this  method 
is  hard  to  apply  in  the  case  of  reverse  biased  diodes  since 
we  do  not  have  an  mdepei.ueni  way  to  establish  the  num¬ 
ber  of  free  carriers  in  these  diodes. 

When  the  diodes  are  forward  biased  the  current  flow  is 
governed  by  a  transit  time  or  lifetime  ( r)  of  the  carriers.  In 
that  case  N  =  Ir/q  and  upon  substitution  into  Eq.  (2)  one 
Ands 

SAf)/l=^a,fi/Tf.  (3) 

Measurements  of  the  ac  impedance  of  the  diode  have  been 
used  before^  to  obtain  this  time  constant  from  resonances 
in  the  real  part  of  the  impedance.  Even  when  the  spectra 
have  an  exponent  unequal  to  —  1.0  can  we  deAne 
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FIG  7  Current  dependence  of  ihe  noise  iniensily  al  I  Hz  for  the  three 
diodes  under  forward  bias  ( no  illumination)  The  slope  of  the  line  IS  -i-  i 
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TABLE  II.  The  measured  average  noise  intensiey  at  I  Hz  [5,(  1  Hz))  of 
the  diodes  under  forward  bias  and  no  illumination,  normalized  by  the 
current  The  slope  of  the  measured  low-frequency  spectra  is  y  The  entries 
in  the  table  are  averaged  over  at  least  five  different  currents  ranging  from 
1  to  14  fjA 


Diode 

No 

5,(IHz)//(IO  "  A/Hz) 

r 

£>, 

5  1*11 

-  0.99  ±0  66 

5.4  ±0.7 

-O«2±0.03 

L. 

5.1  ±  1.5 

-  0  82  ±0  02 

„,r  5,(1  Hllr 

- ft - =  ,/  • 

where  a/f  is  not  dimensionless  anymore,  it  just  defines  the 
noisiness  at  /  =  1  Hz.  The  diodes  have  been  shown  to  be 
very-high-speed  devices;  they  have  been  operated  at  a  few 
GHz.  Therefore,  the  lifetime  in  these  devices  has  to  be  at 
least  10  ~  *  s,  resulting  in  values  for  that  are  larger  than 
3  X  10  “  Note  that  the  values  of  a///r  are  very  close  for 
;he  three  diodes  (see  Table  III).  The  value  for  Of,  (o/y) 
found  here  does  not  coincide  with  the  values  predicted  by 
Handel's  quantum  1//  noise  theory,  apart  from  the  fact 
that  the  frequency  exponent  of  the  low-frequency  parts  of 
the  spectra  are  not  equal  to  —  1.0.  Accordingly,  we  cannot 
attribute  the  noise  we  observe  to  these  mechanisms.  On  the 
other  hand,  the  fact  that  the  values  for  a^ia/f)  are  almost 
identical  for  the  three  diodes  might  suggest  that  apart  from 
a  surpnsing  uniformity  that  we  have  not  observed  as  yet, 
the  noise  is  of  a  fundamental  nature.  It  requires  further 
study  of  what  other  fundamental  noise  sources  can  exist 
besides  quantum  1//  noise. 

In  this  light  it  should  be  considered  that  earlier  exper¬ 
iments  on  Hg,  _  ,Cd,Te  p-i-n  diodes  revealed  values  for 
Off  in  the  range  of  2x  10“‘-2x  10“’;  the  frequency  ex¬ 
ponents  were  smaller  than  —  1.0  in  these  devices  as  well.’ 
The  noise  was  attributed  to  a  distribution  of  generation- 
recombination  centers  having  different  activation  energies. 
This  would  result  in  noise  spectra  that  are  summations  of 
Lorentzians  [l/(  I J,  This  so-called  generation- 
recombination  1//  noise  might  also  be  the  source  of  the 
noise  in  our  InGaAs  diodes. 

IV.  QUANTUM  EFFICIENCY 

In  the  search  for  a  fundamental  source  of  noise  one 
might  want  to  consider  fluctuations  in  the  quantum  effi¬ 
ciency.  One  should  realize  that  this  would  only  be  a  source 


TABLE  HI.  The  vitues  of  a,,  from  the  current  none  spectri  of  the 
forward  biased  diodes  Equation  (4)  is  used  to  calculate  a„  and  il  is 
assumed  that  r  =  1  x  10  ’  s 


Diode 

a.i 

No 

(m  unils  of  10  ') 

Y 

0, 

6  4 

10 

L, 

6  8 

08 

t, 

66 

08 

of  noise  when  light  is  shone  on  the  devices  as  we  did  when 
we  reverse  biased  the  devices.  This  would  not  explain  the 
noise  in  the  devices  under  forward  bias  Let  n  photons 
arrive  at  the  surface  per  second  and  let  the  reflection  coef¬ 
ficient  of  the  surface  be  r,  then  m  =  n(  1  —  r)  of  these  pho¬ 
tons  will  enter  the  semiconductor.  They  will  have  proba¬ 
bility  7j  to  generate  an  electron-hole  pair  resulting  in  the 
production  of  ,V  =  rjzif  1  —  r)  electron-hole  pairs  per  sec¬ 
ond.  It  is  shown  that  the  variance  of  A  equals 

varfA)  =  [var(/i)  -  (zi)  ]( 1  —  r)’)/  -i-  <A).  (5) 

The  angular  brackets  { (  > )  indicate  taking  the  ensemble 
average.  If  the  incandescent  light  has  a  Poissonian  distri¬ 
bution  [var(zi)  =  <«>)  then  the  distribution  of  the  electron- 
hole  pairs  will  be  Poissonian  as  well  and  the  high- 
frequency  shot  noise  will  be  Iql.  If  the  distribution  of  the 
photons  is  super-Poissonian  [varfzi)  >  (n)]  then 
var(A)>(A)  and  the  shot  noise  will  be  enhanced.  We 
have  not  observed  the  enhancement  of  the  shot  noise  in 
either  diode  and  therefore  conclude  that  the  light  source 
gives  a  Poissonian  distribution  of  photons. 

A  possible  source  for  low-frequency  fluctuations  might 
be  the  quantum  efficiency,  resulting  in  a  l//-like  spectrum. 
Since  the  photocurrent  is  proportional  to  the  quantum  ef¬ 
ficiency  one  can  show  that  the  relative  spectrum  of  the 
fluctuations  in  the  current  is  equal  to  the  relative  spectrum 
of  the  fluctuations  in  the  quantum  efficiency; 

S/I^  =  S,/ri\  (6) 

Since  r]  is  practically  independent  of  the  current,  the  same 
is  true  for  so  that  5,  varies  as  /^. 

This  model  would  explain  qualitatively  the  current  de¬ 
pendence  of  the  current  spectra  when  the  diodes  are  re¬ 
verse  biased  and  illuminated.  It  would  not  explain  the  mag¬ 
nitude  of  the  noise  we  observed. 

V.  SUMMARY 

The  following  observations  were  made  from  the  mea¬ 
sured  current  noise  spectra  of  two  types  of  InGaAs/InP 
photodiodes  under  reverse  bias  and  illumination  and  under 
forward  bias  without  illumination; 

( 1 )  For  each  set  of  devices  the  slopes  of  the  low- 
frequency  spectra  where  the  same  under  both  bias  condi¬ 
tions. 

(2)  Under  forward  bias  both  sets  of  devices  showed  a 
current  dependence  of  the  noise  that  is  linear. 

(3)  Under  reverse  bias  both  sets  of  devices  showed  a 
current  dependence  of  the  noise  that  is  quadratic. 

(4)  The  values  of  a,f  obtained  from  the  noise  spectra 
under  forward  bias  are  not  compatible  with  Handel's  quan¬ 
tum  !//  noise;  they  seem  to  be  indicative  of  1//  noise  due 
to  a  distribution  of  recombination  centers. 

From  the  first  observation  one  might  infer  that  the 
mechanism  responsible  for  the  noise  is  the  dominant  one 
under  both  bias  conditions  This  is  not  obvious  since  the 
mechanism  controlling  the  current  flow  under  the  two  con¬ 
ditions  IS  different.  The  second  and  third  observation  sup¬ 
port  the  fact  that  the  current  flow  in  the  diodes  is  con¬ 
trolled  by  different  mechanisms  under  the  two  conditions. 
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We  have  shown  that  fluctuations  in  the  quantum  eth- 
ciency  could  give  rise  to  the  quadratic  current  dependence 
of  the  noise  intensity.  Of  course,  this  mechanism  cannot 
explain  the  noise  when  the  diodes  are  forward  biased  and 
not  illuminated. 

A  distribution  of  traps  with  different  activation  ener¬ 
gies  seems  to  be  a  likely  candidate  to  explain  our  observa¬ 
tions.  It  would  explain  the  identical  slope  of  the  spectra 
under  the  two  bias  conditions.  The  different  current  depen¬ 
dence  of  the  noise  under  these  two  conditions  would  be 
explained  by  the  different  mechanisms  controlling  the  cur¬ 
rent  flow  The  latter  ones  determine  how  a  local  noise 
source  contributes  to  the  noise  as  it  appears  at  the  termi¬ 
nals  We  would  like  to  emphasize  that  for  the  three  diodes 
under  forward  bias  the  value  oi  S,(  IHz)//  is  remarkably 


close  (Table  11).  This  might  be  indicative  of  a  more  luii- 
damental  noise  source. 
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NOTES 

EXTENSION  OF  THE  HOOGE  EQUATION  AND  OF  THE  HOOGE 
PARAMETER  CONCEPT 

(Received  30  June  1988:  in  revised  form  9  June  1990) 


INTRODUCTION 

The  Hooge  equation  is  extended  to  the  case  where  the  noise 
spectrum  is  of  the  Torm  I//’  with  y  slightly  different  from 
unity.  This  leads  to  a  generalization  of  the  Hooge  parameter 
that  corrects  an  earlier  ambiguity. 


DISCUSSION 

In  1969  Hooge  published  his  well-known  equation(l) 

SlU) 

I‘  fN  '  ' 

giving  an  empirical  description  of  the  releative  I  //  noise  of 
a  uniform  semiconductor  sample  carrying  a  current  /.  Here 
5,(/)  is  the  spectral  intensity  of  the  current  fluctuations,  A( 
the  number  of  carriers  in  the  sample  and  Sh  empirical 
constant,  known  as  the  Hooge  parameter,  that  characterizes 
the  noise. 

Later,  eqn  (I)  was  extended  to  include  mobility  and 
diffusion  I// noise  and  applied  to  devices.  Equation  (I)  may 
then  be  writu*n(2,3|  as 

^,(/)  *M  .  •So(/)  «H  . 

^  JS-  ~D^'fN 

The  two  equations  are  connected  by  means  of  the  Einstein 
relation  qD  •kT,. 

In  the  case  of  spatially  distributed  I// noise  sources  eqn 
(I)  can  be  written  as 

»N  . 

n  fNfx)&x'  ^ 

Here  5,(.r,/)  is  the  current  spectrum  in  a  section  Ax  at  x; 
Nfx)  is  the  number  of  carrien  per  unit  length  at  x. 

Equation  (I)  can  be  used  for  evaluating  from  the 
measured  spectrum  5,(/)  If  the  spectrum  is  exactly  \lf,  we 
may  write  S,{f)  •  Sif  hlf:  solving  for  yields 

Nf5,(f)  NS,(\) 

“m  ”  — Jl —  ■  —Jj— 

The  value  of  thus  obtained,  is  independent  of  the 
frequency  /at  which  5,(/)  is  measured. 

At  first  the  first  half  of  eqn  (2)  was  also  used  without 
having  ascertained  that  the  spectrum  was  exactly  Ilf.  That 
is  ixh  was  expressed  as 

irrespective  of  the  frequency  dependence  of  5,(/).  But 
this  led  to  ambiguous  results,  in  that  a„  now  depended  on 
the  frequeiKy  /.  For  if  5,(/)  had  a  I,/'  dependence: 
5,(/)  •  S,{\)lf\  one  finds  by  substitution 


which  depends  on  the  frequency  /  and  is  only  equal  to  eqn 
(2)  at  I  Hz.  Usually  the  error  is  relatively  small,  however. 
Note  that  5,(1)  has  the  dimension  -A’(Hz)y*'. 


To  remedy  the  situation  we  generalize  eqn  (I).  If  the 
spectrum  S,{f)  =  5,(1)//  with  y  #  1.  the  simplest  general¬ 
ized  form  is 


Substituting  for  5,(/)  yields 

/V/’5,(/)  /VS,(1) 

«H - p —  -  -jr-  (#») 

which  is  identical  with  eqn  (2)  The  expression  for  a^  is  now 
independent  of  the  measuring  frequency  and  does  not 
contain  the  exponent  y  explicitly.  Note  that  3^  has  the 
dimension  (Hz)y 

If  we  compare  eqns  (3a)  and  (4a)  we  see  that  the 
correction  factor  is  /'■’.  where  /  is  the  measurement 
frequency  For /  =  I  Hz  the  correction  factor  is  unity  for  all 
values  of  y;  for /  =  10  Hz  and  y  =  0.90  the  correction  factor 
is  1.26.  whereas  for  /=>  100  Hz  and  y  «0.90  it  is  1.58.  For 
reasonable  accuracy  one  should  measure  close  to  I  Hz(4). 
Finally  eqn  (4)  may  be  written 

/'5,(/)-^/'-5,(l)  (5) 

This  is  easily  adapted  to  devices  with  distributed  noise 
sources  such  as  MOSFETs  and  short  p-n  junction  diodes. 

As  a  first  example  consider  a  MOSFET  operating  in  the 
linear  mode  (K.  4  F.,,)|5|  If  the  device  has  a  length  L  a 
drain  voltage  F  and  a  current  /.  and  n  is  the  carrier  mobility, 
then 

■S((/)~  .  or  /'S,(/)  -  5,(1)  - — —  (6) 

As  a  second  example  consider  a  short  n  --p  diode.  The 
spectrum  is  now(3.5) 


,  3„e/  r  A/(0)  1  /V(0)  5„f4'. 

S,lf)^— —  n  — »  I  + -2— ' ,  (7) 

2/T*.  N(W^)  D. 

Here  >♦',  is  the  length  of  the  p-region.  5„  the  electron 
recombination  velocity  at  thi  p<ontact  (6  x  10*  cm  in  n- 
type  Si).  D,  the  diffusion  constants  for  electrons  and 
-  H'‘l2D,  the  diffusion  time  of  the  electrons  through  the 
p-region.  Hence 

/^5,(/)-5,(l)»^ln^l (7a) 

The  extension  does  not  apply  too  well  to  GaAs  MESFETs 
since  the  spectrum  has  g-r  "bumps"  superimposed  upon  a 
l//background|6).  If  one  does  not  require  great  accuracy, 
however,  eqn  (3a)  or  eqn  (4a)  can  be  used  as  approximations 
to  find  "order  of  magnitude"  values  of 

If  one  wants  to  compare  the  quantum-theoretical(7.g)  and 
the  experimental  values  of  the  Hooge  parameter  Xh,  ones 
uses  the  experimental  values  given  by  eqns  (4a),  (6)  or  (7a), 
whereas  the  theoretical  values  are  taken  from  Handel's 
equation(7.g]. 

43,  At’’ 

(8) 
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Notes 


where  1,=  1/(137)  is  the  fine  structure  constant,  c  the 
velocity  of  light  and  dv  the  vectonal  change  in  carrier 
veloaty  in  the  collision  process. 

It  should  furthermore  be  borne  in  mind  that  the  quantum 
theory  of  I//  noise  results  in  a  I//'  spectrum  with  y  very 
close  to  unity.  If  deviations  occur,  either  in  y  or  in  one 
should  investigate  whether  a  superposition  of  spectra  might 
be  present.  In  some  FETs  a  generation-recombination  type 
I//'  spectrum  sometimes  masks  the  wanted  1,/ spectrum; 
this  unwanted  noise  source  must  then  be  corrected  for(9|.  In 
vacuum  photodiodes  a  quantum  I //spectrum  and  a  classical 
surface  I,'/”  spectrum  are  present  simultaneously  and  must 
be  separated  before  an  analysis  can  be  attempted[IO|.  Care 
must  be  taken  that  conclusions  are  not  made  prematurely. 
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NOISE  IN  THE  QUANTUM  EFFICIENCY  i>  OF  p*-n  DIODES  DUE 
TO  FLUCTUATION  IN  THE  SURFACE  GENERATION-RECOMBINATION 

OF  CARRIERS 

(Rtceh-ed  7  April  1990;  m  revised  form  3!  May  1990) 


According  to  Meingailis  and  Harman!  I  ]  the  quantum 
efficiency  ti  in  z  p'’~n  photodiode  with  surface  generation 
and  recombination  of  carriers  is  given  by 


coshfd/f.,)  -t-  («,//..  )sinh(i//f.,) 

Here  r  is  the  photon  reflection  coefficient  at  the  surface. 
T,  the  carrier  lifetime,  d  the  depth  of  the  junction,  is  the 
hole  diffusion  length,  and  s  the  surface  recombination 
velocity. 

We  now  calculate  the  noise  specturm  of  7.  Since  s 
fluctuates  in  a  I//  fashion,  7  will  fluctuate  in  the  same 
fashion.  Hence 

i7 _ (l-0((t,/t,)s.nh(d/X,)i 

[cosh(d/t,)  +  (.fT,/Z,,)sinh(d/f.,)H 
Putting  Q{s,  t,.d)  «(«,/£, )tanh(d/iL,)  yields 

7  I  -t-  0(r.r,,d)  j 

Making  a  Fourier  analysis  and  calculating  the  power  spec¬ 
trum  yields 

For  long  diodes  dlL,>  I  and  hence  tanh(d/t,)»  I.  so 
that 


If  diL,  d  I.  we  have  tanhfd/Z.,)  ^dlL^  and  0^  d  I,  so  that 

=  ,4b) 

7-  5- 

In  either  case  5,(/)  =  const  s//  so  that  5,(/)  has  a 
spectrum  that  is  nearly  independent  of  cuiTent(2j. 

In  modem  p‘~i~n  diodes  recombination-generation 
usually  occurs  in  the  bulk|3).  In  that  case  the  above  theory 
does  not  apply 
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Abstract : 

Reported  here  are  measurements  of  the  spectral  intensity  of  the  current 
fluctuations  in  double-barrier  cjuantum  well  resonant  tunneling  diodes  as  a 
function  of  temperature  and  bias  current.  Two  types  of  devices  were  studied:  one 
with  AlAs  barriers  and  GaAs  well  and  contact  regions,  and  the  other  has 
Alo  jGa,  ,As  barriers.  The  frequency  range  covered  is  1  Hz  to  100  kHz  and  the 
temperature  range  is  78  K  to  400  K.  The  noise  spectra  are  decon.posed  in  a  1/f 
part,  resulting  in  the  magnitude  of  the  1/f  noise,  and  contributions  due  to 
carrier  trapping,  resulting  in  the  activation  energies  of  the  traps.  It  is  found 
that  a  reduction  of  the  A1  contents  of  the  barrier  material  reduces  the  number 
of  traps  and  further  the  magnitude  of  the  1/f  noise  is  practically  independent 
of  the  temperature  and  of  the  A1  content  of  the  barrier. 


Double-barrier  resonant  tunnelling  structures  (DBRTS)  have  attracted  much 
attention  recently  because  of  their  functionality  in  a  wide  variety  of 
applications  that  include  frequency  multipliers,  parity  generators,  multi-state 
memory,  analog-digital  converters',  opto-electronic  devices^,  etc.  While  there 
would  be  variations  of  present  or  potential  applications  of  resonant  tunnelling 
devices,  DBRTS  offer  not  only  the  key  attributes  for  the  applications  mentioned 
above,  but  also  tender  opportunities  for  experimental  studies  on  quantum  effects 
in  carrier  transport^’  .  In  the  last  ten  years  the  impressive  advances  in  the 
molecular  beam  eprtaxy  hive  contributed  greatly  to  the  amelioration  of  the  DBRTS 
making  possible  more  refined  studies  of  the  characteristics  of  the  devices  such 
as  the  noise  performance  as  it  correlates  with  the  presence  of  defects, 
impurities,  and  tunnelling  processes.  We  report  here  on  an  extensive  study  of 
just  this  low-frequency  noise  behavior  of  two  types  of  DBRTS:  one  with  A1 
barriers  and  GaAs  electrodes  and  well,  and  another  type  with  AlGaAs  barriers  and 
GaAs  electrodes  and  well.  The  measured  current  noise  spectra  are  interpreted  in 
terms  of  1/f  noise  and  generation-recombination  (g-r)  noise.  The  noise 
spectroscopy  has  been  applied  before,  to  bulk  semiconductors**  as  well  as  to 
DRTBS’,  and  has  proven  to  be  very  successful  in  identifying  carrier  traps . 

The  noise  behavior  of  two  types  of  devices  was  investigated.  The  symmetric 
tunnel  structures  were  grown  by  MBE  on  an  n'  substrate  and  are  identical 
except  for  the  barriers  which  are  3nm  thick  AlAs  for  one  device  type  and  5  nm 
thick  Alj  jGa, ,A3  for  the  other  type.  The  details  of  the  layer  structure  are 
presented  in  Fig.  1.  In  all  the  doped  layers  Si  was  used  as  a  dopant.  Current- 
voltage  characteristics  of  the  two  devices  both  at  77  K  and  300  K  are  shown  in 
figures  2a  and  2b.  The  peak-to-valley  ratio  of  the  current  of  the  AlAs  barrier 
device  is  about  13  at  77  K  and  decreases  to  about  3  at  room  temperature.  For 
the  AlGaAs  barrier  device  the  current  peak-to-valley  ratio  is  about  5  at  77  K, 
dropping  to  2.2  at  room  temperature,  and  further  decreasing  to  13  at  400  K. 

In  order  to  measure  the  noise  as  a  function  of  temperature  the  devices  were 
mounted  in  a  flow  cryostat.  The  noise  signal  was  fed  into  a  low-noise 
amplifier  and  detected  by  a  Fast  Fourier  Transform  dynamic  signal  analyzer.  In 
calculating  the  spectra  allowance  was  made  for  the  amplifier  noise  as  well  as 
the  non-flatness  of  the  amplifier  gain.  The  resultant  noise  spectra  usually 
consisted  of  a  1/f  component,  Lorentzian  shaped  bumps  due  to  trapping  and  de- 
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trapping  of  carriers,  and  a  white  part  due  to  the  thermal  noise. 

The  experimental  investigation  is  divided  into  three  parts:  (i)  Spectral 
intensity  of  the  current  fluctuations  as  a  function  of  temperature  at  a  fixed 
current,  (ii)  Current  dependence  of  the  noise  for  voltages  smaller  than  the 
peaic-voltage  (Vp)  ,  and  (iii)  Current  dependence  of  the  noise  for  voltages 
larger  than  the  valley-voltage  (V,)  . 

To  study  the  temperature  dependence  of  the  noise  the  diodes  were  biased  at 
a  constant  current  of  0.7  mA,  relatively  close  to  the  peak,  and  the 
temperature  was  varied  between  78  K  and  400  K.  The  spectra  consist  of  a  1/f 
part  and  Lorentzian-shaped  bumps  due  to  generation  and  recombination  (g-r 
noise)  of  the  carriers.  Upon  multiplication  of  the  spectra  by  the  freqpjency 
the  1/f  contribution  will  appear  independent  of  frequency  and  the  bumps 
associated  with  the  g-r  noise  will  appear  as  peaks.  From  these  graphs  one  can 
readily  extract  the  magnitude  of  the  1/f  noise  and  the  frequencies  (fp)  for 
which  the  g-r  peaks  occur.  A  subsequent  plot  of  the  square  of  the  temperature 
T  divided  by  fp  as  a  function  of  the  reciprocal  temperature  permits  the 
calculation  of  the  activation  energies  of  the  traps  from  the  slopes  of  the 
curves . 

In  Fig.  3  we  show  these  Arrhenius  plots  for  the  device  with  the  AlAs 
barriers  and  for  the  device  with  the  Al^  jGa,,  ,As  barriers.  For  the  mentioned 
range  of  temperatures  five  trapping  levels  could  be  detected  in  the  AlAs 
barrier  device  with  activation  energies  of  0.017,  0.15,  0.21,  0.38,  and  0.55 
eV.  Under  the  same  bias  current  and  in  the  same  temperature  range  the  device 
with  the  Alg  jGao  ,A3  barriers  showed  only  two  traps  with  activation  energies  of 
0.16  and  0.55  eV:  these  energies  are  close  to  two  of  the  energies  we  found  in 
the  first  device.  Since  the  aluminum  content  of  the  barrier  material  is  the 
main  difference  in  the  device  structures  we  attribute  the  presence  of  the 
three  other  traps  to  the  larger  A1  mole  fraction.  Apparently,  there  is  some 
trap  assistance  in  the  tunneling  process  and  sweeping  the  temperature  allows 
us  to  probe  these  traps  by  making  noise  measurements. 

In  Fig.  4  we  present  the  magnitude  of  the  1/f  component  we  observed  in  the 
noise  spectra  as  a  function  of  the  temperature.  The  bias  current  (1^)  was  held 
constant  at  0.7  mA  for  both  devices.  We  want  to  mention  here  that  the 
frequency  exponent  of  the  spectra  is  equal  to  -1.00  +  0.05.  Due  to  the 
magnitude  of  the  g-r  noise  in  the  AlGaAs  barrier  device  (solid  triangles)  at 
temperatures  below  190  K  it  was  net  possible  reliably  to  extract  the  magnitude 
of  the  1/f  noise. 

Our  measurements  indicate  that  the  magnitude  of  the  1/f  noise  component  is 
relatively  constant  over  the  entire  temperature  range  for  both  devices  when 
the  diodes  are  biased  close  to  the  peak  (0.7  mA) .  This  strongly  suggests  that 
the  1/f  component  is  associated  with  the  tunneling  current  since  this  current 
is  independent  of  the  temperature  in  first  order.  Because  the  device  is  biased 
relatively  close  to  the  peak  we  expect  the  main  current  component  to  be  due  to 
tunneling.  We  also  observe  that  the  magnitude  of  the  1/t  noise  is  practically 
the  same  in  the  two  different  devices.  This  seems  to  imply  that  the  physics 
(tunnelling  through  the  barrier  at  this  operating  point)  that  governs  the 
transport,  and  therefore  the  noise,  is  the  same  for  these  devices,  although 
the  materials  properties  of  the  barriers  differ. 

The  last  part  of  our  study  of  the  noise  behavior  of  resonant  tunnel  diodes 
involved  the  current  dependence  of  the  1/f  noise  component  both  at  different 
temperatures  as  well  as  for  bias  points  with  0  <  V  <  Vp  and  for  bias  points 
with  V  >  V,  for  the  two  devices.  The  peak  current  occurs  for  V=Vp,  the  valley 
current  occurs  for  V,.  In  Fig.  5  we  show  the  spectral  intensity  of  the  current 
fluctuations  at  1  Hz  versus  the  bias  current  measured  at  room  temperature  for 
the  AlAs  barrier  device.  The  solid  circles  indicate  biaspoints  with  V  <  Vp  and 
the  solid  triangles  represents  the  noise  for  the  biaspoints  for  which  V  >  V„. 
In  the  domain  where  the  current  is  due  mostly  to  tunneling  through  the 
barriers  the  noise  increases  with  tlie  current  and  can  be  described  by  I‘ ‘. 

When  the  voltage  is  larger  than  V,  the  current  dependence  of  the  noise  is  even 
stronger,  close  to  l'.  The  latter  current  dependence  is  expected  for  a 
resistor-like  device.  In  this  bias  regime  the  current  is  due  mostly  to 
carriers  that  are  excited  over  the  barriers.  It  is  important  to  observe  that 
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for  the  same  current  the  device  is  noisier  when  biased  at  voltages  smaller 
than  the  peak  value. 

When  the  temperature  is  lowered  to  78  K  the  current  dependence  of  the  noise 
is  less  strong,  the  device  noise  tends  to  behave  more  like  that  of  a  regular 
diode,  although  the  magnitude  of  the  noise  does  not  change  very  much  (Fig.  5, 
open  circles) . 

For  the  AlGaAs  barrier  device  these  measurements  were  repeated  at 
temperatures  of  400,  300,  and  200  K  (see  Fig.  6) .  At  the  highest  temperature 
we  observe  roughly  the  same  behavior  as  we  observed  for  the  AlAs  barrier 
device;  the  current  dependence  of  the  noise  for  the  smaller  voltages  (circles) 
is  less  strong  than  it  is  for  the  voltages  larger  than  Vv  (triangles) . 

However,  the  data  for  which  V  <  Vp  show  a  knee  close  to  the  maximum  current 
point.  This  is  more  dramatic  for  the  300  K  data.  The  200  K  data  do  not  show 
this  behavior  very  strongly.  A  point  to  note  is  that  as  the  temperature  is 
lowered  from  400  K  to  300  K  to  200  K  the  noise  in  the  after-valley  regime 
increases  with  respect  to  the  noise  when  the  device  is  biased  in  the  before¬ 
peak  regime.  This  is  consistent  with  the  remarks  we  made  earlier  that  the  1/f 
noise  in  the  before-p >ak-regime  is  associated  with  the  tunneling  process 
whereas  the  1/f  noise  in  the  after-valley-regime  is  due  to  other  processes  and 
has  a  stronger  temperature  dependence. 

In  summary  measurements  are  presented  here  of  the  low-frequency  noise  of 
ti*pes  of  resonant  tunnelling  diodes,  one  with  AlAs  barriers  and  one  type 
with  A1o.3  Gao. 7  As  barriers.  The  temperature  dependence  of  the  1/f  noise  is 
investigated  when  the  devices  are  biased  just  below  the  peak  in  the  current - 
voltage  characteristic  and  found  to  be  very  small.  Additionally,  the  magnitude 
of  the  1/f  noise  is  very  similar  in  both  device  types.  Both  these  observations 
are  believed  to  be  consistent  with  the  main  current  coitponent  in  this  bias 
regime:  the  tunneling  through  the  barriers. 

Second,  from  generation-recombination  noise  measurements  trap  activation 
energies  were  extracted.  We  noticed  five  traps  in  the  AlAs  barrier  device  and 
two  in  the  AlGaAs  barrier  device  with  activation  energies  equal  to  two  of  the 
trap  energies  we  found  in  the  AlAs  barrier  device.  This  seems  to  support  the 
idea  that  the  larger  A1  concentration  in  the  barrier  material  introduces 
carrier  traps  in  the  device. 

Third,  we  compared  the  current  dependence  of  the  1/f  noise  in  the  devices  in 
the  two  biasing  regimes:  before  the  peak  and  after  the  valley.  We  observed  a 
stronger  current  dependence  of  the  magnitude  of  the  1/f  noise  when  the  device 
was  biased  in  the  tunnelling  regime  than  when  it  was  biased  in  the  field- 
assisted  thermionic  emission  mode.  The  magnitude  of  the  noise  in  the  latter 
mode  of  operation  is  smaller  consistent  with  expectations. 
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Figure  Captions 

1.  Device  layer  configuration. 

2.  Current-voltage  characteristics  of  (a)  the  AlAs  barrier  and  (b)  the  AlGaAs 
barrier  device  at  different  temperatures. 

3.  Activation  energies  as  deduced  from  the  generation-recombination  spectra 
for  the  AlAs  barrier  (circles)  and  the  AlGaAs  barrier  (triangles)  device.  The 
quiescent  current  during  the  noise  measurements  was  0.7  mA. 

4.  Spectral  intensity  of  the  current  fluctuations  (1/f  component)  at  f  -  1  Hz 
and  an  operating  current  of  0.7  mA  as  a  function  of  the  temperature  for  the 
AlAs  barrier  (circles)  and  the  AlGaAs  barrier  (triangles)  device. 

5.  Spectral  intensity  of  the  1/f  current  fluctuations  at  f  =  1  Hz  versus  the 
operating  current  for  the  AlAs  barrier  device.  The  solid  circles  represent  the 
data  taken  at  room  temperature  and  for  voltages  smaller  than  the  I-V  peak 
voltage.  The  solid  triangles  are  room  temperature  data  measured  as  the  voltage 
was  larger  than  the  valley  voltage.  The  open  circles  specify  the  data  for  T  = 
78  K  and  voltages  smaller  than  the  peak  voltage.  The  letter  B  labels  before 
the  peak,  A  labels  after  the  valley. 

6.  Spectral  intensity  of  the  current  fluctuations  at  f  =  1  Hz  of  the  AlGaAs 
barrier  device  versus  the  operating  current  at  T  -  400  K,  T  =  300  K,  and  T  - 
200  K.  The  circles  represent  data  taken  with  voltages  smaller  than  the  peak 
voltage  and  the  triangle  show  the  data  for  applied  voltages  larger  than  the 
valley  voltage. 
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ABSTRACT 


It  is  shown  that  the  present  model  used  to  describe  the  silicon 
bipolar  Junction  transistors  is  far  from  being  precise.  In  general  it 
demonstrates  the  approximate  locations  of  the  constituents  of  both 
the  base  and  the  collector  noise  current  sources  and  the  relative 
magnitudes,  but  shows  minimally  about  the  shape  and  type  of  the  noise 
sources.  For  the  transistors  under  investigation,  we  found  that  the 
collector  noise  spectra  are  generally  1//^  type  with  y  *  2.0.  This  could 
mean  that  the  dominant  noise  source  in  the  collector  current  is  G-R 
type  noise.  Given  this  fact,  the  Hooge  parameters  obtained  from  these 
collector  noise  spectra  have  either  no  meaning  or  can  only  be 
interpreted  as  an  upper  limit  to  the  an-  The  base  noise,  on  the  other 

hand  conforms  to  the  Klelnpennlng-van  der  Ziel  l/f  diffusion  theory. 
One  can  therefore  determine  whether  the  noise  source  is  of  a 
fundamental  nature  by  comparing  and  contrasting  the  experimentally 
obtained  Hooge  parameters  oh’s  with  theoretically  calculated  values. 

In  Gallium-Arsenide  (GaAs)  heterojunction  transistors,  same 

result  was  obtained  for  the  collector  in  some  preliminary  experiments. 

We  have  again  G-R  type  noise  appearing  in  the  collector  current 

spectra.  As  for  the  base  current  spectra,  it  is  different  than  that  of  the 

Si  transistors,  that  is,  the  near  l/f  base  current  spectra  have  a  base 
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current  dependence  as  Ib:  contrary  to  its  counterpart  in  Si  transistors 
which  obey  the  1  //  diffusion  theory. 


V 


Since  ttH  alone  is  not  always  sufficient  to  determine  the  noise 
source.  New  quantities  fcdcJ  ^ind  feds)  as  1/f  cut-off  frequencies  for 
both  the  collector  and  the  base  are  introduced.  These  generalized 
quantities  can  be  used  to  determine  whether  the  devices  have  l/f 
diffusion  noise  or  not.  Usually  if  the  cut-off  frequency  varies  slowly 
with  the  respective  bias  current  it  is  l/f  diffusion  noise,  otherwise  it 
is  not  likely  to  be. 

Finally.  Oh's  found  in  the  experiments  are  compared  to  the 

theoretical  values  calculated  by  Handel.  Van  Vliet  and  others  to  show 
that  we  have  devices  exhibiting  collision  l/f  noise.  Umklapp  noise 

and  perhaps  intervalley  scattering  noise.  In  most  cases,  we  have 
combinational  effects  therefore  the  values  of  an’s  tend  to  be  larger 

than  the  individual  quantum  limits. 
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The  main  body  of  the  thesis  is  virtually  divided  into  two  parts, 
Chapters  I, II, ill  and  Chapters  IV, V.  In  the  first  three  chapters  an 
introduction  to  noise  phenomena  is  attempted  and  the  basic  concepts 
related  to  noise  are  presented.  The  flicker  (1/f)  noise  is  first 
discusc  in  Chapter  II  along  with  the  Hooge  parameter  an  which  is 
introduced  as  an  empirical  parameter  which  measures  the  noisiness 
of  several  solid-state  devices  exhibiting  1/f  noise.  The  incoherent 
state  1/f  noise  theory  and  the  coherent  state  1/f  noise  theory  are 
mentioned  in  Chapter  III.  The  equations  yielding  the  Hooge 
parameter,  oh.  are  given  in  each  case  along  with  the  theoretical 
predictions  regarding  the  values  of  oh  in  a  number  of  different 
devices. 

Whereas  the  first  part  of  the  thesis  is  devoted  to  the  theory  of 
1/f  noise,  the  second  part  (Chapters  IV  and  V)  is  devoted  to 
experimental  results  taken  from  p-channel  JFETs  and  both  short  and 
long  n-channel  JFETs.  The  measurement  procedure  and  the  model  we 
assume  in  order  to  describe  the  noise  behavior  of  JFETs  are 
presented  in  Chapter  IV.  For  a  saturated  channel  JFET,  we  split  the 
channel  into  two  parts,  the  ohmic  part  in  which  the  mobility  is 
regarded  to  be  constant  and  the  saturated  part  in  which  the  velocity 
equals  the  saturation  velocity.  The  noise  is  accounted  for  by 
assuming  two  current  noise  sources,  one  for  each  part  of  the 
channel,  characterized  by  the  same  oh-  The  equation  which  relates 
the  value  of  oh  to  experimentaly  accessible  variables  is  given  in 
this  chapter. 


Experimental  data  taken  from  JFET  devices  of  three  different 
types  are  presented  in  Chapter  V.  The  noise  behavior  of  2N4867  n- 
channel  JFETs  with  channel  length  of  30|im  is  investigated  first. 
Noise  measurments  at  room  temperature  indicate  that  the  noise 
originates  from  Umkiapp  scattering  processes.  JFETs  with  a  p- 
channel  of  the  .  3  length  (U168)  are  also  measured  at  different 
temperatures  and  below  saturation.  The  product  of  aHxp  (mobility) 
seems  to  be  insensitive  to  temperature  changes  for  a  range  of, 
about,  220  to  350  K  for  these  devices.  Finally,  noise  measurements 
are  carried  out  on  short  channel  (2.6pm)  n-JFETs  (2SK152).  The 
dependence  of  an  on  temperature  and  drain-to-source  voltage,  V^s  is 
investigated.  The  obtain  values  (>5x10'iO)  are  one  order  of 
magnitude  lower  than  the  theoretical  predictions  and  indicate  a 
possible  dependence  of  oh  on  the  electic  field  within  the  channel. 

The  references  used  in  the  course  of  this  research  work  are 
given  together  at  the  end  of  the  thesis. 
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II.  THEORETICAL  RESULTS 


1.  lairQ.du.giion 

Progress  has  been  achieved  during  1989  and  1990  in  the  study  of  nonlinear 
systems  which  generate  chaotic  1/f  fluctuations,  in  the  application  of  the  Quantum 
1/f  Theory^'^’^'l  ^  to  various  materials  used  in  small  and  ultrasmall  electronic 
devices,  and  in  the  application  of  the  Quantum  1/f  Theory  to  electronic  devices. 

The  close  similarity  of  the  classical  and  quantum  1/f  theories,  and  the  initial 
development  of  the  quantum  1/f  theory  by  the  author  out  of  his  efforts  to  quantize 
his  classical  turbulence  theory,  have  led  to  sustained  efforts  of  the  author  aimed  to 
integrate  all  his  theories  as  various  forms  or  realizations  of  a  fundamental  notion  of 
chaos  in  nonlinear  systems.  During  this  grant  period,  these  efforts  finally  beared 
fruits.  A  general  sufficient  criterion  was  formulated,  allowing  to  identify  the 
nonlinear  systems  which  exhibit  1/f  spectra.  This  criterion  is  presented  in  Sec.  3 
below.  It  is  followed  in  Sec.  4  and  in  the  Appendix  by  examples,  in  which  the 
criterion  is  applied  to  the  classical  and  quantum  mechanical  forms  of  the  author's  1/f 
noise  theory.  These  examples  clarify  the  physical  meaning  of  the  new  criterion. 

For  the  practical  application  of  the  Quantum  1/f  Theory  it  is  necessary  to 
derive  the  quantum  1/f  fluctuations  of  various  kinetic  (transport)  coefficients  which 
characterize  the  materials  used  in  electronic  and  microelectronic  applications,  from 
the  author's  fundamental  quantum  1/f  formula.  The  latter  is  applicable  only  to  cross 
sections  and  rates  of  elementary  processes.  Most  important  is  the  calculation  of 
mobility  fluctuations  in  Si,  GaAs  and  Hgi.xCd^Te.  An  earlier  calculation  (Kousik,  Van 
Vliet,  Handel,  1985)  of  mobility  fluctuations  in  Si  and  GaAs  is  replaced  in  Sec.  5  by  a 
more  rigorous  calculation,  based  on  the  new  quantum  1/f  cross-correlations, 
presented  at  the  1989  Conference  on  Noise  in  Physical  Systems.  The  new  calculation 
yields  increased  1/f  noise,  and  is  in  vey  good  agreement  with  the  experiment. 
Although  not  presented  here  yet,  we  also  performed  a  Monte  Carlo  simulation  for 
HgCdTe  during  this  period.  The  simulation  has  yet  to  be  improved  and  compared  with 
the  experiment. 

Primarily,  we  have  tried  to  improve  the  application  of  quantum  1/f  theory  to 
the  collector  noise  of  bipolar  transistors.  This  short  calculation  is  presented  in  Sec.  2. 

In  the  same  time,  many  new  contributions  to  the  quantum  1/f  theory  and 
experiment  were  published  by  other  workers  in  the  field,  considerably  advancing 
the  field  of  infra-quantum  physics  and  quantum  1/f  noise  in  high-technology 
applications.  These  new  contributions,  as  well  as  new  PhD  thesis  work  in  this  field 
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and  contributions  presented  at  the  IV  Conference  on  Quantum  1/f  Noise  and  other 
Low-Frequency  Fluctuations  and  included  in  the  Conference  Proceedings^'^,  are 
included  in  the  updated  General  Quantum  l/f  Bibliography  appended  to  this  Report. 


2.  Collector  quantum  l/f  noise  in  BJTs 

l/f  noise  in  bipolar  junction  transistors  (BJT's)  was  elegantly  treated  by  van 
der  Ziel^'^  who  applied  a  Hooge-type  approach  similar  to  Kleinpenning's  treatment'^ 
of  pn  junctions,  and  used  experimental  data  to  determine  the  Hooge  constant  which 
was  in  turn  compared  with  the  quantum  l/f  theory.  However,  since  the  BIT  is  a 
minority  carrier  device,  it  requires  the  application  of  the  quantum  l/f  equation^'^ 
from  the  beginning,  for  the  correct  interpretation  of  the  number  of  carriers  in  the 
denominator  of  the  Langevin  noise  source. 

In  the  most  elementary  model*  of  a  BJT,  the  collector  current  Ic  arises  from 
minority  carriers  injected  from  the  emitter  into  the  base,  which  diffuse  across  the 
width  Xb  base  and  are  then  all  swept  across  the  reverse-biased  collector  junction  by 
the  built-in  field  of  the  junction.  If  we  neglect  the  usually  small  leakage  current  of 
the  collector  junction  and  the  small  fraction  of  the  carriers  recombining  in  the  base, 
we  get  for  a  n'^pn  BJT 

Ic  *  AqDn[noBCxp(qVBE/kT)/XBl,  (1) 

where  A  is  the  cross  sectional  area  of  the  base,  q=-e  is  the  charge  of  the  minority 
carriers  in  the  base.  Dn  their  diffusion  coefficient  in  the  base, 

nB(0)=noBC  xp(qVBE/kT)  is  the  electron  concentration  at  the  limit  of  the  emitter  space 
charge  region,  Vbe  ‘s  the  applied  base  -  emitter  voltage,  and  Xb  is  the  width  of  the 
base.  The  v  xpression  in  rectangular  brackets  is  the  electron  concentration  gradient 
calculated  with  the  boundary  condition  of  a  vanishing  electron  concentration  at  the 
limit  of  the  collector  space  charge  region.  We  assume  the  base  to  be  much  narrower 
than  the  electron  diffusion  length  Ln=(DnX),  XB<<Ln.  but  sufficiently  wide  to  avoid 
ballistic  electron  transport  across  the  base.  Usually  Xb  is  a  fraction  of  a  micron. 

Quantum  l/f  fluctuations  of  the  coilisional  cross  sections  of  the  electrons  in 
the  base  will  yield  fluctuations  of  the  diffusion  constant,  and  of  the  mobility 

(6Dn/Dn=5ti/u) 

8Ic  =  Aq(8Dn)[noBCxp(qVBE/l^T)/XBl-  (2) 
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The  corresponding  spectral  density  of  fractional  fluctuations  I'^Sic  is 

IC-2<(8Ic)2>f  =  Dn-2<(SDn)2>  =  p-2<(5p)2>  =  On/fN.  ( 3 ) 

In  the  last  step  our  quantum  1/f  equation^'2  was  used,  where  N  is  the  number  of 
carriers  which  define  the  scattered,  or  diffused,  current  leaving  the  base  and 
emerging  in  the  collector,  while  On  =aAn  is  the  effective  quantum  1/f  noise 
coefficient,  or  Hooge  constant.  The  number  of  electrons  N  is  thus  determined  by  the 
effective  lifetime  t  of  the  electrons,  which  will  be  slightly  lower  than  the  lifetime  in 
the  unbounded  collector  material,  due  to  the  collector  lead  contact  processes,  and  due 
to  lateral  surface  recombination.  Indeed,  we  can  write  N  =  tIc/Q-  Thus  we  finally 
obtain  the  spectral  density  of  the  collector  current  fluctuations 

Sic  =  anIq/fT.  (4) 

in  which  t  ii'  the  effective  lifetime  of  the  majority  carriers  in  the  collector.  This 
expression  is  simpler,  but  similar  to  the  expression  derived  earlier,  with  the 
important  difference  that  now  we  have  a  lifetime  of  the  carriers  in  the  denominator, 
while  before  it  was  the  usually  much  smaller  diffusion  time  tD  =  XB2/Dn  of  the 
electrons  in  the  base.  Eq.  (4)  also  implies  that  in  narrow-base  BJTs  of  various  base- 
widths  On  will  be  constant,  as  in  other  devices,  rather  than  an/to- 


3.  A  sufficient  criterion  for  1/f  noise  in  chaotic  nonlinear  systems 

In  spite  of  the  practical  success  of  our  quantum  I/f  theory  in  explaining 

electronic  1/f  noise  in  most  high-tech  devices,  and  in  spite  of  the  conceptual  success 
of  our  earlier  classical  turbulence  approach  to  1/f  noise,  the  question  about  the 

ultimate  origin  of  nature's  omnipresent  1/f  spectra  remained  unanswered.  During 
the  last  three  decades,  we  have  claimed  repeatedly  that  nonlinearity  is  a  general 
cause  of  1/f  noise.  Our  new  result  proves  that  nonlinearity  always  leads  to  a  1/f 

spectrum  if  homogeneity  is  also  present  in  the  equation(s)  of  motion.  Specifically, 
let  ihe  system  be  described  in  terms  of  the  dimensionless  vector  function  Y(x,t)  by 
the  m'**  order  nonlinear  system  of  differential  equation 

<i>[t.  X,  Y.  aY/at,  aY/axi...aY/axn. 
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a2Y/at2.a2Y/axi2 . amY/axnm]=0  (5) 


where  the  vector  function  ^  may  be  nonlinear  in  any  of  its  arguments.  If  a  number 
9  exists  such  that  Eq.  (1)  implies 

<&[x.®t.  xx,  Y.  aY/xfiat,  aY/xaxi...aYAaxn.  b^y/xiW, 

a2YA2axi2 . amYA">axnm)  =  0  (6) 

for  any  real  number  X,  the  power  spectral  density  of  any  chaotic  solution  for  the 
vector  function  Y  defined  by  Eq.  (5)  is  proportional  to  1/f.  We  shall  present  here  a 
simplified  proof,  assuming  Eq.  (5)  to  contain  only  the  first  time  derivative. 

Consider  a  n-dimensional  nonlinear  system  described  in  terms  of  the 
dimensionless  vector  function  Y(x,t)  by  the  m^*'  order  nonlinear  dynamical  equation 

aY/at  +  F(x.  Y.  aY/axi...aY/axn.a2Y/axi2..  .amY/axH™)  =  o  (?) 

Here  F  is  a  nonlinear  vector  function  of  its  arguments  which  include  the  vector 
x(xi....xn).  If 

F[Xx.  Y,  aY/(Xaxi)...aY/(Xaxn).  a2Y/(Xaxi)2 . 9mY/(xaxn)m] 

=  x-PF(x. Y. aY/axi...aY/axn.  a2Y/axi2 . amY/axn"’)  (8) 

for  any  real  number  X,  Eq.  (7)  is  said  to  be  homogeneous  of  order  -p.  Performing  a 
Fourier  transformation  with  respect  to  the  vector  x,  we  get  in  terms  of  the  Fourier- 
transformed  wavevector  k  the  nonlinear  integro-differential  equation 

ay(k,t)/at  +  G  [k,  y(k,t).  kiy(k,t)...kjiy(k.t).  kj  2y(k,t)...kn‘^  y(k,t)]  =  0,  (9) 

where  y(k,t)  is  the  Fourier  transform  of  Y  (x,t).  Due  to  Eq.  (8),  the  nonlinear 
integro-differential  operator  G  satisfies  the  relation 

G[Xk.  y.  Xkiy...Xkny,  (Xki)2y . (Xkn)™y] 

=  XPG(k.y,  k|y...kny.  ki2y . kn^iyl.  (10) 

Eq.  (9)  can  thus  be  rewritten  in  the  form 
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dy/d(tAP)  +  G[Xk,  y,  X.k i y...Xkny.  (Xki)2y . (X.kn)'’^y]  =  0, 


(11) 


Taking  X=l/k,  where  k=lkl=(ki2+,,  .+^^^2)1/2^  setting  kPt=z,  we  notice  that  k  has 

been  eliminated  from  the  dynamical  equation,  and  only  k/k  is  left.  This  means  that 

there  is  no  privileged  scale  left  for  the  system  in  x  or  k  space,  other  than  the  scale 

defined  by  the  given  time  t,  and  expressed  by  the  dependence  on  z.  We  call  this 

property  of  the  dynamical  system  "sliding-scale  invariance" . 

In  certain  conditions,  instabilities  of  a  solution  of  Eq.  (7)  or  (5)  may  generate 
chaos,  or  turbulence.  In  a  sufficiently  large  system  described  by  the  local  dynamical 
equation  (7)  or  (S),  in  which  the  boundary  conditions  become  immaterial, 
homogeneous,  isotropic  turbulence,  (chaos)  can  be  obtained,  with  a  spectral  density 
determined  only  by  Eq.  (7)  or  (5).  The  stationary  autocorrelation  function  A(x)  is 
defined  as  an  average  scalar  product,  the  average  being  over  the  turbulent  ensemble 

A(x)  =  <Y(x,t)Y(x.t+t)>  =  |<y(k,t)y(k.t+T)>d"k  =  Ju(k,z)d"k.  (12) 

Here  we  have  introduced  the  scalar 

u(k,z)  =  <y(k,t)y(k,t+t)>  (13) 

of  homogeneous,  isotropic  chaos  (turbulence),  which  depends  only  on  Ik  I  and  z=kPT. 

All  integrals  are  from  minus  infinity  to  plus  infinity.  The  chain  of  integro- 
differential  equations  for  the  correlation  functions  of  any  order  obeys  the  same 

sliding-scale  invariance  which  we  have  noticed  in  the  fundamental  dynamical 
equation  above.  Therefore,  in  isotropic,  homogeneous,  conditions,  u  can  only  depend 

on  k  and  z.  Furthermore,  the  direct  dependence  on  k  must  reflect  this  sliding-scale 

invariance,  and  is  therefore  of  the  form 

u(k.z)  =  k-"v(z).  (14) 

Indeed,  only  this  form  insures  that  u(k,z)d"k  and  therefore  also  the  corresponding 
integrals  and  multiple  convolutions  in  k  space  have  the  necessary  sliding-scale 
invariance. 
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According  to  the  Wiener-Khintchine  theorem,  the  spectral  density  is  the 
Fourier-transfom  of  A(t), 

Sy(f)  =  Je2nif^A(T)dr  =  (1/f)  Je2««t’jk'-nv(z)d"k'df  =  C/f.  (15) 

where  we  have  set  fx=t',  k''=:fk'",  z=k"t=k'"t’,  and  the  integral 

C  =  I  c2nit'J  k'-nv(z)d"k'dt’  =  J  e2Jtit7  k"-nv(k"")d"k"dt’  (16) 

is  independent  of  f.  We  have  defined  the  vector  k''=t'^/”  k. 

Eq.  (15)  proves  the  criterion  as  stated  in  Eq.  (8)  The  extension  to  the  form  (5-6) 
is  trivial.  In  conclusion,  nonlinearity  +  homogeneity  =  1/f  noise  if  the  laminar 

(nonchaotic)  solution  is  unstable.  The  ultimate  cause  of  the  ubiquitous  1/f  noise  in 
nature  is  the  omnipresence  of  nonlinearities  (no  matter  how  we  and 
homogeneity.  The  latter  is  finally  related  to  rotational  invariance  and  to  the  isotropy 
of  space.  All  our  four  specific  theories  of  1/f  chaos  in  nonlinear  systems  are  just 
special  cases  to  which  our  criterion  is  applicable.  They  include  our  magneto-plasma 
theory  of  turbulence  in  intrinsic  symmetric  semiconductors  (1966),  our  similar 
theory  for  metals  (1971),  the  quantum  1/f  theory  (pure  quantum  electrodynamics, 

1975),  and  the  theory  of  Musha's  traffic  turbulence  (1989).  A  fifth  application 

concerns  a  one-dimensior.al  crystal,  i.e.,  a  chain  of  atoms  with  slightly  anharmonic 

interaction  potentials,  which  is  presented  next. 


4.  Application  of  the  sufficient  criterion  to  a  chain  of  atoms  with  a  harmonic 
coupling 

Consider  a  chain  of  atoms  in  the  x  direction,  with  a  lattice  constant  b  and 
displacements  qi  from  the  equilibrium  position.  The  equations  of  motion 

md2qn/dt2=  A((qn+ 1 -qn)-(qn-qn-l  )1  +  B[(qn+i -qn)^-(qn-qn-l  )^] 

+  C[(qn+i-qn)^-(qn-qn-l)^l  (17) 

contain  anharmonic  terms  as  long  as  B  and  C  are  different  from  zero.  With 
On=qn+l*qn  Pn=qn'qn-1-  neglecting  the  second-order  term,  we  obtain 
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md2qn/dt2=  [(qn+l-qn)-(qii-qn-l)][A  +  B(qn+l+qn-l)  +  CCan^+anPn+Pn^)]  (18) 

Going  over  from  finite  differences  to  a  continuum  description^  we  obtain  a 
differential  equation  for  q(x,t) 

m32q/3t2  =  b^3^q/3x2[A  +  2Bq  +  3Cb2(3q/3x)^].  (19) 

Performing  a  fourier  transform  with  respect  to  x, 
m32qk/3t^  =  -Ab^k^q^  -  2h^B  Jk'^qk'qic-k'dk' 

+  3Cb4/dk’/dk"k’k"(k-k’-k”)2qk’qk"qk-k'-k".  (20) 

All  integrals  are  from  minus  infinity  to  infinity.  Substituting 

qk=u(k.t)cxp[ikbt(A/m)l/2], 

m32u/3t2  +  2ikb(A/m)l/23u/3t  =  -2b2B  Jk'^uk’Uk-k'dlt’ 

+  3Cb4Jdk’/dk"k'k"(k-k’-k'*)2uk'Uk'’Uk-k’-k".  (2 1 ) 

Our  1/f  noise  criterion  requires  both  nonlinearity  and  homogeneity,  as  well  as  the 

presence  of  chaos  or  of  a  quasichaotic  state.  The  nonlinearity  condition  is  satisfied 
unless  B=0  and  C=0,  while  homogeneity  requires  the  existence  of  two  numbers  p  and  9 
such  that  replacing  k  by  A.k  everywhere  except  in  the  integration  differentials,  and 
replacing  t  by  .1^1  leaves  the  equation  multiplied  by  a  general  factor  XP,  i.e.,  formally 
invariant.  In  our  case  we  note  that  this  criterion  is  satisfied  with  p=2  and  9=-l  if  we 
neglect  the  third-order  term  bv  setting  C=Q.  On  the  other  hand,  both  in  the  general 
case  and  in  the  B=Q  case  the  criterion  is  not  satisfied,  except  for  some  low-frequency 
Iimitiilg...case..  in  which  all  k  values  and  frequencies  are  so  small,  that  we  can  neglect 
the  term  with  B  and  one  of  the  left  hand  side  terms. 

To  see  how  the  criterion  works  for  C^Q.  we  set  X.=l/lkl  and  call  lklt=z 

m32u/3z2  +  2ib(A/m)l/23u/3z  =  -2b2B  J(k7k)2uk'Uk-k'dk’  (22) 


Substituting  u(k,t)=k*  *  v(k,t)  we  get  for  v 
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(23) 


nid2v/az2  +  2ib(A/m)l/2av/az  =  -2b2B/(k7k)v(k'.z)v(k-k’.z)dk7(k-k') 

=  -2b2B  /k"  v(kk".z)v(k-kk"  ,z)dk7(  1  -k") 

We  note  that  k  has  disappeared  from  the  equation  and  is  present  only  as  a  scale  factor 
in  the  arguments  of  v  on  the  rhs.  Therefore  we  can  expect  the  existence  of  solutions 
v(k,z)  of  this  equation  which  do  not  depend  on  the  first  argument.  Such  solutions 

exhibit  "sliding  scale  invariance",  because  t  and  k  or  t  and  x  provide  a  scale  for  each 

other,  with  no  other  scale  present. 

In  certain  conditions,  instabilities  of  a  solution  of  Eq.  (20)  may  generate  chaos, 
or  turbulence.  In  a  sufficiently  large  system  described  by  the  local  dynamical 

equation  (6),  in  which  the  boundary  conditions  become  immaterial,  homogeneous, 
isotropic  turbulence,  (chaos)  can  be  obtained,  with  a  spectral  density  determined 
only  by  Eq.  (20).  In  the  absence  of  instability  and  chaos,  a  certain  type  of  random 

stirring  forces  can  generate  a  quasichaotic  stochastic  state  which  can  also  be 

described  with  our  methods  familiar  from  turbulence  theory.  The  stationary 
autocorrelation  function  A(t)  is  defined  as  an  average  scalar  product,  the  average 
being  over  the  turbulent  ensemble 

A(t)  =  <u(x,t)u(x,t+t)>  =  J<uic(t)uic(t+T)>dk  =  /U(k,z)dk.  (24) 

Here  we  have  introduced  the  scalar 

U(k,z)  =  <uk(t)uk(t+x)>  (25) 

of  homogeneous,  isotropic  chaos  (turbulence),  which  depends  only  on  Ikl  and  z=lklPt, 

because  there  is  nothing  else  in  Eqs.  (20)  and  (21).  All  integrals  are  from  minus 

infinity  to  plus  infinity.  The  chain  of  integro-differential  equations  for  the 

correlation  functions  of  any  order  obeys  the  same  sliding-scale  invariance  which  we 
have  noticed  in  the  fundamental  dynamical  equation  (21)  above.  Therefore,  in 
isotropic,  homogeneous,  conditions,  Uk(t)  can  only  depend  on  k  and  z.  Furthermore, 
the  direct  dependence  on  k  must  reflect  this  sliding-scale  invariance,  and  is 
therefore  of  the  form 

Uk(z)  =  Ikl-lv(z).  (26) 
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Indeed,  only  this  form  insures  that  uic(z)dk  and  therefore  also  the  corresponding 
Integrals  and  multiple  convolutions  in  k  space  have  the  necessary  sliding-scale 
invariance. 

According  to  the  Wiener-Khintchine  theorem,  the  spectral  density  is  the 
Fourier-transform  of  A(t), 

Su(f)  =  Je2«if^A(t)dT  =  (1/0  Jdk'Jdt'e2nit’k-lv(z)  =  C/f,  (27) 

where  we  have  set  fT=t',  k=fk',  z=kT=k't',  and  the  integral 

C  =  Jdk’/dt'e2nit'k’-lv(z)  =  |dk"/dt‘e2itit'k-lv(k")  (28) 

is  independent  of  f.  We  have  defined  the  vector  k"=t’  k . 

The  1/f  spectrum  obtained  by  us  for  the  amplitude  u(t)  carries  over  also  for 
the  squared  amplitude  u2(t).  Indeed,  the  autocorrelation  A’('t)=<u2(x,t)u2(x,t+T)>  is 
given  by  2A2(-c)-h a2(0)  if  we  assume  the  amplitude  u(t)  to  be  well  approximated  by  a 
Gaussian  process.  The  Fourier  transform  of  a2(t)  is  the  autoconvolution  of  C/f, 
which  is  C2/f,  if  we  interpret  all  1/f  spectra  as  the  limit  of  spectra  for  arbitrarily 

small  e.  Therefore,  Su2(0  =  2C2/f  +  a2(0)5(0.  where  5(0  is  the  delta  function.  The 
energy  density  and  the  phonon  number  density  are  both  proportional  to  u2.  This 
proves  that  in  this  case  the  energy  density  and  the  phonon  number  density  are  both 
fluctuating  with  a  1/f  spectral  density  if  they  fluctuate  at  all,  i.e.,  if  the  system  is 
either  chaotic,  or  in  a  quasichaotic  state  caused  by  a  suitable  system  of  random 
stirring  forces.  This  confirms  for  the  one-dimensional  case  the  prediction  of  T. 
Musha.  In  three-dimensional  piezoelectric  crystals  similar  I/f  fluctuations  of  the 
phonon  number  are  predicted  by  the  quantum  1/f  theory  as  we  show  next,  and  have 
been  observed  experimentally  in  the  Brillouin  scattering  of  light  by  T.  Musha^^. 


5.  Quantum  1/f  effect  in  quart  resonators 

According  to  the  general  quantum  1/f  formula,  p( f)  =  2a A/f  with  a  =  e2/hc 

=  1/137  and  A=2(Aj/ec)2/3jt  is  the  quantum  1/f  effect  in  any  physical  process  rate  F. 
Setting  j  =  dP/dt=P,  where  P  is  the  vector  of  the  dipole  moment  of  the  quartz  crystal, 
we  obtain  for  the  rate  F  of  phonon  removal  from  the  mam  resonator  oscillation  mode 
by  scattering  on  a  phonon  from  any  other  mode  of  the  crystal  the  spectral  density 
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Sr(f)  =  r24a(AP)2/37rc2c2. 


(29) 


where  (aP)2  is  the  square  of  the  polarization  rate  change  associated  with  the 
removal  of  one  of  the  N  phonons  present  in  the  main  resonator  mode.  To  calculate  it, 
we  write  the  energy  W  of  the  resonator  and  its  change  in  the  form 

W  =  (N+l/2)hto  =  p2/2Vx(o2;  AW  =  ho>  =PAP/Vxo)2.  (30) 


Her  X  is  the  susceptibility  and  V  the  volume  of  the  quartz  crystal.  Solving  the 

last  equation  for  AP,  squaring,  and  multiplying  with  the  first,  we  get 

(AP)2  =  h(i)3Vx/(2N+l),  (31) 


and 


Sr(0  =  r24ali(oVx/3jie2c2(2N+I)f  =  r2(2tD3x/3Jtc3f)(h(DV/W)  (32) 


The  corresponding  frequency  fluctuations  arc  given  by 

a)-2Sa,(0  =  (2o)3x/37tc3f)(h(oV/W).  (33) 


6.  Quantum  1/f  mobility  fluctuations  in  semiconductors 

A  first  principles  calculation  of  quantum  1/f  cross  correlations  performed^ 
for  the  first  time  in  1987  has  yielded  a  slightly  different  result  compared  to  earlier 
expectations.  This  same  new  form  of  the  quantum  1/f  cross  correlations  was  derived 
again  with  a  different  method  by  Van  VUct  in  1989.  It  differs  from  the  old  form  used 
in  the  1985  calculation  of  Kousik  et.  al.  by  a  correction  which  is  zero  when  the 
momentum  changes  of  the  two  current  carriers  involved  in  the  cross  correlation  are 
identical,  but  increases  to  finite  values  when  the  momentum  differences  caused  by 
the  scattering  process  are  different.  The  correction  is  proportional  to  the  squared 
difference  of  the  two  momentum  changes.  We  have  repeated  all  calculations  in  the 
original  paper  by  Kousik  et.al'®.,  obtaining  both  for  impurity  scattering  and  for  the 
various  types  of  phonon  scattering  new  analytical  expressions  which  show  a 
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considerable  increase  of  the  final  quantum  1/f  noise.  The  results  obtained  are 
applicable  both  to  direct  and  indirect  bandgap  semiconductors. 

6.1  Introduction 

We  have  performed  an  analytical  calculation  of  mobility  fluctuations  in 

silicon  and  gallium  arsenide,  using  the  new  quantum  1/f  cross-correlations  formula. 
This  calculation  is  of  major  importance  for  the  l/f  noise-related  optimization  of  the 
two  types  of  materials,  and  of  the  many  devices  constructed  with  them  for  military 
and  civilian  applications  in  the  electronic  and  opto-electronic  industry. 

The  new  cross-correlation  formula  gives  the  cross-spectral  density  which 

describes  the  way  in  which  simultaneous  quantum  1/f  scattering  rate  fluctuations 
AW  observed  in  the  direction  of  the  outgoing  scattered  wave-vector  K’  arc  correlated 

with  those  in  the  K  "  direction,  when  the  two  corresponding  incoming  current 
carriers  have  the  wave  vectors  Ki  and  K2: 

Saw(Ki,K':K2,K";0 

=  (2a/3jtf)(h/m*c)2WKi.K'WK2,K"((K'-Ki)2+(K"-K2)2l5Ki.K2-  (34) 

The  form  conjectured  by  us  earlier  had  2(K ’-K 1  ){K  "-K  2)  in  place  of  the  rectangular 

bracket. 

6.2  Impurity  Scattering 

For  impurity  scattering  of  electrons  in  solids,  fluctuations  At  of  the  collision 
times  T  will  cause  mobility  fluctuations 

Atiband(0  =  [c/m*«v2»]lKVK2Ax(t)nK-  (3  5) 

where  «v^»  is  both  the  average  over  all  states  of  wave-vectors  K ,  with  occupation 
numbers  ok-  in  the  conduction  band,  and  the  thermal  equilibrium  average  of  the 
quadratic  carrier  velocities.  With  the  help  of  the  relation 

1/t(K)  =  (V/8«3)/(1  -  cos07cose)WK,Kd^K’.  (36) 

the  mobility  fluctuations  arc  reduced  to  fluctuations  of  the  elementary  scattering 
rates  Wk.k  -  governed  by  Eq.  (34).  Here  V  is  the  volume  of  the  normalization  box 
which  disappears  in  the  final  result,  0  and  0’  respectively  the  angles  K  and  K'  form 


with  the  direction  of  the  applied  field.  One  finally  obtains  after  tedious  multiple 
integrations 

^■2SA^(f)  =  [256rtaK2e4hl2/3m*8z4e8Ni2l(l/f)ZKKl0[ln(l+a2)-a2/(l+a2)]-3 
[(2a2+a4)/(l+a2)-21n(l+a2)]F(EK)(2KVK2T(K)F(EK)r2,  (3  7) 

where  a=2K/K,  K2=e2n(T)/ekBT,  n(T)  is  the  electron  concentration,  F(EK)=exp(EF-EK) 
for  non-degenerate  semiconductors.  Nj  the  concentration  of  impurities  of  charge  Ze 
and  E  the  dielectric  constant.  The  corresponding  partial  Hooge  parameter  for 
impurity  scattering  is  thus 


Oj  =  [4V2rcaKfi3N^/3tn*2/2(kgX)3/2(;2]  1  l/2g-x 

0 

[In(bx+l)-bx/(bx+l)]-3[(2bx+b2x2)/(bx  +  l)-21n(bx+l)] 

oe 

{  Jdxx3e-*[ln(bx+l)-bx/(bx+l)]*l}-2.  (38) 

0 

6.3  Acoustic  Electron-Phonon  Scattering 

In  this  case  the  calculation  is  similar,  and  leads  to  the  result 

oo 

Oac  =  [32noNcm*C2h3/3c2kBT)^]{(l/R2)  jdxx--^ 

1 

[(x-1)2/7+(R+1)(x-1)6/6+R(x- 1)5/5] 

((x-l)5/5  +  (R+l)(x-l)'*/4+R(x-I)3/3]exp(-x2/4R) 

I 

+  Jdxx-'‘[(x+l)5/5-(x+l)6/6+(x-l)5/5+(x- 1)6/6] 

0 

[(x+l)3/3+(x-l)^/4-t-(x-l)3/3-(x+l)4/4]exp(-x2/4R) 

OO 

+  /dxx-4((x+l)5/5-(x+l)6/6][(x+l)3/3-(x+l)V4]exp(-x2/4R)}.  (39) 

1 

where  R=kBT/2m*C  1 2,  Ci  is  the  deformation  potential,  and  Nc  is  the  effective  density 
of  states  for  the  conduction  band. 

6.4  Non-Polar  Optical  Phonon  Scattering 
This  time  one  obtains 
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oo 


On  o  ph  =  [SW 2  h{OoaNcli^/3m*5/2c2{^^]^  |dxx5/2 

0 

[(F+l)(x-l)l/2e(x-l)+F(x+l)l/2l-4 

[(F+1  )2(x-l  )(2x-l)0(x-l)+F2(x+l)(2x  +  l)]exp(-h(Oox/kBT)} 

{  |dxx3/2[(F+l)(x-l)l/2e(x-l)+F(x+l)l/2]-lexp(-ho)ox/kBT)}-2,  (40) 

0 

where  F=[exp(h(i)o/kBT)-l]'*.  and  (Oq  is  the  optical  phonon  frequency. 

6.5  Polar  Optical  Phonon  Scattering 
Proceeding  as  in  Secs.  2  and  4,  we  obtain 

_  oe 

“p  o  oh  =  [8jtV2  ha)|oNcfi2/3m*5/2(;2(^j]{  /dxx^ 

0 

[F2(x+l)l/2ln(2xl/2+2(x+ 1)1/2) 

+  (F+l)2(x-l)l/2ln(2xl/2+(x-l)l/2)0(x-l)lexp(-ha)|x/kBT) 
[(F+l)arcsinh(x-l)l/20(x-l)+Farcsinh(xl/2)]-'*}.  (41) 

Here  coi  is  the  longitudinal  phonon  frequency. 

6.6  Intervalley  Scattering 

This  type  of  scattering,  present  in  indirect  bandgap  semiconductors,  transfers 
electrons  from  one  of  the  six  minima  (or  valleys)  of  the  conduction  band  energy  in 
k-space  to  one  of  the  other  five  minima.  Transitions  between  a  valley  and  the 
nearest  valley,  which  is  along  the  same  k-space  direction  in  the  next  copy  of  the  first 
Brilloin  zone  in  the  periodic  zone  scheme,  are  of  the  Umklapp  type,  and  are  called  g- 
processcs.  Transitions  to  the  four  valleys  present  in  the  same  zone  along  the  other 
two  k-space  directions  arc  called  f-proccsscs.  Repeating  a  previous  calculation! 0  on 
the  basis  of  the  new  cross-correlation  formula^’  H  (34),  we  obtain  for  g-processes 

_  OO 

Og  =  (8W2  ha)jjaNcfi2/3ni*5/2c2(jjj]^  Jdxx5/2 

((F+l)(x-l)l/20(x-l)+F(x+l)l/2]-4 

[(F+1  )2(x-1)(2x-1)0(x-1)+f2(x+1)(2x+1  )Iexp(-ho)jjx/kBT)} 
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{  Jdxx3/2[(F+l)(x-l)l/29(x-l)+F(x+l)l/2]-lexp(-ha)iix/kBT)}-2,  (42) 

0 

where  hojij  is  the  phonon  energy  corresponding  to  the  momentum  difference 
required  by  the  intervalley  transition.  For  the  corresponding  f-process  we  obtain^  ^ 

«f  =  (ko/qo)^Og.  (43) 

where  ko/qo  is  the  ratio  between  the  position  vector  of  a  conduction  band  energy 
minimum  in  k  space,  and  twice  the  distance  of  the  minimum  from  the  Brillouin  zone 
boundary,  0.85/0.3  for  silicon.  There  are  three  g-typc  (from  LA,  TA  and  LO  phonons) 
and  three  f-type  contributions  (from  TA,  LA  and  TO  phonons). 

The  various  quantum  1/f  contributions  derived  here  can  be  approximately 

superposed  to  yield  the  resultant  quantum  1/f  coefficient  according  to  the  rule 

ttH  =  2i(n/4i)2  a,  (44) 

The  support  of  the  Army  Research  Office  is  thankfully  acknnowledged. 
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APPENDIX  D 


EXAMPLES 

The  general  criterion  developed  in  the  main  text  will  now  be  illustrated  on  the  basis 
of  more  examples. 

D.l  Classical  Turbulence  Theory  for  the  Current  Carriers  in  Semiconductors 

In  the  case  of  homogeneous,  isotropic  turbulence^ caused  in  the  electron-hole 
plasma  of  an  infinite  sample  of  a  symmetric  intrinsic  semiconductor  by  dynamical 
instabilities  of  any  kind,  we  start  from  the  equations 


vv*"  =  (e/2c)vxB  -  (l/n)VP, 

(Dl) 

vv  =  2e[E  +  v+xB/c)  -  (2/n)V(Pp  -  Pn). 

(D2) 

V  v+a  0  (n  =  const). 

(D3) 

V  xE  =  -(I/c)aB/at. 

(D4) 

V  xB  =  2iienv7c, 

(D5) 

VB  =  0. 

(D6) 

Here  n  is  the  total  carrier  concentration  including  an  equal  number  of  electrons 
and  holes,  v/e  their  reciprocal  mobility  assumed  to  be  the  same  for  electrons  and  holes,  Pn 
and  Pp  the  partial  pressures  of  electrons  and  holes,  P  the  total  carrier  pressure.  2V''’  v  the 
sum  and  the  difference  of  the  carrier  drift  velocities.  Inertial  terms  proportional  to  the 
effective  masses  of  the  carriers,  as  well  as  electrostatic  terms  and  compressibility  terms 
have  been  neglected  here  in  a  consistent  way'*’^,  because  we  are  interested  in  the  low- 
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frequency  domain  only.  Although  we  do  not  work  this  out  here,  this  system  of  equations 
can  be  shown  to  admit  an  energy  theorem.  Performing  a  Fourier  expansion,  we  obtain 


vv  +  (k)  =  (e/2c)Zk'V(k')xB(k-k')  -  (i/n)k  P(k),  (D7) 

vv-(k)  =2e[E(k)+  Zkv  +  (k’)xB(k-k')/c]-(2i/n)k(Pp-Pn),  (D8) 

k  v-*-(k)  =  0.  (D9) 

ik  xE(k)  =  -(l/c)8B(k)/3t.  (DIO) 

ik  xB(k)  =  2itenv(k)/c,  (Dll) 

k  B(k)  =  0.  (D12) 


Substituting  E  from  Eq.  (D8)  into  Eq.  (DIO),  we  obtain  with  Eq.  (Dll) 

aB(k)/at  +  qk2B(k)  =  ikx  Ik'  v+(k’)  xB(k-k’).  (D13) 

where  u  =  c^v/4itnc^.  Eqs.  (D7)  and  (Dll)  yield 

v  +  (k)  =  (i/47tvn)  Ik"{  B(k")[k"  B(k-k"))  -  k"[B(k")  B(k-k'’)) 

(D14) 

-(k/k2)(i-5k,o)(l‘"  B(k-k")  k  B(k")  -  B(k")  B(k-k")  (k  k")] }  . 

Substituting  this  into  Eq.  (D13),  we  obtain  the  fundamental  dynamical  field-equation  of 
turbulence  in  the  electron-hole  plasma  of  a  symmetric  intrinsic  semiconductor 

abp(k,t)/at  +  pk2bp(k.t)  =  Ik  k"  bj(k-k',t)  bi(k",t)  bm(k'-k",t) 

(D15) 

(kj5ps-ks5pj)[ks5ijn-kjn  Sjj  +  (kj /k  2)( l -5k',o)(^m  kj-k  k  5im)], 
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in  terms  of  b  s  B/(2rtvn)^/^.  This  dynamical  equation  has  the  form  of  Eq.  (3),  with  p=2  m 
Eq.  (4),  and  with  G  defined  as  the  r.h.s.  minus  the  term  in  on  the  l.h.s..  Our  sufficient 
criterion  thus  tells  us  that  this  nonlinear  system  will  yield  a  1/f  spectrum.  We  present 
below  the  direct  derivation  for  this  example. 

Multiplying  Eq.  (D15)  with  ba*(k,t-T)  and  taking  the  average  over  a  statistical 

ensemble  which  represents  our  notion  of  stationary  turbulence,  we  obtain  in  quasi¬ 
stationary  conditions 

3wap(k,x)/3T  +  pk^wapfk.x) 

=Ik'k"<ba*(k.t-x)bj(k  k’.t)bi(k",t)bm(k'-k",t)>Rjmip.  (D16) 

with  Wa^(k,x)  s  (L/2jt)3<ba*(k,t-x)bp(k,t)>,  where  L  is  the  edge  of  the  cubic  normalization 
box,  and 


Rjmip  ^(^j^ps'ks5pj)[ks5im-km  Sij  +(ks’/k'^)(l-5|j'_o)(kni"ki’-k' k"5]in)]'  (D17) 

Multiplying  Eq.  (D15)  with  more  magnetic  field  components  and  averaging,  we 

obtain  equations  connecting  the  fourth-order  correlation  tensor  to  the  sixiSi-order  tenso-, 
and  so  on‘^‘7.  To  end  this  infinite  chain  of  equations  for  the  correlation  functions,  w>: 
make  a  quasinormality  assumption  which  expresses  the  fourth-order  moment  appearing  in 

Eq.  (D16)  according  to  the  scheme 

<ABCEt>=<AB><a3>+<AC><BDC>+<AEt><BC>,  (D18) 

valid  if  the  four  field  components  would  approximate  a  joint  normal  distribution.  This 

approximation  does  not  alter  the  homogeneity  of  the  system,  which  ultimately  causes  the 

1/f  spectrum.  This  approximation  yields  the  closed  equation 

^Wa3(k,r)/3x  +  pk2wa0(k,x) 

=  (27t/L)3wotj(k,x)Zk'W|n,(k’,0)Rjnil3(k,k'),  (D19) 
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Isotropic  turbulence  requires  Wap  =  Ai(k)5ap  +  A2(k)kakp,  with  coefficients  Ai  and  A2 
related  through  Eq.  (D12),  yielding 


Wap(k.'t)=  (l/2)[5ap  -  kakp/k2]u(k.t). 


(D20) 


where  u(k,T)  =  laWaa(kit).  Therefore,  the  scalar  correlation  function  u(k,x)  satisfies  the 
dynamical  equation  of  homogeneous,  isotropic,  stationary  turbulence 


dixl 


+  v(k,x)  =  -  jv(k 


,x)J- 


d^k’  k^+k  k'  k  k'  - 


(D21) 


where  v(k,x)s3k'3u(k,T),  and  x  =  ptk^  is  a  dimensionless  variable  replacing  x.  We  convince 
ourselves  that  the  integral  is  independent  of  k,  provided  v(k,x)  does  not  depend  on  its  first 
argument,  by  setting  k’/ksK.  This  yields  a  solution.  However,  with  v=e'"''’^'  we  get  a 
logarithmic  divergence  at  k®0.  We  look  for  an  exact  solution  of  the  form'*'^ 

v(k,x)  =  hkEe-'’^'"’^''),  or  u(k,x)  =  (h/k^  £)  (D22) 


where  m(k)  is  very  close  to  a  constant,  almost  independent  of  k,  and  h  a  constant 
proportional  to  the  intensity  of  the  turbulence,  or  the  turbulence  level.  Substituting  this 
into  Eq.  (Dll)  and  performing  the  integration,  we  obtain  a  finite  result  only  for  0<e<2; 

m(k)  =  I  +  h  r(e)ke,  with  r(E)=[27c2cotan(EJt/2)]/[{l-e2)(3-E2)].  (D23) 

We  notice  that  m(k)  is  indeed  practically  constant  when  0<e«1  is  very  small,  arbitrarily 
small.  The  value  e=0  leads  to  a  loganinmic  divergence,  but  we  can  set  e=0  for  practical 
purposes. 
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which  is  proportional  to  interesting  to  note  that  for  e«l  cotan  e)X/2  »  llt%  in 

Eq.  (23),  and  that  the  value  of  e  calculated  from  Eq.  (23)  is  therefore  proportional  to  h,  or  to 
the  intensity  of  the  turbulence.  This  feature  of  the  classical  theory*'*’  is  expressed  with 
fascinating  clarity  in  the  quantum  form  of  the  theory,  where  e  is  replaced  by  2a A  which 
also  appears  as  an  intensity  factor  multiplying  the  quantum  1/f  noise. 

The  essential  element  which  led  to  the  1/f  spectrum  in  the  classical  turbulence 
theory  is  the  nonlinearity  of  the  equations  of  motion,  caused  by  the  reaction  of  the  electric 
currents  back  on  themselves  via  the  generated  electromagnetic  field.  The  same  feedback 
reaction,  via  the  electromagnetic  field,  also  caused  the  nonlinearity  in  the  quantum  1/f 
theory,  and  in  QED  in  general,  leading  in  the  same  way  to  an  identical  1/f  spectrum,  this 
time  with  a  physically  more  meaningful  E=2aA.  This  nonlinearity  induces  the  coupling 
between  various  scales  of  turbulence  and  leads  to  the  dynamical  equilibrium  between 
eddies  of  all  sizes,  expressed  by  the  1/f  spectrum.  In  the  e=0,  or  aA=0,  limit,  this  dynamical 
equilibrium  assumes  both  for  the  quantum  case  and  for  homogeneous,  isotropic, 
turbulence  in  the  unbounded  semiconductor  sample  the  simplest  form,  characterized  by 
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scale-homogeneity,  or  scale  invariance.  Indeed,  replacing  for  e=0  in  Eqs.  (D21)  and  (D22)  k 
and  k’  by  Xk  and  Xk',  while  also  replacing  t  by  t/X^,  (or  to  by  X^to),  Eq.  (D21)  is  not  affected, 
and  X  drops  out.  We  conclude  that  in  the  weak  turbulence  limit  (e=0)  we  obtain  perfect 


at  all  scales  in  space  and  time,  classically  and 


quantum-mechanically .  The  implied 


is  caused  by  the  absence  of  any 


characteristic  length  or  time  scale,  or  by  the  presence  of  a  sliding  scale.  Indeed,  the 
frequency  scale  pk^  is  a  function  of  the  size  of  the  eddies,  given  by  the  wave  number  k 
which  can  have  any  value.  The  actual  frequency  and  wave-number  spectra  are  closely 


L,  but  in  the  weak-turbulence  limit  they  approach  an  exact  1/f  and  IDs} 


spectrum  respectively.  In  fact,  we  are  here  understanding  the  nonlinear  dynamics  which 
shapes  this  fractal  for  the  first  time. 


D.2  Turbulence  Theory  for  Orude  Electrons  in  Metals 


Our  classical  turbulence  theory  can  be  extended  to  the  case  of  metals  or  degenerate 
extrinsic  semiconductors^  in  the  Drude  model.  The  system  of  integro-differential 
equations  is  quite  different. 


vv-*-  =  -eE  -  (e/c)vxB  -  (l/n)VP 

(D26) 

VxE  =  -(l/c)aB/at 

(D27) 

VxB  =  -(4xen/c)v 

(D28) 

V.B  =  0, 

(D29) 

and  leads  to  a  third-order  nonlinearity^  in  the  resulting  closed  equation  of  motion,  or 
nonlinear  field-equation,  which  replaces  Eq.  (D14): 

9B(k,t)/3t  +  vk^Bfk.t)  =  -(c/47ine)kxj  d3k'B(k-k',t)x[k'xB(k',t)].  (D30) 
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This  is  again  in  the  form  of  Eq.  (3),  with  p=2  in  Eq.  (4).  We  thus  expect  a  1/f  spectrum  in 
this  system  as  well.  This  time  we  only  sketch  the  derivation.  The  corresponding  infinite 
chain  of  equations  for  the  correlation  tensors  now  goes  in  steps  of  one.  As  was  shown 
above,  it  went  in  steps  of  two  for  semiconductors.  The  third-order  correlation  can  be 
eliminated  between  the  first  and  second  equations  in  the  chain.  The  resulting  dynamical 
equation^  for  homogeneous,  isotropic,  stationary  turbulence,  which  replaces  Eq.  (D21), 
with  the  same  notations,  using  63  as  the  unit  vector  of  the  third  axis,  is 


3^v(k,x) 

0x2 


+ 


v(k.x) 

1  r  d^K 

<3 


Ie3+K|3l‘‘  J 


v(Kk,K2x)v(kle3+Kl,xle3+ic|2). 


(D31) 


This  also  admits,  in  the  e=0  limit  of  weak  turbulence,  a  solution  v(k,x)  which  does  not 
depend  on  the  first  argument,  and  u(k,x)  *  j(-3e-xm(x)^  t^jjs  time  with  an  x-dependent  m. 
With  the  change  of  variables  t=(i)T  and  k'=k/-Va)  in  the  second  .'middle,  invol'-ing  u)  form  of 
Eq.  (D24),  X  remains  invariant,  and  a  factor  l/o)  will  appear  in  front  of  the  integrals  which 
themselves  will  just  yield  a  constant  factor  independent  of  to,  t  or  k.  We  thus  obtain  again  a 
universal  1/f  spectrum.  As  is  shown  in  detail  elsewhere^,  this  1/f  spectrum  is  expressed  in 
the  corresponding  current  and  voltage  fluctuations  which  can  be  observed  in  the 
semiconductor  or  metallic  medium.  We  conclude  that  the  1/f  spectrum  is  a  general 
property  of  electrically  conducting  systems  in  interaction  with  the  electromagnetic  field,  a 
property  which  is  caused  by  the  nonlinearity  of  the  system  of  carriers  and  field  in  mutual 
interaction  due  to  the  absence  of  a  characteristic  scale  in  the  nonlinear  equation  of 
motion,  and  which  finds  its  clearest  expression  in  the  Quantum  1/f  Effect. 
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